home *** CD-ROM | disk | FTP | other *** search
open in:
MacOS 8.1
|
Win98
|
DOS
browse contents |
view JSON data
|
view as text
This file was processed as: LaTeX Document
(document/latex).
Confidence | Program | Detection | Match Type | Support
|
---|
100%
| dexvert
| LaTeX Document (document/latex)
| magic
| Supported |
1%
| dexvert
| Text File (text/txt)
| fallback
| Supported |
100%
| file
| LaTeX document text
| default
| |
99%
| file
| LaTeX document, ASCII text
| default
| |
100%
| checkBytes
| Printable ASCII
| default
| |
100%
| perlTextCheck
| Likely Text (Perl)
| default
| |
100%
| detectItEasy
| Format: plain text[LF]
| default (weak)
|
|
hex view+--------+-------------------------+-------------------------+--------+--------+
|00000000| 5c 64 6f 63 75 6d 65 6e | 74 73 74 79 6c 65 5b 74 |\documen|tstyle[t|
|00000010| 77 6f 73 69 64 65 5d 7b | 73 69 61 6d 6c 74 65 78 |woside]{|siamltex|
|00000020| 7d 0a 0a 25 20 64 65 66 | 69 6e 69 74 69 6f 6e 20 |}..% def|inition |
|00000030| 6f 66 20 7b 64 65 66 69 | 6e 69 74 69 6f 6e 7d 20 |of {defi|nition} |
|00000040| 74 68 65 6f 72 65 6d 2d | 6c 69 6b 65 20 65 6e 76 |theorem-|like env|
|00000050| 69 72 6f 6e 65 6d 65 6e | 74 2c 0a 25 20 77 68 69 |ironemen|t,.% whi|
|00000060| 63 68 20 69 73 20 6e 6f | 74 20 69 6e 63 6c 75 64 |ch is no|t includ|
|00000070| 65 64 20 61 73 20 61 20 | 73 74 61 6e 64 61 72 64 |ed as a |standard|
|00000080| 20 70 61 72 74 20 6f 66 | 20 74 68 65 20 53 49 41 | part of| the SIA|
|00000090| 4d 20 6d 61 63 72 6f 73 | 2e 0a 0a 5c 6e 65 77 74 |M macros|...\newt|
|000000a0| 68 65 6f 72 65 6d 7b 64 | 65 66 69 6e 69 74 69 6f |heorem{d|efinitio|
|000000b0| 6e 7d 5b 74 68 65 6f 72 | 65 6d 5d 7b 44 65 66 69 |n}[theor|em]{Defi|
|000000c0| 6e 69 74 69 6f 6e 7d 0a | 0a 25 20 64 65 66 69 6e |nition}.|.% defin|
|000000d0| 69 74 69 6f 6e 73 20 75 | 73 65 64 20 62 79 20 69 |itions u|sed by i|
|000000e0| 6e 63 6c 75 64 65 64 20 | 61 72 74 69 63 6c 65 73 |ncluded |articles|
|000000f0| 2c 20 72 65 70 72 6f 64 | 75 63 65 64 20 68 65 72 |, reprod|uced her|
|00000100| 65 20 66 6f 72 20 0a 25 | 20 65 64 75 63 61 74 69 |e for .%| educati|
|00000110| 6f 6e 61 6c 20 62 65 6e | 65 66 69 74 2c 20 61 6e |onal ben|efit, an|
|00000120| 64 20 74 6f 20 6d 69 6e | 69 6d 69 7a 65 20 61 6c |d to min|imize al|
|00000130| 74 65 72 61 74 69 6f 6e | 73 20 6e 65 65 64 65 64 |teration|s needed|
|00000140| 20 74 6f 20 62 65 20 6d | 61 64 65 0a 25 20 69 6e | to be m|ade.% in|
|00000150| 20 64 65 76 65 6c 6f 70 | 69 6e 67 20 74 68 69 73 | develop|ing this|
|00000160| 20 73 61 6d 70 6c 65 20 | 66 69 6c 65 2e 0a 0a 5c | sample |file...\|
|00000170| 6e 65 77 63 6f 6d 6d 61 | 6e 64 7b 5c 70 65 7d 7b |newcomma|nd{\pe}{|
|00000180| 5c 70 73 69 7d 0a 5c 64 | 65 66 5c 64 7b 5c 64 65 |\psi}.\d|ef\d{\de|
|00000190| 6c 74 61 7d 20 0a 5c 64 | 65 66 5c 64 73 7b 5c 64 |lta} .\d|ef\ds{\d|
|000001a0| 69 73 70 6c 61 79 73 74 | 79 6c 65 7d 20 0a 5c 64 |isplayst|yle} .\d|
|000001b0| 65 66 5c 65 7b 7b 5c 65 | 70 73 69 6c 6f 6e 7d 7d |ef\e{{\e|psilon}}|
|000001c0| 20 0a 5c 64 65 66 5c 65 | 62 7b 5c 62 61 72 7b 5c | .\def\e|b{\bar{\|
|000001d0| 65 74 61 7d 7d 20 20 0a | 5c 64 65 66 5c 65 6e 6f |eta}} .|\def\eno|
|000001e0| 72 6d 23 31 7b 5c 7c 23 | 31 5c 7c 5f 32 7d 20 0a |rm#1{\|#|1\|_2} .|
|000001f0| 5c 64 65 66 5c 46 70 7b | 46 5e 5c 70 72 69 6d 65 |\def\Fp{|F^\prime|
|00000200| 7d 20 20 0a 5c 64 65 66 | 5c 66 69 73 68 70 61 63 |} .\def|\fishpac|
|00000210| 6b 7b 7b 46 49 53 48 50 | 41 43 4b 7d 7d 20 0a 5c |k{{FISHP|ACK}} .\|
|00000220| 64 65 66 5c 66 6f 72 74 | 72 61 6e 7b 7b 46 4f 52 |def\fort|ran{{FOR|
|00000230| 54 52 41 4e 7d 7d 20 0a | 5c 64 65 66 5c 67 6d 72 |TRAN}} .|\def\gmr|
|00000240| 65 73 7b 7b 47 4d 52 45 | 53 7d 7d 20 0a 5c 64 65 |es{{GMRE|S}} .\de|
|00000250| 66 5c 67 6d 72 65 73 6d | 7b 7b 5c 72 6d 20 47 4d |f\gmresm|{{\rm GM|
|00000260| 52 45 53 28 24 6d 24 29 | 7d 7d 20 0a 5c 64 65 66 |RES($m$)|}} .\def|
|00000270| 5c 4b 63 7b 7b 5c 63 61 | 6c 20 4b 7d 7d 20 0a 5c |\Kc{{\ca|l K}} .\|
|00000280| 64 65 66 5c 6e 6f 72 6d | 23 31 7b 5c 7c 23 31 5c |def\norm|#1{\|#1\|
|00000290| 7c 7d 20 0a 5c 64 65 66 | 5c 77 62 7b 7b 5c 62 61 ||} .\def|\wb{{\ba|
|000002a0| 72 20 77 7d 7d 20 0a 5c | 64 65 66 5c 7a 62 7b 7b |r w}} .\|def\zb{{|
|000002b0| 5c 62 61 72 20 7a 7d 7d | 20 0a 0a 25 20 73 6f 6d |\bar z}}| ..% som|
|000002c0| 65 20 64 65 66 69 6e 69 | 74 69 6f 6e 73 20 6f 66 |e defini|tions of|
|000002d0| 20 62 6f 6c 64 20 6d 61 | 74 68 20 69 74 61 6c 69 | bold ma|th itali|
|000002e0| 63 73 20 74 6f 20 6d 61 | 6b 65 20 74 79 70 69 6e |cs to ma|ke typin|
|000002f0| 67 20 65 61 73 69 65 72 | 2e 0a 25 20 54 68 65 79 |g easier|..% They|
|00000300| 20 61 72 65 20 75 73 65 | 64 20 69 6e 20 74 68 65 | are use|d in the|
|00000310| 20 63 6f 72 6f 6c 6c 61 | 72 79 2e 0a 0a 5c 64 65 | corolla|ry...\de|
|00000320| 66 5c 62 66 45 7b 5c 6d | 62 6f 78 7b 5c 62 6f 6c |f\bfE{\m|box{\bol|
|00000330| 64 6d 61 74 68 24 45 24 | 7d 7d 0a 5c 64 65 66 5c |dmath$E$|}}.\def\|
|00000340| 62 66 47 7b 5c 6d 62 6f | 78 7b 5c 62 6f 6c 64 6d |bfG{\mbo|x{\boldm|
|00000350| 61 74 68 24 47 24 7d 7d | 0a 0a 5c 74 69 74 6c 65 |ath$G$}}|..\title|
|00000360| 7b 53 61 6d 70 6c 65 20 | 66 69 6c 65 20 66 6f 72 |{Sample |file for|
|00000370| 20 53 49 41 4d 20 5c 4c | 61 54 65 58 5c 20 6d 61 | SIAM \L|aTeX\ ma|
|00000380| 63 72 6f 20 70 61 63 6b | 61 67 65 5c 74 68 61 6e |cro pack|age\than|
|00000390| 6b 73 7b 54 68 69 73 20 | 0a 20 20 20 20 20 20 20 |ks{This |. |
|000003a0| 20 77 6f 72 6b 20 77 61 | 73 20 73 75 70 70 6f 72 | work wa|s suppor|
|000003b0| 74 65 64 20 62 79 20 74 | 68 65 20 53 6f 63 69 65 |ted by t|he Socie|
|000003c0| 74 79 20 66 6f 72 20 49 | 6e 64 75 73 74 72 69 61 |ty for I|ndustria|
|000003d0| 6c 20 61 6e 64 0a 20 20 | 20 20 20 20 20 20 41 70 |l and. | Ap|
|000003e0| 70 6c 69 65 64 20 4d 61 | 74 68 65 6d 61 74 69 63 |plied Ma|thematic|
|000003f0| 73 2c 20 50 68 69 6c 61 | 64 65 6c 70 68 69 61 2c |s, Phila|delphia,|
|00000400| 20 50 65 6e 6e 73 79 6c | 76 61 6e 69 61 2e 7d 7d | Pennsyl|vania.}}|
|00000410| 0a 0a 25 20 54 68 65 20 | 74 68 61 6e 6b 73 20 6c |..% The |thanks l|
|00000420| 69 6e 65 20 69 6e 20 74 | 68 65 20 74 69 74 6c 65 |ine in t|he title|
|00000430| 20 73 68 6f 75 6c 64 20 | 62 65 20 66 69 6c 6c 65 | should |be fille|
|00000440| 64 20 69 6e 20 69 66 20 | 74 68 65 72 65 20 69 73 |d in if |there is|
|00000450| 0a 25 20 61 6e 79 20 73 | 75 70 70 6f 72 74 20 61 |.% any s|upport a|
|00000460| 63 6b 6e 6f 77 6c 65 64 | 67 65 6d 65 6e 74 20 66 |cknowled|gement f|
|00000470| 6f 72 20 74 68 65 20 6f | 76 65 72 61 6c 6c 20 77 |or the o|verall w|
|00000480| 6f 72 6b 20 74 6f 20 62 | 65 20 69 6e 63 6c 75 64 |ork to b|e includ|
|00000490| 65 64 0a 25 20 54 68 69 | 73 20 5c 74 68 61 6e 6b |ed.% Thi|s \thank|
|000004a0| 73 20 69 73 20 61 6c 73 | 6f 20 75 73 65 64 20 66 |s is als|o used f|
|000004b0| 6f 72 20 74 68 65 20 72 | 65 63 65 69 76 65 64 20 |or the r|eceived |
|000004c0| 62 79 20 64 61 74 65 20 | 69 6e 66 6f 2c 20 62 75 |by date |info, bu|
|000004d0| 74 0a 25 20 61 75 74 68 | 6f 72 73 20 61 72 65 20 |t.% auth|ors are |
|000004e0| 6e 6f 74 20 65 78 70 65 | 63 74 65 64 20 74 6f 20 |not expe|cted to |
|000004f0| 70 72 6f 76 69 64 65 20 | 74 68 69 73 2e 0a 0a 5c |provide |this...\|
|00000500| 61 75 74 68 6f 72 7b 50 | 61 75 6c 20 44 75 67 67 |author{P|aul Dugg|
|00000510| 61 6e 5c 74 68 61 6e 6b | 73 7b 43 6f 6d 70 6f 73 |an\thank|s{Compos|
|00000520| 69 74 69 6f 6e 20 44 65 | 70 61 72 74 6d 65 6e 74 |ition De|partment|
|00000530| 2c 20 53 6f 63 69 65 74 | 79 0a 20 20 20 20 20 20 |, Societ|y. |
|00000540| 20 20 66 6f 72 20 49 6e | 64 75 73 74 72 69 61 6c | for In|dustrial|
|00000550| 20 61 6e 64 20 41 70 70 | 6c 69 65 64 20 4d 61 74 | and App|lied Mat|
|00000560| 68 65 6d 61 74 69 63 73 | 2c 20 33 36 30 30 20 55 |hematics|, 3600 U|
|00000570| 6e 69 76 65 72 69 73 74 | 79 20 0a 20 20 20 20 20 |niverist|y . |
|00000580| 20 20 20 43 69 74 79 20 | 53 63 69 65 6e 63 65 20 | City |Science |
|00000590| 43 65 6e 74 65 72 2c 20 | 50 68 69 6c 61 64 65 6c |Center, |Philadel|
|000005a0| 70 68 69 61 2c 20 50 65 | 6e 6e 73 79 6c 76 61 6e |phia, Pe|nnsylvan|
|000005b0| 69 61 2c 20 0a 20 20 20 | 20 20 20 20 20 31 39 31 |ia, . | 191|
|000005c0| 30 34 2d 32 36 38 38 20 | 28 7b 5c 74 74 20 64 75 |04-2688 |({\tt du|
|000005d0| 67 67 61 6e 40 73 69 61 | 6d 2e 6f 72 67 7d 29 2e |ggan@sia|m.org}).|
|000005e0| 7d 0a 20 20 20 20 20 20 | 20 20 5c 61 6e 64 20 56 |}. | \and V|
|000005f0| 61 72 69 6f 75 73 20 41 | 2e 7e 55 2e 20 54 68 6f |arious A|.~U. Tho|
|00000600| 72 73 5c 74 68 61 6e 6b | 73 7b 56 61 72 69 6f 75 |rs\thank|s{Variou|
|00000610| 73 20 41 66 66 69 6c 69 | 61 74 69 6f 6e 73 2c 20 |s Affili|ations, |
|00000620| 0a 20 20 20 20 20 20 20 | 20 73 75 70 70 6f 72 74 |. | support|
|00000630| 65 64 20 62 79 20 76 61 | 72 69 6f 75 73 20 66 6f |ed by va|rious fo|
|00000640| 75 6e 64 61 74 69 6f 6e | 20 67 72 61 6e 74 73 2e |undation| grants.|
|00000650| 7d 7d 0a 0a 5c 62 65 67 | 69 6e 7b 64 6f 63 75 6d |}}..\beg|in{docum|
|00000660| 65 6e 74 7d 0a 0a 5c 6d | 61 6b 65 74 69 74 6c 65 |ent}..\m|aketitle|
|00000670| 0a 0a 5c 62 65 67 69 6e | 7b 61 62 73 74 72 61 63 |..\begin|{abstrac|
|00000680| 74 7d 0a 41 6e 20 65 78 | 61 6d 70 6c 65 20 6f 66 |t}.An ex|ample of|
|00000690| 20 53 49 41 4d 20 5c 4c | 61 54 65 58 5c 20 6d 61 | SIAM \L|aTeX\ ma|
|000006a0| 63 72 6f 73 20 69 73 20 | 70 72 65 73 65 6e 74 65 |cros is |presente|
|000006b0| 64 2e 20 56 61 72 69 6f | 75 73 0a 61 73 70 65 63 |d. Vario|us.aspec|
|000006c0| 74 73 20 6f 66 20 63 6f | 6d 70 6f 73 69 6e 67 20 |ts of co|mposing |
|000006d0| 6d 61 6e 75 73 63 72 69 | 70 74 73 20 66 6f 72 20 |manuscri|pts for |
|000006e0| 53 49 41 4d 27 73 20 6a | 6f 75 72 6e 61 6c 20 73 |SIAM's j|ournal s|
|000006f0| 65 72 69 65 73 0a 61 72 | 65 20 69 6c 6c 75 73 74 |eries.ar|e illust|
|00000700| 72 61 74 65 64 20 77 69 | 74 68 20 61 63 74 75 61 |rated wi|th actua|
|00000710| 6c 20 65 78 61 6d 70 6c | 65 73 20 66 72 6f 6d 20 |l exampl|es from |
|00000720| 61 63 63 65 70 74 65 64 | 0a 6d 61 6e 75 73 63 72 |accepted|.manuscr|
|00000730| 69 70 74 73 2e 20 53 49 | 41 4d 27 73 20 73 74 79 |ipts. SI|AM's sty|
|00000740| 6c 69 73 74 69 63 20 73 | 74 61 6e 64 61 72 64 73 |listic s|tandards|
|00000750| 20 61 72 65 20 61 64 68 | 65 72 65 64 20 74 6f 0a | are adh|ered to.|
|00000760| 74 68 72 6f 75 67 68 6f | 75 74 2c 20 61 6e 64 20 |througho|ut, and |
|00000770| 69 6c 6c 75 73 74 72 61 | 74 65 64 2e 0a 5c 65 6e |illustra|ted..\en|
|00000780| 64 7b 61 62 73 74 72 61 | 63 74 7d 0a 0a 5c 62 65 |d{abstra|ct}..\be|
|00000790| 67 69 6e 7b 6b 65 79 77 | 6f 72 64 73 7d 20 0a 73 |gin{keyw|ords} .s|
|000007a0| 69 67 6e 2d 6e 6f 6e 73 | 69 6e 67 75 6c 61 72 20 |ign-nons|ingular |
|000007b0| 6d 61 74 72 69 78 2c 20 | 4c 55 2d 66 61 63 74 6f |matrix, |LU-facto|
|000007c0| 72 69 7a 61 74 69 6f 6e | 2c 20 69 6e 64 69 63 61 |rization|, indica|
|000007d0| 74 6f 72 0a 70 6f 6c 79 | 6e 6f 6d 69 61 6c 20 0a |tor.poly|nomial .|
|000007e0| 5c 65 6e 64 7b 6b 65 79 | 77 6f 72 64 73 7d 0a 0a |\end{key|words}..|
|000007f0| 5c 62 65 67 69 6e 7b 41 | 4d 53 7d 0a 31 35 41 31 |\begin{A|MS}.15A1|
|00000800| 35 2c 20 31 35 41 30 39 | 2c 20 31 35 41 32 33 0a |5, 15A09|, 15A23.|
|00000810| 5c 65 6e 64 7b 41 4d 53 | 7d 0a 0a 5c 70 61 67 65 |\end{AMS|}..\page|
|00000820| 73 74 79 6c 65 7b 6d 79 | 68 65 61 64 69 6e 67 73 |style{my|headings|
|00000830| 7d 0a 5c 74 68 69 73 70 | 61 67 65 73 74 79 6c 65 |}.\thisp|agestyle|
|00000840| 7b 70 6c 61 69 6e 7d 0a | 5c 6d 61 72 6b 62 6f 74 |{plain}.|\markbot|
|00000850| 68 7b 50 2e 20 44 55 47 | 47 41 4e 20 41 4e 44 20 |h{P. DUG|GAN AND |
|00000860| 56 2e 20 41 2e 20 55 2e | 20 54 48 4f 52 53 7d 7b |V. A. U.| THORS}{|
|00000870| 53 49 41 4d 20 4d 41 43 | 52 4f 20 45 58 41 4d 50 |SIAM MAC|RO EXAMP|
|00000880| 4c 45 53 7d 0a 0a 0a 5c | 73 65 63 74 69 6f 6e 7b |LES}...\|section{|
|00000890| 49 6e 74 72 6f 64 75 63 | 74 69 6f 6e 20 61 6e 64 |Introduc|tion and|
|000008a0| 20 65 78 61 6d 70 6c 65 | 73 7d 0a 54 68 69 73 20 | example|s}.This |
|000008b0| 70 61 70 65 72 20 70 72 | 65 73 65 6e 74 73 20 61 |paper pr|esents a|
|000008c0| 20 73 61 6d 70 6c 65 20 | 66 69 6c 65 20 66 6f 72 | sample |file for|
|000008d0| 20 74 68 65 20 75 73 65 | 20 6f 66 20 53 49 41 4d | the use| of SIAM|
|000008e0| 27 73 0a 5c 4c 61 54 65 | 58 5c 20 6d 61 63 72 6f |'s.\LaTe|X\ macro|
|000008f0| 20 70 61 63 6b 61 67 65 | 2e 20 49 74 20 69 6c 6c | package|. It ill|
|00000900| 75 73 74 72 61 74 65 73 | 20 74 68 65 20 66 65 61 |ustrates| the fea|
|00000910| 74 75 72 65 73 20 6f 66 | 20 74 68 65 0a 6d 61 63 |tures of| the.mac|
|00000920| 72 6f 20 70 61 63 6b 61 | 67 65 2c 20 75 73 69 6e |ro packa|ge, usin|
|00000930| 67 20 61 63 74 75 61 6c | 20 65 78 61 6d 70 6c 65 |g actual| example|
|00000940| 73 20 63 75 6c 6c 65 64 | 20 66 72 6f 6d 20 76 61 |s culled| from va|
|00000950| 72 69 6f 75 73 0a 70 61 | 70 65 72 73 20 70 75 62 |rious.pa|pers pub|
|00000960| 6c 69 73 68 65 64 20 69 | 6e 20 53 49 41 4d 27 73 |lished i|n SIAM's|
|00000970| 20 6a 6f 75 72 6e 61 6c | 73 2e 20 49 74 20 69 73 | journal|s. It is|
|00000980| 20 74 6f 20 62 65 20 65 | 78 70 65 63 74 65 64 0a | to be e|xpected.|
|00000990| 74 68 61 74 20 74 68 69 | 73 20 73 61 6d 70 6c 65 |that thi|s sample|
|000009a0| 20 77 69 6c 6c 20 70 72 | 6f 76 69 64 65 20 65 78 | will pr|ovide ex|
|000009b0| 61 6d 70 6c 65 73 20 6f | 66 20 68 6f 77 20 74 6f |amples o|f how to|
|000009c0| 20 75 73 65 20 74 68 65 | 0a 6d 61 63 72 6f 73 20 | use the|.macros |
|000009d0| 74 6f 20 67 65 6e 65 72 | 61 74 65 20 73 74 61 6e |to gener|ate stan|
|000009e0| 64 61 72 64 20 65 6c 65 | 6d 65 6e 74 73 20 6f 66 |dard ele|ments of|
|000009f0| 20 6a 6f 75 72 6e 61 6c | 20 70 61 70 65 72 73 2c | journal| papers,|
|00000a00| 0a 65 2e 67 2e 2c 20 74 | 68 65 6f 72 65 6d 73 2c |.e.g., t|heorems,|
|00000a10| 20 64 65 66 69 6e 69 74 | 69 6f 6e 73 2c 20 6f 72 | definit|ions, or|
|00000a20| 20 66 69 67 75 72 65 73 | 2e 20 54 68 69 73 20 70 | figures|. This p|
|00000a30| 61 70 65 72 20 61 6c 73 | 6f 0a 73 65 72 76 65 73 |aper als|o.serves|
|00000a40| 20 61 73 20 61 6e 20 65 | 78 61 6d 70 6c 65 20 6f | as an e|xample o|
|00000a50| 66 20 53 49 41 4d 27 73 | 20 73 74 79 6c 69 73 74 |f SIAM's| stylist|
|00000a60| 69 63 20 70 72 65 66 65 | 72 65 6e 63 65 73 20 66 |ic prefe|rences f|
|00000a70| 6f 72 0a 74 68 65 20 66 | 6f 72 6d 61 74 74 69 6e |or.the f|ormattin|
|00000a80| 67 20 6f 66 20 73 75 63 | 68 20 65 6c 65 6d 65 6e |g of suc|h elemen|
|00000a90| 74 73 20 61 73 20 62 69 | 62 6c 69 6f 67 72 61 70 |ts as bi|bliograp|
|00000aa0| 68 69 63 20 72 65 66 65 | 72 65 6e 63 65 73 2c 0a |hic refe|rences,.|
|00000ab0| 64 69 73 70 6c 61 79 65 | 64 20 65 71 75 61 74 69 |displaye|d equati|
|00000ac0| 6f 6e 73 2c 20 61 6e 64 | 20 65 71 75 61 74 69 6f |ons, and| equatio|
|00000ad0| 6e 20 61 72 72 61 79 73 | 2c 20 61 6d 6f 6e 67 20 |n arrays|, among |
|00000ae0| 6f 74 68 65 72 73 2e 0a | 53 6f 6d 65 20 73 70 65 |others..|Some spe|
|00000af0| 63 69 61 6c 20 63 69 72 | 63 75 6d 73 74 61 6e 63 |cial cir|cumstanc|
|00000b00| 65 73 20 61 72 65 20 6e | 6f 74 20 64 65 61 6c 74 |es are n|ot dealt|
|00000b10| 20 77 69 74 68 20 69 6e | 20 74 68 69 73 0a 73 61 | with in| this.sa|
|00000b20| 6d 70 6c 65 20 66 69 6c | 65 3b 20 66 6f 72 20 73 |mple fil|e; for s|
|00000b30| 75 63 68 20 69 6e 66 6f | 72 6d 61 74 69 6f 6e 20 |uch info|rmation |
|00000b40| 6f 6e 65 20 73 68 6f 75 | 6c 64 20 73 65 65 20 74 |one shou|ld see t|
|00000b50| 68 65 0a 69 6e 63 6c 75 | 64 65 64 20 64 6f 63 75 |he.inclu|ded docu|
|00000b60| 6d 65 6e 74 61 74 69 6f | 6e 20 66 69 6c 65 2e 0a |mentatio|n file..|
|00000b70| 0a 7b 5c 65 6d 20 4e 6f | 74 65 3a 7d 20 54 68 69 |.{\em No|te:} Thi|
|00000b80| 73 20 70 61 70 65 72 20 | 69 73 20 6e 6f 74 20 74 |s paper |is not t|
|00000b90| 6f 20 62 65 20 72 65 61 | 64 20 69 6e 20 61 6e 79 |o be rea|d in any|
|00000ba0| 20 66 6f 72 6d 20 66 6f | 72 20 63 6f 6e 74 65 6e | form fo|r conten|
|00000bb0| 74 2e 20 0a 54 68 65 20 | 63 6f 6e 67 6c 6f 6d 65 |t. .The |conglome|
|00000bc0| 72 61 74 69 6f 6e 20 6f | 66 20 65 71 75 61 74 69 |ration o|f equati|
|00000bd0| 6f 6e 73 2c 20 6c 65 6d | 6d 61 73 2c 20 61 6e 64 |ons, lem|mas, and|
|00000be0| 20 6f 74 68 65 72 20 74 | 65 78 74 20 65 6c 65 6d | other t|ext elem|
|00000bf0| 65 6e 74 73 20 77 65 72 | 65 20 0a 70 75 74 20 74 |ents wer|e .put t|
|00000c00| 6f 67 65 74 68 65 72 20 | 73 6f 6c 65 6c 79 20 66 |ogether |solely f|
|00000c10| 6f 72 20 74 79 70 6f 67 | 72 61 70 68 69 63 20 69 |or typog|raphic i|
|00000c20| 6c 6c 75 73 74 72 61 74 | 69 76 65 20 70 75 72 70 |llustrat|ive purp|
|00000c30| 6f 73 65 73 20 61 6e 64 | 20 64 6f 6e 27 74 20 0a |oses and| don't .|
|00000c40| 6d 61 6b 65 20 61 6e 79 | 20 73 65 6e 73 65 20 61 |make any| sense a|
|00000c50| 73 20 6c 65 6d 6d 61 73 | 2c 20 65 71 75 61 74 69 |s lemmas|, equati|
|00000c60| 6f 6e 73 2c 20 65 74 63 | 2e 0a 0a 5c 73 75 62 73 |ons, etc|...\subs|
|00000c70| 65 63 74 69 6f 6e 7b 53 | 61 6d 70 6c 65 20 74 65 |ection{S|ample te|
|00000c80| 78 74 7d 0a 4c 65 74 20 | 24 53 3d 5b 73 5f 7b 69 |xt}.Let |$S=[s_{i|
|00000c90| 6a 7d 5d 24 20 28 24 31 | 5c 6c 65 71 20 69 2c 6a |j}]$ ($1|\leq i,j|
|00000ca0| 5c 6c 65 71 20 6e 24 29 | 20 62 65 20 61 20 24 28 |\leq n$)| be a $(|
|00000cb0| 30 2c 31 2c 2d 31 29 24 | 2d 6d 61 74 72 69 78 0a |0,1,-1)$|-matrix.|
|00000cc0| 6f 66 20 6f 72 64 65 72 | 20 24 6e 24 2e 20 54 68 |of order| $n$. Th|
|00000cd0| 65 6e 20 24 53 24 20 69 | 73 20 61 20 7b 5c 65 6d |en $S$ i|s a {\em|
|00000ce0| 20 73 69 67 6e 2d 6e 6f | 6e 73 69 6e 67 75 6c 61 | sign-no|nsingula|
|00000cf0| 72 20 6d 61 74 72 69 78 | 7d 0a 28 53 4e 53 2d 6d |r matrix|}.(SNS-m|
|00000d00| 61 74 72 69 78 29 20 70 | 72 6f 76 69 64 65 64 20 |atrix) p|rovided |
|00000d10| 74 68 61 74 20 65 61 63 | 68 20 72 65 61 6c 20 6d |that eac|h real m|
|00000d20| 61 74 72 69 78 20 77 69 | 74 68 20 74 68 65 20 73 |atrix wi|th the s|
|00000d30| 61 6d 65 0a 73 69 67 6e | 20 70 61 74 74 65 72 6e |ame.sign| pattern|
|00000d40| 20 61 73 20 24 53 24 20 | 69 73 20 6e 6f 6e 73 69 | as $S$ |is nonsi|
|00000d50| 6e 67 75 6c 61 72 2e 20 | 54 68 65 72 65 20 68 61 |ngular. |There ha|
|00000d60| 73 20 62 65 65 6e 0a 63 | 6f 6e 73 69 64 65 72 61 |s been.c|onsidera|
|00000d70| 62 6c 65 20 72 65 63 65 | 6e 74 20 69 6e 74 65 72 |ble rece|nt inter|
|00000d80| 65 73 74 20 69 6e 20 63 | 6f 6e 73 74 72 75 63 74 |est in c|onstruct|
|00000d90| 69 6e 67 20 61 6e 64 0a | 63 68 61 72 61 63 74 65 |ing and.|characte|
|00000da0| 72 69 7a 69 6e 67 20 53 | 4e 53 2d 6d 61 74 72 69 |rizing S|NS-matri|
|00000db0| 63 65 73 20 5c 63 69 74 | 65 7b 62 73 7d 2c 20 5c |ces \cit|e{bs}, \|
|00000dc0| 63 69 74 65 7b 6b 6c 6d | 7d 2e 20 54 68 65 72 65 |cite{klm|}. There|
|00000dd0| 0a 68 61 73 20 61 6c 73 | 6f 20 62 65 65 6e 20 69 |.has als|o been i|
|00000de0| 6e 74 65 72 65 73 74 20 | 69 6e 20 73 74 72 6f 6e |nterest |in stron|
|00000df0| 67 20 66 6f 72 6d 73 20 | 6f 66 0a 73 69 67 6e 2d |g forms |of.sign-|
|00000e00| 6e 6f 6e 73 69 6e 67 75 | 6c 61 72 69 74 79 20 5c |nonsingu|larity \|
|00000e10| 63 69 74 65 7b 64 6a 64 | 7d 2e 20 49 6e 20 74 68 |cite{djd|}. In th|
|00000e20| 69 73 20 70 61 70 65 72 | 20 77 65 20 67 69 76 65 |is paper| we give|
|00000e30| 20 61 20 6e 65 77 0a 67 | 65 6e 65 72 61 6c 69 7a | a new.g|eneraliz|
|00000e40| 61 74 69 6f 6e 20 6f 66 | 20 53 4e 53 2d 6d 61 74 |ation of| SNS-mat|
|00000e50| 72 69 63 65 73 20 61 6e | 64 20 69 6e 76 65 73 74 |rices an|d invest|
|00000e60| 69 67 61 74 65 20 73 6f | 6d 65 20 6f 66 0a 74 68 |igate so|me of.th|
|00000e70| 65 69 72 20 62 61 73 69 | 63 20 70 72 6f 70 65 72 |eir basi|c proper|
|00000e80| 74 69 65 73 2e 0a 20 0a | 4c 65 74 20 24 53 3d 5b |ties.. .|Let $S=[|
|00000e90| 73 5f 7b 69 6a 7d 5d 24 | 20 62 65 20 61 20 24 28 |s_{ij}]$| be a $(|
|00000ea0| 30 2c 31 2c 2d 31 29 24 | 2d 6d 61 74 72 69 78 20 |0,1,-1)$|-matrix |
|00000eb0| 6f 66 20 6f 72 64 65 72 | 20 24 6e 24 20 61 6e 64 |of order| $n$ and|
|00000ec0| 0a 6c 65 74 20 24 43 3d | 5b 63 5f 7b 69 6a 7d 5d |.let $C=|[c_{ij}]|
|00000ed0| 24 20 62 65 20 61 20 72 | 65 61 6c 20 6d 61 74 72 |$ be a r|eal matr|
|00000ee0| 69 78 20 6f 66 20 6f 72 | 64 65 72 20 24 6e 24 2e |ix of or|der $n$.|
|00000ef0| 20 54 68 65 20 70 61 69 | 72 0a 24 28 53 2c 43 29 | The pai|r.$(S,C)|
|00000f00| 24 20 69 73 20 63 61 6c | 6c 65 64 20 61 20 7b 5c |$ is cal|led a {\|
|00000f10| 65 6d 20 6d 61 74 72 69 | 78 20 70 61 69 72 20 6f |em matri|x pair o|
|00000f20| 66 20 6f 72 64 65 72 7d | 20 24 6e 24 2e 0a 54 68 |f order}| $n$..Th|
|00000f30| 72 6f 75 67 68 6f 75 74 | 2c 20 24 58 3d 5b 78 5f |roughout|, $X=[x_|
|00000f40| 7b 69 6a 7d 5d 24 20 64 | 65 6e 6f 74 65 73 20 61 |{ij}]$ d|enotes a|
|00000f50| 20 6d 61 74 72 69 78 20 | 6f 66 20 6f 72 64 65 72 | matrix |of order|
|00000f60| 20 24 6e 24 0a 77 68 6f | 73 65 20 65 6e 74 72 69 | $n$.who|se entri|
|00000f70| 65 73 20 61 72 65 20 61 | 6c 67 65 62 72 61 69 63 |es are a|lgebraic|
|00000f80| 61 6c 6c 79 20 69 6e 64 | 65 70 65 6e 64 65 6e 74 |ally ind|ependent|
|00000f90| 20 69 6e 64 65 74 65 72 | 6d 69 6e 61 74 65 73 0a | indeter|minates.|
|00000fa0| 6f 76 65 72 20 74 68 65 | 20 72 65 61 6c 20 66 69 |over the| real fi|
|00000fb0| 65 6c 64 2e 20 4c 65 74 | 20 24 53 5c 63 69 72 63 |eld. Let| $S\circ|
|00000fc0| 20 58 24 20 64 65 6e 6f | 74 65 20 74 68 65 20 48 | X$ deno|te the H|
|00000fd0| 61 64 61 6d 61 72 64 0a | 70 72 6f 64 75 63 74 20 |adamard.|product |
|00000fe0| 28 65 6e 74 72 79 77 69 | 73 65 20 70 72 6f 64 75 |(entrywi|se produ|
|00000ff0| 63 74 29 20 6f 66 20 24 | 53 24 20 61 6e 64 20 24 |ct) of $|S$ and $|
|00001000| 58 24 2e 20 57 65 20 73 | 61 79 20 74 68 61 74 20 |X$. We s|ay that |
|00001010| 74 68 65 0a 70 61 69 72 | 20 24 28 53 2c 43 29 24 |the.pair| $(S,C)$|
|00001020| 20 69 73 20 61 20 7b 5c | 65 6d 20 73 69 67 6e 2d | is a {\|em sign-|
|00001030| 6e 6f 6e 73 69 6e 67 75 | 6c 61 72 20 6d 61 74 72 |nonsingu|lar matr|
|00001040| 69 78 20 70 61 69 72 20 | 6f 66 0a 6f 72 64 65 72 |ix pair |of.order|
|00001050| 7d 20 24 6e 24 2c 20 61 | 62 62 72 65 76 69 61 74 |} $n$, a|bbreviat|
|00001060| 65 64 20 53 4e 53 2d 7b | 5c 65 6d 20 6d 61 74 72 |ed SNS-{|\em matr|
|00001070| 69 78 20 70 61 69 72 20 | 6f 66 20 6f 72 64 65 72 |ix pair |of order|
|00001080| 7d 20 24 6e 24 2c 0a 70 | 72 6f 76 69 64 65 64 20 |} $n$,.p|rovided |
|00001090| 74 68 61 74 20 74 68 65 | 20 6d 61 74 72 69 78 20 |that the| matrix |
|000010a0| 5c 5b 41 3d 53 5c 63 69 | 72 63 20 58 2b 43 5c 5d |\[A=S\ci|rc X+C\]|
|000010b0| 20 69 73 20 6e 6f 6e 73 | 69 6e 67 75 6c 61 72 0a | is nons|ingular.|
|000010c0| 66 6f 72 20 61 6c 6c 20 | 70 6f 73 69 74 69 76 65 |for all |positive|
|000010d0| 20 72 65 61 6c 20 76 61 | 6c 75 65 73 20 6f 66 20 | real va|lues of |
|000010e0| 74 68 65 20 24 78 5f 7b | 69 6a 7d 24 2e 20 20 49 |the $x_{|ij}$. I|
|000010f0| 66 20 24 43 3d 4f 24 0a | 74 68 65 6e 20 74 68 65 |f $C=O$.|then the|
|00001100| 20 70 61 69 72 20 24 28 | 53 2c 4f 29 24 20 69 73 | pair $(|S,O)$ is|
|00001110| 20 61 20 53 4e 53 2d 6d | 61 74 72 69 78 20 70 61 | a SNS-m|atrix pa|
|00001120| 69 72 20 69 66 20 61 6e | 64 20 6f 6e 6c 79 20 69 |ir if an|d only i|
|00001130| 66 0a 24 53 24 20 69 73 | 20 61 20 53 4e 53 2d 6d |f.$S$ is| a SNS-m|
|00001140| 61 74 72 69 78 2e 20 20 | 49 66 20 24 53 3d 4f 24 |atrix. |If $S=O$|
|00001150| 20 74 68 65 6e 20 74 68 | 65 20 70 61 69 72 20 24 | then th|e pair $|
|00001160| 28 4f 2c 43 29 24 20 69 | 73 20 61 0a 53 4e 53 2d |(O,C)$ i|s a.SNS-|
|00001170| 6d 61 74 72 69 78 20 70 | 61 69 72 20 69 66 20 61 |matrix p|air if a|
|00001180| 6e 64 20 6f 6e 6c 79 20 | 69 66 20 24 43 24 20 69 |nd only |if $C$ i|
|00001190| 73 20 6e 6f 6e 73 69 6e | 67 75 6c 61 72 2e 20 54 |s nonsin|gular. T|
|000011a0| 68 75 73 0a 53 4e 53 2d | 6d 61 74 72 69 78 20 70 |hus.SNS-|matrix p|
|000011b0| 61 69 72 73 20 69 6e 63 | 6c 75 64 65 20 62 6f 74 |airs inc|lude bot|
|000011c0| 68 20 6e 6f 6e 73 69 6e | 67 75 6c 61 72 20 6d 61 |h nonsin|gular ma|
|000011d0| 74 72 69 63 65 73 20 61 | 6e 64 0a 73 69 67 6e 2d |trices a|nd.sign-|
|000011e0| 6e 6f 6e 73 69 6e 67 75 | 6c 61 72 20 6d 61 74 72 |nonsingu|lar matr|
|000011f0| 69 63 65 73 20 61 73 20 | 73 70 65 63 69 61 6c 20 |ices as |special |
|00001200| 63 61 73 65 73 2e 0a 20 | 0a 54 68 65 20 70 61 69 |cases.. |.The pai|
|00001210| 72 73 20 24 28 53 2c 43 | 29 24 20 77 69 74 68 0a |rs $(S,C|)$ with.|
|00001220| 5c 5b 53 3d 5c 6c 65 66 | 74 5b 5c 62 65 67 69 6e |\[S=\lef|t[\begin|
|00001230| 7b 61 72 72 61 79 7d 7b | 63 63 7d 31 26 30 5c 5c |{array}{|cc}1&0\\|
|00001240| 30 26 30 5c 65 6e 64 7b | 61 72 72 61 79 7d 5c 72 |0&0\end{|array}\r|
|00001250| 69 67 68 74 5d 2c 5c 71 | 71 75 61 64 20 0a 43 3d |ight],\q|quad .C=|
|00001260| 5c 6c 65 66 74 5b 5c 62 | 65 67 69 6e 7b 61 72 72 |\left[\b|egin{arr|
|00001270| 61 79 7d 7b 63 63 7d 31 | 26 31 5c 5c 31 26 31 5c |ay}{cc}1|&1\\1&1\|
|00001280| 65 6e 64 7b 61 72 72 61 | 79 7d 5c 72 69 67 68 74 |end{arra|y}\right|
|00001290| 5d 5c 5d 20 61 6e 64 20 | 0a 5c 5b 53 3d 5c 6c 65 |]\] and |.\[S=\le|
|000012a0| 66 74 5b 5c 62 65 67 69 | 6e 7b 61 72 72 61 79 7d |ft[\begi|n{array}|
|000012b0| 7b 63 63 63 7d 31 26 31 | 26 30 5c 5c 31 26 31 26 |{ccc}1&1|&0\\1&1&|
|000012c0| 30 5c 5c 30 26 30 26 30 | 5c 65 6e 64 7b 61 72 72 |0\\0&0&0|\end{arr|
|000012d0| 61 79 7d 5c 72 69 67 68 | 74 5d 2c 5c 71 71 75 61 |ay}\righ|t],\qqua|
|000012e0| 64 20 0a 43 3d 5c 6c 65 | 66 74 5b 5c 62 65 67 69 |d .C=\le|ft[\begi|
|000012f0| 6e 7b 61 72 72 61 79 7d | 7b 63 63 63 7d 30 26 30 |n{array}|{ccc}0&0|
|00001300| 26 31 5c 5c 30 26 32 26 | 30 5c 5c 0a 33 26 30 26 |&1\\0&2&|0\\.3&0&|
|00001310| 30 5c 65 6e 64 7b 61 72 | 72 61 79 7d 5c 72 69 67 |0\end{ar|ray}\rig|
|00001320| 68 74 5d 5c 5d 20 61 72 | 65 20 65 78 61 6d 70 6c |ht]\] ar|e exampl|
|00001330| 65 73 20 6f 66 20 53 4e | 53 2d 6d 61 74 72 69 78 |es of SN|S-matrix|
|00001340| 20 70 61 69 72 73 2e 0a | 0a 5c 73 75 62 73 65 63 | pairs..|.\subsec|
|00001350| 74 69 6f 6e 7b 41 20 72 | 65 6d 75 6e 65 72 61 74 |tion{A r|emunerat|
|00001360| 69 6f 6e 20 6c 69 73 74 | 7d 0a 49 6e 20 74 68 69 |ion list|}.In thi|
|00001370| 73 20 70 61 70 65 72 20 | 77 65 20 63 6f 6e 73 69 |s paper |we consi|
|00001380| 64 65 72 20 74 68 65 20 | 65 76 61 6c 75 61 74 69 |der the |evaluati|
|00001390| 6f 6e 20 6f 66 20 69 6e | 74 65 67 72 61 6c 73 20 |on of in|tegrals |
|000013a0| 6f 66 20 74 68 65 20 0a | 66 6f 6c 6c 6f 77 69 6e |of the .|followin|
|000013b0| 67 20 66 6f 72 6d 73 3a | 0a 5c 62 65 67 69 6e 7b |g forms:|.\begin{|
|000013c0| 65 71 75 61 74 69 6f 6e | 7d 0a 5c 69 6e 74 5f 61 |equation|}.\int_a|
|000013d0| 5e 62 20 5c 6c 65 66 74 | 28 20 5c 73 75 6d 5f 69 |^b \left|( \sum_i|
|000013e0| 20 45 5f 69 20 42 5f 7b | 69 2c 6b 2c 78 7d 28 74 | E_i B_{|i,k,x}(t|
|000013f0| 29 20 5c 72 69 67 68 74 | 29 0a 20 20 20 20 20 20 |) \right|). |
|00001400| 20 20 20 5c 6c 65 66 74 | 28 20 5c 73 75 6d 5f 6a | \left|( \sum_j|
|00001410| 20 46 5f 6a 20 42 5f 7b | 6a 2c 6c 2c 79 7d 28 74 | F_j B_{|j,l,y}(t|
|00001420| 29 20 5c 72 69 67 68 74 | 29 20 64 74 2c 5c 6c 61 |) \right|) dt,\la|
|00001430| 62 65 6c 7b 70 72 6f 62 | 6c 65 6d 7d 0a 5c 65 6e |bel{prob|lem}.\en|
|00001440| 64 7b 65 71 75 61 74 69 | 6f 6e 7d 0a 5c 62 65 67 |d{equati|on}.\beg|
|00001450| 69 6e 7b 65 71 75 61 74 | 69 6f 6e 7d 0a 5c 69 6e |in{equat|ion}.\in|
|00001460| 74 5f 61 5e 62 20 66 28 | 74 29 20 5c 6c 65 66 74 |t_a^b f(|t) \left|
|00001470| 28 20 5c 73 75 6d 5f 69 | 20 45 5f 69 20 42 5f 7b |( \sum_i| E_i B_{|
|00001480| 69 2c 6b 2c 78 7d 28 74 | 29 20 5c 72 69 67 68 74 |i,k,x}(t|) \right|
|00001490| 29 20 64 74 2c 5c 6c 61 | 62 65 6c 7b 70 72 6f 62 |) dt,\la|bel{prob|
|000014a0| 6c 65 6d 32 7d 0a 5c 65 | 6e 64 7b 65 71 75 61 74 |lem2}.\e|nd{equat|
|000014b0| 69 6f 6e 7d 0a 77 68 65 | 72 65 20 24 42 5f 7b 69 |ion}.whe|re $B_{i|
|000014c0| 2c 6b 2c 78 7d 24 20 69 | 73 20 74 68 65 20 24 69 |,k,x}$ i|s the $i|
|000014d0| 24 74 68 20 42 2d 73 70 | 6c 69 6e 65 20 6f 66 20 |$th B-sp|line of |
|000014e0| 6f 72 64 65 72 20 24 6b | 24 20 64 65 66 69 6e 65 |order $k|$ define|
|000014f0| 64 20 6f 76 65 72 20 74 | 68 65 0a 6b 6e 6f 74 73 |d over t|he.knots|
|00001500| 20 24 78 5f 69 2c 20 78 | 5f 7b 69 2b 31 7d 2c 20 | $x_i, x|_{i+1}, |
|00001510| 5c 6c 64 6f 74 73 2c 20 | 78 5f 7b 69 2b 6b 7d 24 |\ldots, |x_{i+k}$|
|00001520| 2e 0a 57 65 20 77 69 6c | 6c 20 63 6f 6e 73 69 64 |..We wil|l consid|
|00001530| 65 72 20 42 2d 73 70 6c | 69 6e 65 73 20 6e 6f 72 |er B-spl|ines nor|
|00001540| 6d 61 6c 69 7a 65 64 20 | 73 6f 20 74 68 61 74 20 |malized |so that |
|00001550| 74 68 65 69 72 20 69 6e | 74 65 67 72 61 6c 20 69 |their in|tegral i|
|00001560| 73 20 6f 6e 65 2e 0a 54 | 68 65 20 73 70 6c 69 6e |s one..T|he splin|
|00001570| 65 73 20 6d 61 79 20 62 | 65 20 6f 66 20 64 69 66 |es may b|e of dif|
|00001580| 66 65 72 65 6e 74 20 6f | 72 64 65 72 73 20 61 6e |ferent o|rders an|
|00001590| 64 0a 64 65 66 69 6e 65 | 64 20 6f 6e 20 64 69 66 |d.define|d on dif|
|000015a0| 66 65 72 65 6e 74 20 6b | 6e 6f 74 20 73 65 71 75 |ferent k|not sequ|
|000015b0| 65 6e 63 65 73 20 24 78 | 24 20 61 6e 64 20 24 79 |ences $x|$ and $y|
|000015c0| 24 2e 0a 4f 66 74 65 6e | 20 74 68 65 20 6c 69 6d |$..Often| the lim|
|000015d0| 69 74 73 20 6f 66 20 69 | 6e 74 65 67 72 61 74 69 |its of i|ntegrati|
|000015e0| 6f 6e 20 77 69 6c 6c 20 | 62 65 20 74 68 65 20 65 |on will |be the e|
|000015f0| 6e 74 69 72 65 20 72 65 | 61 6c 20 6c 69 6e 65 2c |ntire re|al line,|
|00001600| 20 24 2d 5c 69 6e 66 74 | 79 24 0a 74 6f 20 24 2b | $-\inft|y$.to $+|
|00001610| 5c 69 6e 66 74 79 24 2e | 20 4e 6f 74 65 20 74 68 |\infty$.| Note th|
|00001620| 61 74 20 28 5c 72 65 66 | 7b 70 72 6f 62 6c 65 6d |at (\ref|{problem|
|00001630| 7d 29 20 69 73 20 61 20 | 73 70 65 63 69 61 6c 20 |}) is a |special |
|00001640| 63 61 73 65 20 6f 66 20 | 28 5c 72 65 66 7b 70 72 |case of |(\ref{pr|
|00001650| 6f 62 6c 65 6d 32 7d 29 | 0a 77 68 65 72 65 20 24 |oblem2})|.where $|
|00001660| 66 28 74 29 24 20 69 73 | 20 61 20 73 70 6c 69 6e |f(t)$ is| a splin|
|00001670| 65 2e 0a 0a 0a 54 68 65 | 72 65 20 61 72 65 20 66 |e....The|re are f|
|00001680| 69 76 65 20 64 69 66 66 | 65 72 65 6e 74 20 6d 65 |ive diff|erent me|
|00001690| 74 68 6f 64 73 20 66 6f | 72 20 63 61 6c 63 75 6c |thods fo|r calcul|
|000016a0| 61 74 69 6e 67 20 28 5c | 72 65 66 7b 70 72 6f 62 |ating (\|ref{prob|
|000016b0| 6c 65 6d 7d 29 0a 74 68 | 61 74 20 77 69 6c 6c 20 |lem}).th|at will |
|000016c0| 62 65 20 63 6f 6e 73 69 | 64 65 72 65 64 3a 0a 5c |be consi|dered:.\|
|000016d0| 62 65 67 69 6e 7b 72 65 | 6d 75 6e 65 72 61 74 65 |begin{re|munerate|
|000016e0| 7d 0a 5c 69 74 65 6d 20 | 55 73 65 20 47 61 75 73 |}.\item |Use Gaus|
|000016f0| 73 20 71 75 61 64 72 61 | 74 75 72 65 20 6f 6e 20 |s quadra|ture on |
|00001700| 65 61 63 68 20 69 6e 74 | 65 72 76 61 6c 2e 0a 5c |each int|erval..\|
|00001710| 69 74 65 6d 20 43 6f 6e | 76 65 72 74 20 74 68 65 |item Con|vert the|
|00001720| 20 69 6e 74 65 67 72 61 | 6c 20 74 6f 20 61 20 6c | integra|l to a l|
|00001730| 69 6e 65 61 72 20 63 6f | 6d 62 69 6e 61 74 69 6f |inear co|mbinatio|
|00001740| 6e 20 6f 66 0a 20 20 20 | 20 20 20 69 6e 74 65 67 |n of. | integ|
|00001750| 72 61 6c 73 20 6f 66 20 | 70 72 6f 64 75 63 74 73 |rals of |products|
|00001760| 20 6f 66 20 42 2d 73 70 | 6c 69 6e 65 73 20 61 6e | of B-sp|lines an|
|00001770| 64 20 70 72 6f 76 69 64 | 65 20 61 20 72 65 63 75 |d provid|e a recu|
|00001780| 72 72 65 6e 63 65 20 66 | 6f 72 0a 20 20 20 20 20 |rrence f|or. |
|00001790| 20 69 6e 74 65 67 72 61 | 74 69 6e 67 20 74 68 65 | integra|ting the|
|000017a0| 20 70 72 6f 64 75 63 74 | 20 6f 66 20 61 20 70 61 | product| of a pa|
|000017b0| 69 72 20 6f 66 20 42 2d | 73 70 6c 69 6e 65 73 2e |ir of B-|splines.|
|000017c0| 0a 5c 69 74 65 6d 20 43 | 6f 6e 76 65 72 74 20 74 |.\item C|onvert t|
|000017d0| 68 65 20 73 75 6d 73 20 | 6f 66 20 42 2d 73 70 6c |he sums |of B-spl|
|000017e0| 69 6e 65 73 20 74 6f 20 | 70 69 65 63 65 77 69 73 |ines to |piecewis|
|000017f0| 65 0a 20 20 20 20 20 20 | 42 5c 27 7b 65 7d 7a 69 |e. |B\'{e}zi|
|00001800| 65 72 20 66 6f 72 6d 61 | 74 20 61 6e 64 20 69 6e |er forma|t and in|
|00001810| 74 65 67 72 61 74 65 20 | 73 65 67 6d 65 6e 74 0a |tegrate |segment.|
|00001820| 20 20 20 20 20 20 62 79 | 20 73 65 67 6d 65 6e 74 | by| segment|
|00001830| 20 75 73 69 6e 67 20 74 | 68 65 20 70 72 6f 70 65 | using t|he prope|
|00001840| 72 74 69 65 73 20 6f 66 | 20 74 68 65 20 42 65 72 |rties of| the Ber|
|00001850| 6e 73 74 65 69 6e 20 70 | 6f 6c 79 6e 6f 6d 69 61 |nstein p|olynomia|
|00001860| 6c 73 2e 0a 5c 69 74 65 | 6d 20 45 78 70 72 65 73 |ls..\ite|m Expres|
|00001870| 73 20 74 68 65 20 70 72 | 6f 64 75 63 74 20 6f 66 |s the pr|oduct of|
|00001880| 20 61 20 70 61 69 72 20 | 6f 66 20 42 2d 73 70 6c | a pair |of B-spl|
|00001890| 69 6e 65 73 20 61 73 20 | 61 20 6c 69 6e 65 61 72 |ines as |a linear|
|000018a0| 20 63 6f 6d 62 69 6e 61 | 74 69 6f 6e 0a 20 20 20 | combina|tion. |
|000018b0| 20 20 20 6f 66 20 42 2d | 73 70 6c 69 6e 65 73 2e | of B-|splines.|
|000018c0| 0a 20 20 20 20 20 20 55 | 73 65 20 74 68 69 73 20 |. U|se this |
|000018d0| 74 6f 20 72 65 66 6f 72 | 6d 75 6c 61 74 65 20 74 |to refor|mulate t|
|000018e0| 68 65 20 69 6e 74 65 67 | 72 61 6e 64 20 61 73 20 |he integ|rand as |
|000018f0| 61 20 6c 69 6e 65 61 72 | 20 63 6f 6d 62 69 6e 61 |a linear| combina|
|00001900| 74 69 6f 6e 0a 20 20 20 | 20 20 20 6f 66 20 42 2d |tion. | of B-|
|00001910| 73 70 6c 69 6e 65 73 2c | 20 61 6e 64 20 69 6e 74 |splines,| and int|
|00001920| 65 67 72 61 74 65 20 74 | 65 72 6d 20 62 79 20 74 |egrate t|erm by t|
|00001930| 65 72 6d 2e 0a 5c 69 74 | 65 6d 20 49 6e 74 65 67 |erm..\it|em Integ|
|00001940| 72 61 74 65 20 62 79 20 | 70 61 72 74 73 2e 0a 5c |rate by |parts..\|
|00001950| 65 6e 64 7b 72 65 6d 75 | 6e 65 72 61 74 65 7d 0a |end{remu|nerate}.|
|00001960| 4f 66 20 74 68 65 73 65 | 20 66 69 76 65 2c 20 6f |Of these| five, o|
|00001970| 6e 6c 79 20 6d 65 74 68 | 6f 64 73 20 31 20 61 6e |nly meth|ods 1 an|
|00001980| 64 20 35 20 61 72 65 20 | 73 75 69 74 61 62 6c 65 |d 5 are |suitable|
|00001990| 20 66 6f 72 20 63 61 6c | 63 75 6c 61 74 69 6e 67 | for cal|culating|
|000019a0| 20 0a 28 5c 72 65 66 7b | 70 72 6f 62 6c 65 6d 32 | .(\ref{|problem2|
|000019b0| 7d 29 2e 20 54 68 65 20 | 66 69 72 73 74 20 66 6f |}). The |first fo|
|000019c0| 75 72 20 6d 65 74 68 6f | 64 73 20 77 69 6c 6c 20 |ur metho|ds will |
|000019d0| 62 65 20 74 6f 75 63 68 | 65 64 20 6f 6e 20 61 6e |be touch|ed on an|
|000019e0| 64 20 74 68 65 20 0a 6c | 61 73 74 20 77 69 6c 6c |d the .l|ast will|
|000019f0| 20 62 65 20 64 69 73 63 | 75 73 73 65 64 20 61 74 | be disc|ussed at|
|00001a00| 20 6c 65 6e 67 74 68 2e | 0a 0a 0a 5c 73 75 62 73 | length.|...\subs|
|00001a10| 65 63 74 69 6f 6e 7b 53 | 6f 6d 65 20 64 69 73 70 |ection{S|ome disp|
|00001a20| 6c 61 79 65 64 20 65 71 | 75 61 74 69 6f 6e 73 20 |layed eq|uations |
|00001a30| 61 6e 64 20 7b 5c 74 74 | 20 5c 7b 65 71 6e 61 72 |and {\tt| \{eqnar|
|00001a40| 72 61 79 5c 7d 7d 73 7d | 0a 20 20 20 20 20 42 79 |ray\}}s}|. By|
|00001a50| 20 69 6e 74 72 6f 64 75 | 63 69 6e 67 20 74 68 65 | introdu|cing the|
|00001a60| 20 70 72 6f 64 75 63 74 | 20 74 6f 70 6f 6c 6f 67 | product| topolog|
|00001a70| 79 20 6f 6e 20 20 24 52 | 5e 7b 6d 20 5c 74 69 6d |y on $R|^{m \tim|
|00001a80| 65 73 20 6d 7d 20 5c 74 | 69 6d 65 73 0a 52 5e 7b |es m} \t|imes.R^{|
|00001a90| 6e 20 5c 74 69 6d 65 73 | 20 6e 7d 24 20 20 77 69 |n \times| n}$ wi|
|00001aa0| 74 68 20 74 68 65 20 69 | 6e 64 75 63 65 64 20 69 |th the i|nduced i|
|00001ab0| 6e 6e 65 72 20 70 72 6f | 64 75 63 74 0a 5c 62 65 |nner pro|duct.\be|
|00001ac0| 67 69 6e 7b 65 71 75 61 | 74 69 6f 6e 7d 0a 5c 6c |gin{equa|tion}.\l|
|00001ad0| 61 6e 67 6c 65 20 28 41 | 5f 7b 31 7d 2c 42 5f 7b |angle (A|_{1},B_{|
|00001ae0| 31 7d 29 2c 20 28 41 5f | 7b 32 7d 2c 42 5f 7b 32 |1}), (A_|{2},B_{2|
|00001af0| 7d 29 5c 72 61 6e 67 6c | 65 20 3a 3d 20 5c 6c 61 |})\rangl|e := \la|
|00001b00| 6e 67 6c 65 20 41 5f 7b | 31 7d 2c 41 5f 7b 32 7d |ngle A_{|1},A_{2}|
|00001b10| 5c 72 61 6e 67 6c 65 20 | 0a 2b 20 5c 6c 61 6e 67 |\rangle |.+ \lang|
|00001b20| 6c 65 20 42 5f 7b 31 7d | 2c 42 5f 7b 32 7d 5c 72 |le B_{1}|,B_{2}\r|
|00001b30| 61 6e 67 6c 65 2c 5c 6c | 61 62 65 6c 7b 65 71 32 |angle,\l|abel{eq2|
|00001b40| 2e 31 30 7d 0a 5c 65 6e | 64 7b 65 71 75 61 74 69 |.10}.\en|d{equati|
|00001b50| 6f 6e 7d 0a 77 65 20 63 | 61 6c 63 75 6c 61 74 65 |on}.we c|alculate|
|00001b60| 20 74 68 65 20 46 72 5c | 27 7b 65 7d 63 68 65 74 | the Fr\|'{e}chet|
|00001b70| 20 64 65 72 69 76 61 74 | 69 76 65 20 6f 66 20 20 | derivat|ive of |
|00001b80| 24 46 24 20 20 61 73 20 | 66 6f 6c 6c 6f 77 73 3a |$F$ as |follows:|
|00001b90| 0a 5c 62 65 67 69 6e 7b | 65 71 6e 61 72 72 61 79 |.\begin{|eqnarray|
|00001ba0| 7d 0a 20 46 27 28 55 2c | 56 29 28 48 2c 4b 29 20 |}. F'(U,|V)(H,K) |
|00001bb0| 26 3d 26 20 5c 6c 61 6e | 67 6c 65 20 52 28 55 2c |&=& \lan|gle R(U,|
|00001bc0| 56 29 2c 48 5c 53 69 67 | 6d 61 20 56 5e 7b 54 7d |V),H\Sig|ma V^{T}|
|00001bd0| 20 2b 20 55 5c 53 69 67 | 6d 61 20 4b 5e 7b 54 7d | + U\Sig|ma K^{T}|
|00001be0| 20 2d 0a 50 28 48 5c 53 | 69 67 6d 61 20 56 5e 7b | -.P(H\S|igma V^{|
|00001bf0| 54 7d 20 2b 20 55 5c 53 | 69 67 6d 61 20 4b 5e 7b |T} + U\S|igma K^{|
|00001c00| 54 7d 29 5c 72 61 6e 67 | 6c 65 20 5c 6e 6f 6e 75 |T})\rang|le \nonu|
|00001c10| 6d 62 65 72 20 5c 5c 0a | 20 20 20 20 20 20 20 20 |mber \\.| |
|00001c20| 20 26 3d 26 20 5c 6c 61 | 6e 67 6c 65 20 52 28 55 | &=& \la|ngle R(U|
|00001c30| 2c 56 29 2c 48 5c 53 69 | 67 6d 61 20 56 5e 7b 54 |,V),H\Si|gma V^{T|
|00001c40| 7d 20 2b 20 55 5c 53 69 | 67 6d 61 20 4b 5e 7b 54 |} + U\Si|gma K^{T|
|00001c50| 7d 5c 72 61 6e 67 6c 65 | 20 5c 6c 61 62 65 6c 7b |}\rangle| \label{|
|00001c60| 65 71 32 2e 31 31 7d 20 | 5c 5c 0a 26 3d 26 20 5c |eq2.11} |\\.&=& \|
|00001c70| 6c 61 6e 67 6c 65 20 52 | 28 55 2c 56 29 56 5c 53 |langle R|(U,V)V\S|
|00001c80| 69 67 6d 61 5e 7b 54 7d | 2c 48 5c 72 61 6e 67 6c |igma^{T}|,H\rangl|
|00001c90| 65 20 2b 20 5c 6c 61 6e | 67 6c 65 20 5c 53 69 67 |e + \lan|gle \Sig|
|00001ca0| 6d 61 5e 7b 54 7d 55 5e | 7b 54 7d 52 28 55 2c 56 |ma^{T}U^|{T}R(U,V|
|00001cb0| 29 2c 4b 5e 7b 54 7d 5c | 72 61 6e 67 6c 65 2e 20 |),K^{T}\|rangle. |
|00001cc0| 20 20 20 20 5c 6e 6f 6e | 75 6d 62 65 72 0a 5c 65 | \non|umber.\e|
|00001cd0| 6e 64 7b 65 71 6e 61 72 | 72 61 79 7d 0a 49 6e 20 |nd{eqnar|ray}.In |
|00001ce0| 74 68 65 20 6d 69 64 64 | 6c 65 20 6c 69 6e 65 20 |the midd|le line |
|00001cf0| 6f 66 20 28 5c 72 65 66 | 7b 65 71 32 2e 31 31 7d |of (\ref|{eq2.11}|
|00001d00| 29 20 77 65 20 68 61 76 | 65 20 75 73 65 64 20 74 |) we hav|e used t|
|00001d10| 68 65 20 66 61 63 74 20 | 74 68 61 74 20 74 68 65 |he fact |that the|
|00001d20| 20 72 61 6e 67 65 20 6f | 66 0a 24 52 24 20 69 73 | range o|f.$R$ is|
|00001d30| 20 61 6c 77 61 79 73 20 | 70 65 72 70 65 6e 64 69 | always |perpendi|
|00001d40| 63 75 6c 61 72 20 74 6f | 20 74 68 65 20 72 61 6e |cular to| the ran|
|00001d50| 67 65 20 6f 66 20 24 50 | 24 2e 20 20 54 68 65 20 |ge of $P|$. The |
|00001d60| 67 72 61 64 69 65 6e 74 | 20 24 5c 6e 61 62 6c 61 |gradient| $\nabla|
|00001d70| 20 46 24 20 20 6f 66 0a | 24 46 24 2c 20 74 68 65 | F$ of.|$F$, the|
|00001d80| 72 65 66 6f 72 65 2c 20 | 20 6d 61 79 20 62 65 20 |refore, | may be |
|00001d90| 69 6e 74 65 72 70 72 65 | 74 65 64 20 61 73 20 74 |interpre|ted as t|
|00001da0| 68 65 0a 70 61 69 72 20 | 6f 66 20 6d 61 74 72 69 |he.pair |of matri|
|00001db0| 63 65 73 3a 0a 5c 62 65 | 67 69 6e 7b 65 71 75 61 |ces:.\be|gin{equa|
|00001dc0| 74 69 6f 6e 7d 0a 20 5c | 6e 61 62 6c 61 20 46 28 |tion}. \|nabla F(|
|00001dd0| 55 2c 56 29 20 3d 20 28 | 52 28 55 2c 56 29 56 5c |U,V) = (|R(U,V)V\|
|00001de0| 53 69 67 6d 61 5e 7b 54 | 7d 2c 52 28 55 2c 56 29 |Sigma^{T|},R(U,V)|
|00001df0| 5e 7b 54 7d 55 5c 53 69 | 67 6d 61 20 29 20 5c 69 |^{T}U\Si|gma ) \i|
|00001e00| 6e 0a 52 5e 7b 6d 20 5c | 74 69 6d 65 73 20 6d 7d |n.R^{m \|times m}|
|00001e10| 20 5c 74 69 6d 65 73 20 | 52 5e 7b 6e 20 5c 74 69 | \times |R^{n \ti|
|00001e20| 6d 65 73 20 6e 7d 2e 20 | 20 20 20 20 20 20 20 20 |mes n}. | |
|00001e30| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001e40| 20 20 5c 6c 61 62 65 6c | 7b 65 71 32 2e 31 32 7d | \label|{eq2.12}|
|00001e50| 0a 5c 65 6e 64 7b 65 71 | 75 61 74 69 6f 6e 7d 0a |.\end{eq|uation}.|
|00001e60| 42 65 63 61 75 73 65 20 | 6f 66 20 74 68 65 20 70 |Because |of the p|
|00001e70| 72 6f 64 75 63 74 20 74 | 6f 70 6f 6c 6f 67 79 2c |roduct t|opology,|
|00001e80| 20 77 65 20 6b 6e 6f 77 | 0a 5c 62 65 67 69 6e 7b | we know|.\begin{|
|00001e90| 65 71 75 61 74 69 6f 6e | 7d 0a 20 7b 5c 63 61 6c |equation|}. {\cal|
|00001ea0| 20 54 7d 5f 7b 28 55 2c | 56 29 7d 28 7b 5c 63 61 | T}_{(U,|V)}({\ca|
|00001eb0| 6c 20 4f 7d 20 28 6d 29 | 20 5c 74 69 6d 65 73 20 |l O} (m)| \times |
|00001ec0| 7b 5c 63 61 6c 20 4f 7d | 20 28 6e 29 29 20 3d 0a |{\cal O}| (n)) =.|
|00001ed0| 7b 5c 63 61 6c 20 54 7d | 5f 7b 55 7d 7b 5c 63 61 |{\cal T}|_{U}{\ca|
|00001ee0| 6c 20 4f 7d 20 28 6d 29 | 20 5c 74 69 6d 65 73 20 |l O} (m)| \times |
|00001ef0| 7b 5c 63 61 6c 20 54 7d | 5f 7b 56 7d 7b 5c 63 61 |{\cal T}|_{V}{\ca|
|00001f00| 6c 20 4f 7d 20 28 6e 29 | 2c 20 20 20 20 20 20 20 |l O} (n)|, |
|00001f10| 20 20 20 20 20 20 20 20 | 5c 6c 61 62 65 6c 7b 65 | |\label{e|
|00001f20| 71 32 2e 31 33 7d 0a 5c | 65 6e 64 7b 65 71 75 61 |q2.13}.\|end{equa|
|00001f30| 74 69 6f 6e 7d 0a 77 68 | 65 72 65 20 20 24 7b 5c |tion}.wh|ere ${\|
|00001f40| 63 61 6c 20 54 7d 5f 7b | 28 55 2c 56 29 7d 28 7b |cal T}_{|(U,V)}({|
|00001f50| 5c 63 61 6c 20 4f 7d 20 | 28 6d 29 20 5c 74 69 6d |\cal O} |(m) \tim|
|00001f60| 65 73 20 7b 5c 63 61 6c | 20 4f 7d 20 28 6e 29 29 |es {\cal| O} (n))|
|00001f70| 24 20 20 73 74 61 6e 64 | 73 20 66 6f 72 20 74 68 |$ stand|s for th|
|00001f80| 65 0a 74 61 6e 67 65 6e | 74 20 73 70 61 63 65 20 |e.tangen|t space |
|00001f90| 74 6f 20 74 68 65 20 6d | 61 6e 69 66 6f 6c 64 20 |to the m|anifold |
|00001fa0| 20 24 7b 5c 63 61 6c 20 | 4f 7d 20 28 6d 29 20 5c | ${\cal |O} (m) \|
|00001fb0| 74 69 6d 65 73 20 7b 5c | 63 61 6c 20 4f 7d 20 28 |times {\|cal O} (|
|00001fc0| 6e 29 24 20 20 61 74 20 | 20 24 28 55 2c 56 29 0a |n)$ at | $(U,V).|
|00001fd0| 5c 69 6e 20 7b 5c 63 61 | 6c 20 4f 7d 20 28 6d 29 |\in {\ca|l O} (m)|
|00001fe0| 20 5c 74 69 6d 65 73 20 | 7b 5c 63 61 6c 20 4f 7d | \times |{\cal O}|
|00001ff0| 20 28 6e 29 24 20 20 61 | 6e 64 20 73 6f 20 6f 6e | (n)$ a|nd so on|
|00002000| 2e 20 20 54 68 65 20 70 | 72 6f 6a 65 63 74 69 6f |. The p|rojectio|
|00002010| 6e 20 6f 66 0a 24 5c 6e | 61 62 6c 61 20 46 28 55 |n of.$\n|abla F(U|
|00002020| 2c 56 29 24 20 20 6f 6e | 74 6f 20 20 24 7b 5c 63 |,V)$ on|to ${\c|
|00002030| 61 6c 20 54 7d 5f 7b 28 | 55 2c 56 29 7d 28 7b 5c |al T}_{(|U,V)}({\|
|00002040| 63 61 6c 20 4f 7d 20 28 | 6d 29 20 5c 74 69 6d 65 |cal O} (|m) \time|
|00002050| 73 20 7b 5c 63 61 6c 20 | 4f 7d 20 28 6e 29 29 24 |s {\cal |O} (n))$|
|00002060| 2c 0a 74 68 65 72 65 66 | 6f 72 65 2c 20 69 73 20 |,.theref|ore, is |
|00002070| 74 68 65 20 70 72 6f 64 | 75 63 74 20 6f 66 20 74 |the prod|uct of t|
|00002080| 68 65 20 70 72 6f 6a 65 | 63 74 69 6f 6e 20 6f 66 |he proje|ction of|
|00002090| 20 74 68 65 20 66 69 72 | 73 74 20 63 6f 6d 70 6f | the fir|st compo|
|000020a0| 6e 65 6e 74 20 6f 66 0a | 24 5c 6e 61 62 6c 61 20 |nent of.|$\nabla |
|000020b0| 46 28 55 2c 56 29 24 20 | 20 6f 6e 74 6f 20 20 24 |F(U,V)$ | onto $|
|000020c0| 7b 5c 63 61 6c 20 54 7d | 5f 7b 55 7d 7b 5c 63 61 |{\cal T}|_{U}{\ca|
|000020d0| 6c 20 4f 7d 20 28 6d 29 | 24 20 20 61 6e 64 20 74 |l O} (m)|$ and t|
|000020e0| 68 65 20 70 72 6f 6a 65 | 63 74 69 6f 6e 20 6f 66 |he proje|ction of|
|000020f0| 20 74 68 65 0a 73 65 63 | 6f 6e 64 20 63 6f 6d 70 | the.sec|ond comp|
|00002100| 6f 6e 65 6e 74 20 6f 66 | 20 20 24 5c 6e 61 62 6c |onent of| $\nabl|
|00002110| 61 20 46 28 55 2c 56 29 | 24 20 20 6f 6e 74 6f 20 |a F(U,V)|$ onto |
|00002120| 20 24 7b 5c 63 61 6c 20 | 54 7d 5f 7b 56 7d 7b 5c | ${\cal |T}_{V}{\|
|00002130| 63 61 6c 20 4f 7d 20 28 | 6e 29 24 2e 20 0a 49 6e |cal O} (|n)$. .In|
|00002140| 20 70 61 72 74 69 63 75 | 6c 61 72 2c 20 77 65 20 | particu|lar, we |
|00002150| 63 6c 61 69 6d 20 74 68 | 61 74 20 74 68 65 0a 70 |claim th|at the.p|
|00002160| 72 6f 6a 65 63 74 69 6f | 6e 20 24 20 67 28 55 2c |rojectio|n $ g(U,|
|00002170| 56 29 24 20 20 6f 66 20 | 74 68 65 20 67 72 61 64 |V)$ of |the grad|
|00002180| 69 65 6e 74 20 20 24 5c | 6e 61 62 6c 61 20 46 28 |ient $\|nabla F(|
|00002190| 55 2c 56 29 24 20 20 6f | 6e 74 6f 0a 24 7b 5c 63 |U,V)$ o|nto.${\c|
|000021a0| 61 6c 20 54 7d 5f 7b 28 | 55 2c 56 29 7d 28 7b 5c |al T}_{(|U,V)}({\|
|000021b0| 63 61 6c 20 4f 7d 20 28 | 6d 29 20 5c 74 69 6d 65 |cal O} (|m) \time|
|000021c0| 73 20 7b 5c 63 61 6c 20 | 4f 7d 20 28 6e 29 29 24 |s {\cal |O} (n))$|
|000021d0| 20 20 69 73 20 67 69 76 | 65 6e 20 62 79 20 74 68 | is giv|en by th|
|000021e0| 65 20 70 61 69 72 20 6f | 66 0a 6d 61 74 72 69 63 |e pair o|f.matric|
|000021f0| 65 73 3a 0a 5c 62 65 67 | 69 6e 7b 65 71 6e 61 72 |es:.\beg|in{eqnar|
|00002200| 72 61 79 7d 0a 67 28 55 | 2c 56 29 20 3d 20 26 26 |ray}.g(U|,V) = &&|
|00002210| 20 5c 6c 65 66 74 28 20 | 5c 66 72 61 63 7b 52 28 | \left( |\frac{R(|
|00002220| 55 2c 56 29 56 5c 53 69 | 67 6d 61 5e 7b 54 7d 55 |U,V)V\Si|gma^{T}U|
|00002230| 5e 7b 54 7d 2d 55 5c 53 | 69 67 6d 61 20 56 5e 7b |^{T}-U\S|igma V^{|
|00002240| 54 7d 52 28 55 2c 56 29 | 5e 7b 54 7d 7d 7b 32 7d |T}R(U,V)|^{T}}{2}|
|00002250| 55 2c 0a 5c 72 69 67 68 | 74 2e 20 20 20 20 20 20 |U,.\righ|t. |
|00002260| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002270| 20 20 20 5c 6e 6f 6e 75 | 6d 62 65 72 20 5c 5c 5b | \nonu|mber \\[|
|00002280| 2d 31 2e 35 65 78 5d 0a | 5c 6c 61 62 65 6c 7b 65 |-1.5ex].|\label{e|
|00002290| 71 32 2e 31 34 7d 5c 5c | 5b 2d 31 2e 35 65 78 5d |q2.14}\\|[-1.5ex]|
|000022a0| 0a 26 26 5c 71 75 61 64 | 20 5c 6c 65 66 74 2e 20 |.&&\quad| \left. |
|000022b0| 5c 66 72 61 63 7b 52 28 | 55 2c 56 29 5e 7b 54 7d |\frac{R(|U,V)^{T}|
|000022c0| 55 5c 53 69 67 6d 61 20 | 56 5e 7b 54 7d 2d 56 0a |U\Sigma |V^{T}-V.|
|000022d0| 20 20 20 5c 53 69 67 6d | 61 5e 7b 54 7d 55 5e 7b | \Sigm|a^{T}U^{|
|000022e0| 54 7d 52 28 55 2c 56 29 | 7d 7b 32 7d 56 20 5c 72 |T}R(U,V)|}{2}V \r|
|000022f0| 69 67 68 74 29 2e 5c 6e | 6f 6e 75 6d 62 65 72 0a |ight).\n|onumber.|
|00002300| 5c 65 6e 64 7b 65 71 6e | 61 72 72 61 79 7d 0a 54 |\end{eqn|array}.T|
|00002310| 68 75 73 2c 20 74 68 65 | 20 76 65 63 74 6f 72 20 |hus, the| vector |
|00002320| 66 69 65 6c 64 0a 5c 62 | 65 67 69 6e 7b 65 71 75 |field.\b|egin{equ|
|00002330| 61 74 69 6f 6e 7d 0a 5c | 66 72 61 63 7b 64 28 55 |ation}.\|frac{d(U|
|00002340| 2c 56 29 7d 7b 64 74 7d | 20 3d 20 2d 67 28 55 2c |,V)}{dt}| = -g(U,|
|00002350| 56 29 20 20 20 20 20 5c | 6c 61 62 65 6c 7b 65 71 |V) \|label{eq|
|00002360| 32 2e 31 35 7d 0a 5c 65 | 6e 64 7b 65 71 75 61 74 |2.15}.\e|nd{equat|
|00002370| 69 6f 6e 7d 0a 64 65 66 | 69 6e 65 73 20 61 20 73 |ion}.def|ines a s|
|00002380| 74 65 65 70 65 73 74 20 | 64 65 73 63 65 6e 74 20 |teepest |descent |
|00002390| 66 6c 6f 77 20 6f 6e 20 | 74 68 65 20 6d 61 6e 69 |flow on |the mani|
|000023a0| 66 6f 6c 64 20 20 24 7b | 5c 63 61 6c 20 4f 7d 20 |fold ${|\cal O} |
|000023b0| 28 6d 29 20 5c 74 69 6d | 65 73 0a 7b 5c 63 61 6c |(m) \tim|es.{\cal|
|000023c0| 20 4f 7d 20 28 6e 29 24 | 20 66 6f 72 20 74 68 65 | O} (n)$| for the|
|000023d0| 20 6f 62 6a 65 63 74 69 | 76 65 20 66 75 6e 63 74 | objecti|ve funct|
|000023e0| 69 6f 6e 20 20 24 46 28 | 55 2c 56 29 24 2e 0a 0a |ion $F(|U,V)$...|
|000023f0| 20 0a 5c 73 65 63 74 69 | 6f 6e 7b 4d 61 69 6e 20 | .\secti|on{Main |
|00002400| 72 65 73 75 6c 74 73 7d | 20 0a 20 0a 4c 65 74 20 |results}| . .Let |
|00002410| 24 28 53 2c 43 29 24 20 | 62 65 20 61 20 6d 61 74 |$(S,C)$ |be a mat|
|00002420| 72 69 78 20 70 61 69 72 | 20 6f 66 20 6f 72 64 65 |rix pair| of orde|
|00002430| 72 20 24 6e 24 2e 20 20 | 54 68 65 20 64 65 74 65 |r $n$. |The dete|
|00002440| 72 6d 69 6e 61 6e 74 0a | 5c 5b 5c 64 65 74 20 28 |rminant.|\[\det (|
|00002450| 53 5c 63 69 72 63 20 58 | 2b 43 29 5c 5d 20 0a 69 |S\circ X|+C)\] .i|
|00002460| 73 20 61 20 70 6f 6c 79 | 6e 6f 6d 69 61 6c 20 69 |s a poly|nomial i|
|00002470| 6e 20 74 68 65 20 69 6e | 64 65 74 65 72 6d 69 6e |n the in|determin|
|00002480| 61 74 65 73 20 6f 66 20 | 24 58 24 20 6f 66 20 64 |ates of |$X$ of d|
|00002490| 65 67 72 65 65 20 61 74 | 0a 6d 6f 73 74 20 24 6e |egree at|.most $n|
|000024a0| 24 20 6f 76 65 72 20 74 | 68 65 20 72 65 61 6c 20 |$ over t|he real |
|000024b0| 66 69 65 6c 64 2e 20 57 | 65 20 63 61 6c 6c 20 74 |field. W|e call t|
|000024c0| 68 69 73 20 70 6f 6c 79 | 6e 6f 6d 69 61 6c 20 74 |his poly|nomial t|
|000024d0| 68 65 0a 7b 5c 65 6d 20 | 69 6e 64 69 63 61 74 6f |he.{\em |indicato|
|000024e0| 72 20 70 6f 6c 79 6e 6f | 6d 69 61 6c 7d 20 6f 66 |r polyno|mial} of|
|000024f0| 20 74 68 65 20 6d 61 74 | 72 69 78 20 70 61 69 72 | the mat|rix pair|
|00002500| 20 24 28 53 2c 43 29 24 | 0a 62 65 63 61 75 73 65 | $(S,C)$|.because|
|00002510| 20 6f 66 20 74 68 65 20 | 66 6f 6c 6c 6f 77 69 6e | of the |followin|
|00002520| 67 20 70 72 6f 70 6f 73 | 69 74 69 6f 6e 2e 0a 20 |g propos|ition.. |
|00002530| 0a 5c 62 65 67 69 6e 7b | 74 68 65 6f 72 65 6d 7d |.\begin{|theorem}|
|00002540| 20 0a 5c 6c 61 62 65 6c | 7b 74 68 3a 70 72 6f 70 | .\label|{th:prop|
|00002550| 7d 20 0a 54 68 65 20 6d | 61 74 72 69 78 20 70 61 |} .The m|atrix pa|
|00002560| 69 72 20 24 28 53 2c 43 | 29 24 20 69 73 20 61 20 |ir $(S,C|)$ is a |
|00002570| 7b 5c 72 6d 20 53 4e 53 | 7d 2d 6d 61 74 72 69 78 |{\rm SNS|}-matrix|
|00002580| 20 70 61 69 72 20 69 66 | 20 61 6e 64 0a 6f 6e 6c | pair if| and.onl|
|00002590| 79 20 69 66 20 61 6c 6c | 20 74 68 65 20 6e 6f 6e |y if all| the non|
|000025a0| 7a 65 72 6f 20 63 6f 65 | 66 66 69 63 69 65 6e 74 |zero coe|fficient|
|000025b0| 73 20 69 6e 20 69 74 73 | 20 69 6e 64 69 63 61 74 |s in its| indicat|
|000025c0| 6f 72 0a 70 6f 6c 79 6e | 6f 6d 69 61 6c 20 68 61 |or.polyn|omial ha|
|000025d0| 76 65 20 74 68 65 20 73 | 61 6d 65 20 73 69 67 6e |ve the s|ame sign|
|000025e0| 20 61 6e 64 20 74 68 65 | 72 65 20 69 73 20 61 74 | and the|re is at|
|000025f0| 20 6c 65 61 73 74 20 6f | 6e 65 0a 6e 6f 6e 7a 65 | least o|ne.nonze|
|00002600| 72 6f 20 63 6f 65 66 66 | 69 63 69 65 6e 74 2e 0a |ro coeff|icient..|
|00002610| 5c 65 6e 64 7b 74 68 65 | 6f 72 65 6d 7d 20 0a 20 |\end{the|orem} . |
|00002620| 0a 5c 62 65 67 69 6e 7b | 70 72 6f 6f 66 7d 0a 41 |.\begin{|proof}.A|
|00002630| 73 73 75 6d 65 20 74 68 | 61 74 20 24 28 53 2c 43 |ssume th|at $(S,C|
|00002640| 29 24 20 69 73 20 61 20 | 53 4e 53 2d 6d 61 74 72 |)$ is a |SNS-matr|
|00002650| 69 78 20 70 61 69 72 2e | 20 20 43 6c 65 61 72 6c |ix pair.| Clearl|
|00002660| 79 20 74 68 65 0a 69 6e | 64 69 63 61 74 6f 72 20 |y the.in|dicator |
|00002670| 70 6f 6c 79 6e 6f 6d 69 | 61 6c 20 68 61 73 20 61 |polynomi|al has a|
|00002680| 20 6e 6f 6e 7a 65 72 6f | 20 63 6f 65 66 66 69 63 | nonzero| coeffic|
|00002690| 69 65 6e 74 2e 20 20 43 | 6f 6e 73 69 64 65 72 20 |ient. C|onsider |
|000026a0| 61 0a 6d 6f 6e 6f 6d 69 | 61 6c 0a 5c 62 65 67 69 |a.monomi|al.\begi|
|000026b0| 6e 7b 65 71 75 61 74 69 | 6f 6e 7d 20 0a 5c 6c 61 |n{equati|on} .\la|
|000026c0| 62 65 6c 7b 65 71 3a 6d | 6f 6e 6f 7d 20 0a 62 5f |bel{eq:m|ono} .b_|
|000026d0| 7b 69 5f 7b 31 7d 2c 5c | 6c 64 6f 74 73 2c 69 5f |{i_{1},\|ldots,i_|
|000026e0| 7b 6b 7d 3b 6a 5f 7b 31 | 7d 2c 5c 6c 64 6f 74 73 |{k};j_{1|},\ldots|
|000026f0| 2c 6a 5f 7b 6b 7d 7d 78 | 5f 7b 69 5f 7b 31 7d 6a |,j_{k}}x|_{i_{1}j|
|00002700| 5f 7b 31 7d 7d 5c 63 64 | 6f 74 73 0a 78 5f 7b 69 |_{1}}\cd|ots.x_{i|
|00002710| 5f 7b 6b 7d 6a 5f 7b 6b | 7d 7d 0a 5c 65 6e 64 7b |_{k}j_{k|}}.\end{|
|00002720| 65 71 75 61 74 69 6f 6e | 7d 20 0a 6f 63 63 75 72 |equation|} .occur|
|00002730| 72 69 6e 67 20 69 6e 20 | 74 68 65 20 69 6e 64 69 |ring in |the indi|
|00002740| 63 61 74 6f 72 20 70 6f | 6c 79 6e 6f 6d 69 61 6c |cator po|lynomial|
|00002750| 20 77 69 74 68 20 61 20 | 6e 6f 6e 7a 65 72 6f 0a | with a |nonzero.|
|00002760| 63 6f 65 66 66 69 63 69 | 65 6e 74 2e 20 20 42 79 |coeffici|ent. By|
|00002770| 20 74 61 6b 69 6e 67 20 | 74 68 65 20 24 78 5f 7b | taking |the $x_{|
|00002780| 69 6a 7d 24 20 74 68 61 | 74 20 6f 63 63 75 72 20 |ij}$ tha|t occur |
|00002790| 69 6e 0a 28 5c 72 65 66 | 7b 65 71 3a 6d 6f 6e 6f |in.(\ref|{eq:mono|
|000027a0| 7d 29 20 6c 61 72 67 65 | 20 61 6e 64 20 61 6c 6c |}) large| and all|
|000027b0| 20 6f 74 68 65 72 73 20 | 73 6d 61 6c 6c 2c 20 77 | others |small, w|
|000027c0| 65 20 73 65 65 20 74 68 | 61 74 20 61 6e 79 0a 6d |e see th|at any.m|
|000027d0| 6f 6e 6f 6d 69 61 6c 20 | 74 68 61 74 20 6f 63 63 |onomial |that occ|
|000027e0| 75 72 73 20 69 6e 20 74 | 68 65 20 69 6e 64 69 63 |urs in t|he indic|
|000027f0| 61 74 6f 72 20 70 6f 6c | 79 6e 6f 6d 69 61 6c 20 |ator pol|ynomial |
|00002800| 77 69 74 68 20 61 0a 6e | 6f 6e 7a 65 72 6f 20 63 |with a.n|onzero c|
|00002810| 6f 65 66 66 69 63 69 65 | 6e 74 20 63 61 6e 20 62 |oefficie|nt can b|
|00002820| 65 20 6d 61 64 65 20 74 | 6f 20 64 6f 6d 69 6e 61 |e made t|o domina|
|00002830| 74 65 20 61 6c 6c 20 6f | 74 68 65 72 73 2e 0a 48 |te all o|thers..H|
|00002840| 65 6e 63 65 20 61 6c 6c | 20 74 68 65 20 6e 6f 6e |ence all| the non|
|00002850| 7a 65 72 6f 20 63 6f 65 | 66 66 69 63 69 65 6e 74 |zero coe|fficient|
|00002860| 73 20 68 61 76 65 20 74 | 68 65 20 73 61 6d 65 20 |s have t|he same |
|00002870| 73 69 67 6e 2e 20 54 68 | 65 0a 63 6f 6e 76 65 72 |sign. Th|e.conver|
|00002880| 73 65 20 69 73 20 69 6d | 2d 5c 6c 69 6e 65 62 72 |se is im|-\linebr|
|00002890| 65 61 6b 20 6d 65 64 69 | 61 74 65 2e 20 5c 71 71 |eak medi|ate. \qq|
|000028a0| 75 61 64 5c 65 6e 64 7b | 70 72 6f 6f 66 7d 20 0a |uad\end{|proof} .|
|000028b0| 20 0a 0a 46 6f 72 20 53 | 4e 53 2d 6d 61 74 72 69 | ..For S|NS-matri|
|000028c0| 78 20 70 61 69 72 73 20 | 24 28 53 2c 43 29 24 20 |x pairs |$(S,C)$ |
|000028d0| 77 69 74 68 20 24 43 3d | 4f 24 20 74 68 65 20 69 |with $C=|O$ the i|
|000028e0| 6e 64 69 63 61 74 6f 72 | 0a 70 6f 6c 79 6e 6f 6d |ndicator|.polynom|
|000028f0| 69 61 6c 20 69 73 20 61 | 20 68 6f 6d 6f 67 65 6e |ial is a| homogen|
|00002900| 65 6f 75 73 20 70 6f 6c | 79 6e 6f 6d 69 61 6c 20 |eous pol|ynomial |
|00002910| 6f 66 20 64 65 67 72 65 | 65 20 24 6e 24 2e 20 49 |of degre|e $n$. I|
|00002920| 6e 0a 74 68 69 73 20 63 | 61 73 65 20 54 68 65 6f |n.this c|ase Theo|
|00002930| 72 65 6d 20 5c 72 65 66 | 7b 74 68 3a 70 72 6f 70 |rem \ref|{th:prop|
|00002940| 7d 20 69 73 20 61 20 73 | 74 61 6e 64 61 72 64 20 |} is a s|tandard |
|00002950| 66 61 63 74 20 61 62 6f | 75 74 0a 53 4e 53 2d 6d |fact abo|ut.SNS-m|
|00002960| 61 74 72 69 63 65 73 2e | 0a 20 0a 5c 62 65 67 69 |atrices.|. .\begi|
|00002970| 6e 7b 6c 65 6d 6d 61 7d | 5b 7b 5c 72 6d 20 53 74 |n{lemma}|[{\rm St|
|00002980| 61 62 69 6c 69 74 79 7d | 5d 0a 5c 6c 61 62 65 6c |ability}|].\label|
|00002990| 7b 73 74 61 62 69 6c 69 | 74 79 7d 0a 47 69 76 65 |{stabili|ty}.Give|
|000029a0| 6e 20 24 54 3e 30 24 2c | 20 73 75 70 70 6f 73 65 |n $T>0$,| suppose|
|000029b0| 20 74 68 61 74 20 24 5c | 7c 20 5c 65 70 73 69 6c | that $\|| \epsil|
|000029c0| 6f 6e 20 28 74 29 20 5c | 7c 5f 7b 31 2c 32 7d 20 |on (t) \||_{1,2} |
|000029d0| 5c 6c 65 71 20 68 5e 7b | 71 2d 32 7d 24 0a 66 6f |\leq h^{|q-2}$.fo|
|000029e0| 72 20 24 30 20 5c 6c 65 | 71 20 74 20 5c 6c 65 71 |r $0 \le|q t \leq|
|000029f0| 20 54 24 20 61 6e 64 20 | 24 71 20 5c 67 65 71 20 | T$ and |$q \geq |
|00002a00| 36 24 2e 20 0a 54 68 65 | 6e 20 74 68 65 72 65 20 |6$. .The|n there |
|00002a10| 65 78 69 73 74 73 20 61 | 20 70 6f 73 69 74 69 76 |exists a| positiv|
|00002a20| 65 20 6e 75 6d 62 65 72 | 20 24 42 24 20 74 68 61 |e number| $B$ tha|
|00002a30| 74 20 64 65 70 65 6e 64 | 73 20 6f 6e 0a 24 54 24 |t depend|s on.$T$|
|00002a40| 20 61 6e 64 20 74 68 65 | 20 65 78 61 63 74 20 73 | and the| exact s|
|00002a50| 6f 6c 75 74 69 6f 6e 20 | 24 5c 70 65 24 20 6f 6e |olution |$\pe$ on|
|00002a60| 6c 79 20 73 75 63 68 20 | 74 68 61 74 20 66 6f 72 |ly such |that for|
|00002a70| 20 61 6c 6c 20 24 30 20 | 5c 6c 65 71 20 74 20 5c | all $0 |\leq t \|
|00002a80| 6c 65 71 20 54 24 2c 0a | 5c 62 65 67 69 6e 7b 65 |leq T$,.|\begin{e|
|00002a90| 71 75 61 74 69 6f 6e 7d | 0a 5c 6c 61 62 65 6c 7b |quation}|.\label{|
|00002aa0| 47 72 6f 6e 7d 0a 5c 66 | 72 61 63 20 7b 64 7d 7b |Gron}.\f|rac {d}{|
|00002ab0| 64 74 7d 20 5c 7c 20 5c | 65 70 73 69 6c 6f 6e 20 |dt} \| \|epsilon |
|00002ac0| 28 74 29 20 5c 7c 20 5f | 7b 31 2c 32 7d 20 20 5c |(t) \| _|{1,2} \|
|00002ad0| 6c 65 71 20 42 0a 20 20 | 20 28 20 68 5e 7b 71 2d |leq B. | ( h^{q-|
|00002ae0| 33 2f 32 7d 20 2b 20 5c | 7c 20 5c 65 70 73 69 6c |3/2} + \|| \epsil|
|00002af0| 6f 6e 20 28 74 29 20 5c | 7c 5f 7b 31 2c 32 7d 29 |on (t) \||_{1,2})|
|00002b00| 5c 3b 2e 0a 5c 65 6e 64 | 7b 65 71 75 61 74 69 6f |\;..\end|{equatio|
|00002b10| 6e 7d 0a 54 68 65 20 66 | 75 6e 63 74 69 6f 6e 20 |n}.The f|unction |
|00002b20| 24 42 28 54 29 24 20 63 | 61 6e 20 62 65 20 63 68 |$B(T)$ c|an be ch|
|00002b30| 6f 73 65 6e 20 74 6f 20 | 62 65 20 6e 6f 6e 64 65 |osen to |be nonde|
|00002b40| 63 72 65 61 73 69 6e 67 | 20 69 6e 20 74 69 6d 65 |creasing| in time|
|00002b50| 2e 0a 5c 65 6e 64 7b 6c | 65 6d 6d 61 7d 0a 20 0a |..\end{l|emma}. .|
|00002b60| 20 0a 5c 62 65 67 69 6e | 7b 74 68 65 6f 72 65 6d | .\begin|{theorem|
|00002b70| 7d 20 0a 5c 6c 61 62 65 | 6c 7b 74 68 3a 67 69 62 |} .\labe|l{th:gib|
|00002b80| 73 6f 6e 7d 20 0a 54 68 | 65 20 6d 61 78 69 6d 75 |son} .Th|e maximu|
|00002b90| 6d 20 6e 75 6d 62 65 72 | 20 6f 66 20 6e 6f 6e 7a |m number| of nonz|
|00002ba0| 65 72 6f 20 65 6e 74 72 | 69 65 73 20 69 6e 20 61 |ero entr|ies in a|
|00002bb0| 20 7b 5c 72 6d 20 53 4e | 53 7d 2d 6d 61 74 72 69 | {\rm SN|S}-matri|
|00002bc0| 78 0a 24 53 24 20 6f 66 | 20 6f 72 64 65 72 20 24 |x.$S$ of| order $|
|00002bd0| 6e 24 20 65 71 75 61 6c | 73 20 5c 5b 5c 66 72 61 |n$ equal|s \[\fra|
|00002be0| 63 7b 6e 5e 7b 32 7d 2b | 33 6e 2d 32 7d 7b 32 7d |c{n^{2}+|3n-2}{2}|
|00002bf0| 5c 5d 20 77 69 74 68 0a | 65 71 75 61 6c 69 74 79 |\] with.|equality|
|00002c00| 20 69 66 20 61 6e 64 20 | 6f 6e 6c 79 20 69 66 20 | if and |only if |
|00002c10| 74 68 65 72 65 20 65 78 | 69 73 74 20 70 65 72 6d |there ex|ist perm|
|00002c20| 75 74 61 74 69 6f 6e 20 | 6d 61 74 72 69 63 65 73 |utation |matrices|
|00002c30| 0a 73 75 63 68 20 74 68 | 61 74 20 24 50 7c 53 7c |.such th|at $P|S||
|00002c40| 51 3d 54 5f 7b 6e 7d 24 | 20 77 68 65 72 65 0a 5c |Q=T_{n}$| where.\|
|00002c50| 62 65 67 69 6e 7b 65 71 | 75 61 74 69 6f 6e 7d 20 |begin{eq|uation} |
|00002c60| 0a 5c 6c 61 62 65 6c 7b | 65 71 3a 67 69 62 73 6f |.\label{|eq:gibso|
|00002c70| 6e 7d 20 0a 54 5f 7b 6e | 7d 3d 5c 6c 65 66 74 5b |n} .T_{n|}=\left[|
|00002c80| 5c 62 65 67 69 6e 7b 61 | 72 72 61 79 7d 7b 63 63 |\begin{a|rray}{cc|
|00002c90| 63 63 63 63 7d 20 31 26 | 31 26 5c 63 64 6f 74 73 |cccc} 1&|1&\cdots|
|00002ca0| 26 31 26 31 26 31 5c 5c | 0a 31 26 31 26 5c 63 64 |&1&1&1\\|.1&1&\cd|
|00002cb0| 6f 74 73 26 31 26 31 26 | 31 5c 5c 20 30 26 31 26 |ots&1&1&|1\\ 0&1&|
|00002cc0| 5c 63 64 6f 74 73 26 31 | 26 31 26 31 5c 5c 20 0a |\cdots&1|&1&1\\ .|
|00002cd0| 5c 76 64 6f 74 73 26 5c | 76 64 6f 74 73 26 5c 64 |\vdots&\|vdots&\d|
|00002ce0| 64 6f 74 73 26 5c 76 64 | 6f 74 73 26 5c 76 64 6f |dots&\vd|ots&\vdo|
|00002cf0| 74 73 26 5c 76 64 6f 74 | 73 5c 5c 20 0a 30 26 30 |ts&\vdot|s\\ .0&0|
|00002d00| 26 5c 63 64 6f 74 73 26 | 31 26 31 26 31 5c 5c 20 |&\cdots&|1&1&1\\ |
|00002d10| 30 26 30 26 5c 63 64 6f | 74 73 26 30 26 31 26 31 |0&0&\cdo|ts&0&1&1|
|00002d20| 5c 65 6e 64 7b 61 72 72 | 61 79 7d 5c 72 69 67 68 |\end{arr|ay}\righ|
|00002d30| 74 5d 2e 20 0a 5c 65 6e | 64 7b 65 71 75 61 74 69 |t]. .\en|d{equati|
|00002d40| 6f 6e 7d 20 0a 5c 65 6e | 64 7b 74 68 65 6f 72 65 |on} .\en|d{theore|
|00002d50| 6d 7d 20 0a 20 0a 57 65 | 20 6e 6f 74 65 20 66 6f |m} . .We| note fo|
|00002d60| 72 20 6c 61 74 65 72 20 | 75 73 65 20 74 68 61 74 |r later |use that|
|00002d70| 20 65 61 63 68 20 73 75 | 62 6d 61 74 72 69 78 20 | each su|bmatrix |
|00002d80| 6f 66 20 24 54 5f 7b 6e | 7d 24 20 6f 66 0a 6f 72 |of $T_{n|}$ of.or|
|00002d90| 64 65 72 20 24 6e 2d 31 | 24 20 68 61 73 20 61 6c |der $n-1|$ has al|
|00002da0| 6c 20 31 73 20 6f 6e 20 | 69 74 73 20 6d 61 69 6e |l 1s on |its main|
|00002db0| 20 64 69 61 67 6f 6e 61 | 6c 2e 20 0a 20 0a 57 65 | diagona|l. . .We|
|00002dc0| 20 6e 6f 77 20 6f 62 74 | 61 69 6e 20 61 20 62 6f | now obt|ain a bo|
|00002dd0| 75 6e 64 20 6f 6e 20 74 | 68 65 20 6e 75 6d 62 65 |und on t|he numbe|
|00002de0| 72 20 6f 66 20 6e 6f 6e | 7a 65 72 6f 20 65 6e 74 |r of non|zero ent|
|00002df0| 72 69 65 73 20 6f 66 0a | 24 53 24 20 69 6e 20 61 |ries of.|$S$ in a|
|00002e00| 20 53 4e 53 2d 6d 61 74 | 72 69 78 20 70 61 69 72 | SNS-mat|rix pair|
|00002e10| 20 24 28 53 2c 43 29 24 | 20 69 6e 20 74 65 72 6d | $(S,C)$| in term|
|00002e20| 73 20 6f 66 20 74 68 65 | 20 64 65 67 72 65 65 20 |s of the| degree |
|00002e30| 6f 66 0a 74 68 65 20 69 | 6e 64 69 63 61 74 6f 72 |of.the i|ndicator|
|00002e40| 20 70 6f 6c 79 6e 6f 6d | 69 61 6c 2e 20 57 65 20 | polynom|ial. We |
|00002e50| 64 65 6e 6f 74 65 20 74 | 68 65 20 73 74 72 69 63 |denote t|he stric|
|00002e60| 74 6c 79 20 75 70 70 65 | 72 0a 74 72 69 61 6e 67 |tly uppe|r.triang|
|00002e70| 75 6c 61 72 20 28 30 2c | 31 29 2d 6d 61 74 72 69 |ular (0,|1)-matri|
|00002e80| 78 20 6f 66 20 6f 72 64 | 65 72 20 24 6d 24 20 77 |x of ord|er $m$ w|
|00002e90| 69 74 68 20 61 6c 6c 20 | 31 73 20 61 62 6f 76 65 |ith all |1s above|
|00002ea0| 20 74 68 65 0a 6d 61 69 | 6e 20 64 69 61 67 6f 6e | the.mai|n diagon|
|00002eb0| 61 6c 20 62 79 20 24 55 | 5f 7b 6d 7d 24 2e 20 54 |al by $U|_{m}$. T|
|00002ec0| 68 65 20 61 6c 6c 20 31 | 73 20 6d 61 74 72 69 78 |he all 1|s matrix|
|00002ed0| 20 6f 66 20 73 69 7a 65 | 20 24 6d 24 20 62 79 0a | of size| $m$ by.|
|00002ee0| 24 70 24 20 69 73 20 64 | 65 6e 6f 74 65 64 20 62 |$p$ is d|enoted b|
|00002ef0| 79 20 24 4a 5f 7b 6d 2c | 70 7d 24 2e 0a 20 0a 20 |y $J_{m,|p}$.. . |
|00002f00| 0a 5c 62 65 67 69 6e 7b | 70 72 6f 70 6f 73 69 74 |.\begin{|proposit|
|00002f10| 69 6f 6e 7d 5b 7b 5c 72 | 6d 20 43 6f 6e 76 6f 6c |ion}[{\r|m Convol|
|00002f20| 75 74 69 6f 6e 20 74 68 | 65 6f 72 65 6d 7d 5d 0a |ution th|eorem}].|
|00002f30| 5c 6c 61 62 65 6c 7b 70 | 72 6f 3a 32 2e 31 7d 20 |\label{p|ro:2.1} |
|00002f40| 20 4c 65 74 0a 5c 62 65 | 67 69 6e 7b 65 71 6e 61 | Let.\be|gin{eqna|
|00002f50| 72 72 61 79 2a 7d 0a 61 | 5c 61 73 74 20 75 28 74 |rray*}.a|\ast u(t|
|00002f60| 29 20 3d 20 5c 69 6e 74 | 5f 30 5e 74 20 61 28 74 |) = \int|_0^t a(t|
|00002f70| 2d 20 5c 74 61 75 29 20 | 75 28 5c 74 61 75 29 20 |- \tau) |u(\tau) |
|00002f80| 64 5c 74 61 75 2c 20 5c | 68 73 70 61 63 65 7b 2e |d\tau, \|hspace{.|
|00002f90| 32 69 6e 7d 20 74 20 5c | 69 6e 0a 28 30 2c 20 5c |2in} t \|in.(0, \|
|00002fa0| 69 6e 66 74 79 29 2e 0a | 5c 65 6e 64 7b 65 71 6e |infty)..|\end{eqn|
|00002fb0| 61 72 72 61 79 2a 7d 0a | 54 68 65 6e 0a 5c 62 65 |array*}.|Then.\be|
|00002fc0| 67 69 6e 7b 65 71 6e 61 | 72 72 61 79 2a 7d 0a 5c |gin{eqna|rray*}.\|
|00002fd0| 77 69 64 65 68 61 74 7b | 61 5c 61 73 74 20 75 7d |widehat{|a\ast u}|
|00002fe0| 28 73 29 20 3d 20 5c 77 | 69 64 65 68 61 74 7b 61 |(s) = \w|idehat{a|
|00002ff0| 7d 28 73 29 5c 77 69 64 | 65 68 61 74 7b 75 7d 28 |}(s)\wid|ehat{u}(|
|00003000| 73 29 2e 0a 5c 65 6e 64 | 7b 65 71 6e 61 72 72 61 |s)..\end|{eqnarra|
|00003010| 79 2a 7d 0a 5c 65 6e 64 | 7b 70 72 6f 70 6f 73 69 |y*}.\end|{proposi|
|00003020| 74 69 6f 6e 7d 0a 0a 5c | 62 65 67 69 6e 7b 6c 65 |tion}..\|begin{le|
|00003030| 6d 6d 61 7d 0a 5c 6c 61 | 62 65 6c 7b 6c 65 6d 3a |mma}.\la|bel{lem:|
|00003040| 33 2e 31 7d 0a 46 6f 72 | 20 24 73 5f 30 20 3e 30 |3.1}.For| $s_0 >0|
|00003050| 24 2c 20 69 66 0a 24 24 | 0a 5c 69 6e 74 5f 30 5e |$, if.$$|.\int_0^|
|00003060| 7b 5c 69 6e 66 74 79 7d | 20 65 5e 7b 2d 32 73 5f |{\infty}| e^{-2s_|
|00003070| 30 20 74 7d 76 5e 7b 28 | 31 29 7d 28 74 29 20 76 |0 t}v^{(|1)}(t) v|
|00003080| 28 74 29 20 64 74 20 5c | 3b 20 5c 6c 65 71 20 30 |(t) dt \|; \leq 0|
|00003090| 20 5c 3b 2c 0a 24 24 0a | 74 68 65 6e 0a 5c 62 65 | \;,.$$.|then.\be|
|000030a0| 67 69 6e 7b 65 71 6e 61 | 72 72 61 79 2a 7d 0a 5c |gin{eqna|rray*}.\|
|000030b0| 69 6e 74 5f 30 5e 7b 5c | 69 6e 66 74 79 7d 20 65 |int_0^{\|infty} e|
|000030c0| 5e 7b 2d 32 73 5f 30 20 | 74 7d 20 76 5e 32 28 74 |^{-2s_0 |t} v^2(t|
|000030d0| 29 20 64 74 20 5c 3b 20 | 5c 6c 65 71 20 5c 3b 20 |) dt \; |\leq \; |
|000030e0| 5c 66 72 61 63 7b 31 7d | 7b 32 73 5f 30 7d 20 76 |\frac{1}|{2s_0} v|
|000030f0| 5e 32 28 30 29 2e 0a 5c | 65 6e 64 7b 65 71 6e 61 |^2(0)..\|end{eqna|
|00003100| 72 72 61 79 2a 7d 0a 5c | 65 6e 64 7b 6c 65 6d 6d |rray*}.\|end{lemm|
|00003110| 61 7d 0a 0a 7b 5c 65 6d | 20 50 72 6f 6f 66 7d 2e |a}..{\em| Proof}.|
|00003120| 20 41 70 70 6c 79 69 6e | 67 20 69 6e 74 65 67 72 | Applyin|g integr|
|00003130| 61 74 69 6f 6e 20 62 79 | 20 70 61 72 74 73 2c 20 |ation by| parts, |
|00003140| 77 65 20 6f 62 74 61 69 | 6e 0a 5c 62 65 67 69 6e |we obtai|n.\begin|
|00003150| 7b 65 71 6e 61 72 72 61 | 79 2a 7d 0a 5c 69 6e 74 |{eqnarra|y*}.\int|
|00003160| 5f 30 5e 7b 5c 69 6e 66 | 74 79 7d 20 65 5e 7b 2d |_0^{\inf|ty} e^{-|
|00003170| 32 73 5f 30 20 74 7d 20 | 5b 76 5e 32 28 74 29 2d |2s_0 t} |[v^2(t)-|
|00003180| 76 5e 32 28 30 29 5d 20 | 64 74 0a 26 3d 26 5c 6c |v^2(0)] |dt.&=&\l|
|00003190| 69 6d 5f 7b 74 5c 72 69 | 67 68 74 61 72 72 6f 77 |im_{t\ri|ghtarrow|
|000031a0| 20 5c 69 6e 66 74 79 7d | 5c 6c 65 66 74 20 28 0a | \infty}|\left (.|
|000031b0| 2d 5c 66 72 61 63 7b 31 | 7d 7b 32 73 5f 30 7d 65 |-\frac{1|}{2s_0}e|
|000031c0| 5e 7b 2d 32 73 5f 30 20 | 74 7d 76 5e 32 28 74 29 |^{-2s_0 |t}v^2(t)|
|000031d0| 20 5c 72 69 67 68 74 20 | 29 20 2b 5c 66 72 61 63 | \right |) +\frac|
|000031e0| 7b 31 7d 7b 73 5f 30 7d | 0a 5c 69 6e 74 5f 30 5e |{1}{s_0}|.\int_0^|
|000031f0| 7b 5c 69 6e 66 74 79 7d | 20 65 5e 7b 2d 32 73 5f |{\infty}| e^{-2s_|
|00003200| 30 20 74 7d 76 5e 7b 28 | 31 29 7d 28 74 29 76 28 |0 t}v^{(|1)}(t)v(|
|00003210| 74 29 64 74 5c 5c 0a 26 | 5c 6c 65 71 26 20 5c 66 |t)dt\\.&|\leq& \f|
|00003220| 72 61 63 7b 31 7d 7b 73 | 5f 30 7d 20 5c 69 6e 74 |rac{1}{s|_0} \int|
|00003230| 5f 30 5e 7b 5c 69 6e 66 | 74 79 7d 20 65 5e 7b 2d |_0^{\inf|ty} e^{-|
|00003240| 32 73 5f 30 20 74 7d 20 | 76 5e 7b 28 31 29 7d 28 |2s_0 t} |v^{(1)}(|
|00003250| 74 29 76 28 74 29 20 64 | 74 20 5c 3b 5c 3b 0a 5c |t)v(t) d|t \;\;.\|
|00003260| 6c 65 71 20 5c 3b 5c 3b | 20 30 2e 0a 5c 65 6e 64 |leq \;\;| 0..\end|
|00003270| 7b 65 71 6e 61 72 72 61 | 79 2a 7d 0a 54 68 75 73 |{eqnarra|y*}.Thus|
|00003280| 0a 24 24 0a 5c 69 6e 74 | 5f 30 5e 7b 5c 69 6e 66 |.$$.\int|_0^{\inf|
|00003290| 74 79 7d 20 65 5e 7b 2d | 32 73 5f 30 20 74 7d 20 |ty} e^{-|2s_0 t} |
|000032a0| 76 5e 32 28 74 29 20 64 | 74 20 5c 3b 5c 3b 5c 6c |v^2(t) d|t \;\;\l|
|000032b0| 65 71 20 76 5e 32 28 30 | 29 20 5c 69 6e 74 5f 30 |eq v^2(0|) \int_0|
|000032c0| 5e 7b 5c 69 6e 66 74 79 | 7d 0a 5c 3b 5c 3b 65 5e |^{\infty|}.\;\;e^|
|000032d0| 7b 2d 32 73 5f 30 20 74 | 7d 20 64 74 5c 3b 5c 3b |{-2s_0 t|} dt\;\;|
|000032e0| 3d 5c 3b 5c 3b 5c 66 72 | 61 63 7b 31 7d 7b 32 73 |=\;\;\fr|ac{1}{2s|
|000032f0| 5f 30 7d 20 76 5e 32 28 | 30 29 2e 5c 65 71 6e 6f |_0} v^2(|0).\eqno|
|00003300| 5c 65 6e 64 70 72 6f 6f | 66 0a 24 24 0a 0a 5c 62 |\endproo|f.$$..\b|
|00003310| 65 67 69 6e 7b 63 6f 72 | 6f 6c 6c 61 72 79 7d 5c |egin{cor|ollary}\|
|00003320| 6c 61 62 65 6c 7b 63 34 | 2e 31 7d 0a 4c 65 74 20 |label{c4|.1}.Let |
|00003330| 24 20 5c 62 66 45 20 24 | 20 73 61 74 69 73 66 79 |$ \bfE $| satisfy|
|00003340| 20 24 28 35 29 24 2d 2d | 24 28 36 29 24 20 61 6e | $(5)$--|$(6)$ an|
|00003350| 64 0a 73 75 70 70 6f 73 | 65 20 24 20 5c 62 66 45 |d.suppos|e $ \bfE|
|00003360| 5e 68 20 24 20 73 61 74 | 69 73 66 69 65 73 20 24 |^h $ sat|isfies $|
|00003370| 28 37 29 24 20 61 6e 64 | 20 24 28 38 29 24 0a 77 |(7)$ and| $(8)$.w|
|00003380| 69 74 68 20 61 20 67 65 | 6e 65 72 61 6c 20 24 20 |ith a ge|neral $ |
|00003390| 5c 62 66 47 20 24 2e 20 | 20 4c 65 74 20 24 20 5c |\bfG $. | Let $ \|
|000033a0| 62 66 47 3d 20 5c 6e 61 | 62 6c 61 20 5c 74 69 6d |bfG= \na|bla \tim|
|000033b0| 65 73 20 7b 5c 62 66 20 | 5c 50 68 69 7d 20 2b 20 |es {\bf |\Phi} + |
|000033c0| 5c 6e 61 62 6c 61 20 70 | 2c 24 0a 24 70 20 5c 69 |\nabla p|,$.$p \i|
|000033d0| 6e 20 48 5f 30 5e 31 20 | 28 5c 4f 6d 65 67 61 29 |n H_0^1 |(\Omega)|
|000033e0| 20 24 2e 20 53 75 70 70 | 6f 73 65 20 74 68 61 74 | $. Supp|ose that|
|000033f0| 20 24 5c 6e 61 62 6c 61 | 20 70 24 20 61 6e 64 20 | $\nabla| p$ and |
|00003400| 24 20 5c 6e 61 62 6c 61 | 20 5c 74 69 6d 65 73 20 |$ \nabla| \times |
|00003410| 0a 7b 5c 62 66 20 5c 50 | 68 69 7d 20 24 20 73 61 |.{\bf \P|hi} $ sa|
|00003420| 74 69 73 66 79 20 61 6c | 6c 20 74 68 65 20 61 73 |tisfy al|l the as|
|00003430| 73 75 6d 70 74 69 6f 6e | 73 20 6f 66 20 54 68 65 |sumption|s of The|
|00003440| 6f 72 65 6d 73 20 24 34 | 2e 31 24 20 61 6e 64 20 |orems $4|.1$ and |
|00003450| 20 0a 24 34 2e 32 24 2c | 20 72 65 73 70 65 63 74 | .$4.2$,| respect|
|00003460| 69 76 65 6c 79 2e 20 49 | 6e 20 61 64 64 69 74 69 |ively. I|n additi|
|00003470| 6f 6e 20 73 75 70 70 6f | 73 65 20 61 6c 6c 20 74 |on suppo|se all t|
|00003480| 68 65 20 72 65 67 75 6c | 61 72 69 74 79 0a 61 73 |he regul|arity.as|
|00003490| 73 75 6d 70 74 69 6f 6e | 73 20 6f 66 20 54 68 65 |sumption|s of The|
|000034a0| 6f 72 65 6d 73 20 24 34 | 2e 31 24 2d 2d 24 34 2e |orems $4|.1$--$4.|
|000034b0| 32 24 20 61 72 65 20 73 | 61 74 69 73 66 69 65 64 |2$ are s|atisfied|
|000034c0| 2e 20 20 54 68 65 6e 20 | 0a 66 6f 72 20 24 20 30 |. Then |.for $ 0|
|000034d0| 20 5c 6c 65 20 74 20 5c | 6c 65 20 54 20 24 20 61 | \le t \|le T $ a|
|000034e0| 6e 64 20 24 20 30 20 3c | 20 5c 65 70 73 69 6c 6f |nd $ 0 <| \epsilo|
|000034f0| 6e 20 5c 6c 65 20 5c 65 | 70 73 69 6c 6f 6e 5f 30 |n \le \e|psilon_0|
|00003500| 20 24 20 74 68 65 72 65 | 20 65 78 69 73 74 73 20 | $ there| exists |
|00003510| 61 20 0a 63 6f 6e 73 74 | 61 6e 74 20 24 20 43 20 |a .const|ant $ C |
|00003520| 3d 20 43 28 5c 65 70 73 | 69 6c 6f 6e 2c 20 54 29 |= C(\eps|ilon, T)|
|00003530| 20 24 20 73 75 63 68 20 | 74 68 61 74 0a 24 24 0a | $ such |that.$$.|
|00003540| 5c 56 65 72 74 20 28 5c | 62 66 45 20 2d 20 5c 62 |\Vert (\|bfE - \b|
|00003550| 66 45 5e 68 29 28 74 29 | 20 5c 56 65 72 74 5f 30 |fE^h)(t)| \Vert_0|
|00003560| 20 5c 6c 65 20 43 20 68 | 5e 7b 6b 2b 31 2d 20 5c | \le C h|^{k+1- \|
|00003570| 65 70 73 69 6c 6f 6e 7d | 2c 0a 24 24 0a 77 68 65 |epsilon}|,.$$.whe|
|00003580| 72 65 20 24 20 43 20 24 | 20 61 6c 73 6f 20 64 65 |re $ C $| also de|
|00003590| 70 65 6e 64 73 20 6f 6e | 20 74 68 65 20 63 6f 6e |pends on| the con|
|000035a0| 73 74 61 6e 74 73 20 67 | 69 76 65 6e 20 69 6e 20 |stants g|iven in |
|000035b0| 54 68 65 6f 72 65 6d 73 | 20 0a 24 34 2e 31 24 20 |Theorems| .$4.1$ |
|000035c0| 61 6e 64 20 24 34 2e 32 | 24 2e 0a 5c 65 6e 64 7b |and $4.2|$..\end{|
|000035d0| 63 6f 72 6f 6c 6c 61 72 | 79 7d 0a 0a 0a 5c 62 65 |corollar|y}...\be|
|000035e0| 67 69 6e 7b 64 65 66 69 | 6e 69 74 69 6f 6e 7d 0a |gin{defi|nition}.|
|000035f0| 7b 5c 72 6d 20 4c 65 74 | 20 24 53 24 20 62 65 20 |{\rm Let| $S$ be |
|00003600| 61 6e 20 69 73 6f 6c 61 | 74 65 64 20 69 6e 76 61 |an isola|ted inva|
|00003610| 72 69 61 6e 74 20 73 65 | 74 20 77 69 74 68 20 69 |riant se|t with i|
|00003620| 73 6f 6c 61 74 69 6e 67 | 20 6e 65 69 67 68 62 6f |solating| neighbo|
|00003630| 72 68 6f 6f 64 20 24 4e | 24 2e 0a 41 6e 20 7b 5c |rhood $N|$..An {\|
|00003640| 65 6d 20 69 6e 64 65 78 | 20 70 61 69 72 7d 20 66 |em index| pair} f|
|00003650| 6f 72 20 24 53 24 20 69 | 73 20 61 20 70 61 69 72 |or $S$ i|s a pair|
|00003660| 20 6f 66 20 63 6f 6d 70 | 61 63 74 20 73 65 74 73 | of comp|act sets|
|00003670| 20 24 28 4e 5f 7b 31 7d | 2c 4e 5f 7b 30 7d 29 24 | $(N_{1}|,N_{0})$|
|00003680| 0a 77 69 74 68 20 24 4e | 5f 7b 30 7d 20 5c 73 75 |.with $N|_{0} \su|
|00003690| 62 73 65 74 20 4e 5f 7b | 31 7d 20 5c 73 75 62 73 |bset N_{|1} \subs|
|000036a0| 65 74 20 4e 24 20 73 75 | 63 68 20 74 68 61 74 3a |et N$ su|ch that:|
|000036b0| 0a 5c 62 65 67 69 6e 7b | 72 6f 6d 61 6e 6e 75 6d |.\begin{|romannum|
|000036c0| 7d 0a 5c 69 74 65 6d 20 | 24 63 6c 28 4e 5f 7b 31 |}.\item |$cl(N_{1|
|000036d0| 7d 20 5c 62 61 63 6b 73 | 6c 61 73 68 20 4e 5f 7b |} \backs|lash N_{|
|000036e0| 30 7d 29 24 0a 69 73 20 | 61 6e 20 69 73 6f 6c 61 |0})$.is |an isola|
|000036f0| 74 69 6e 67 20 6e 65 69 | 67 68 62 6f 72 68 6f 6f |ting nei|ghborhoo|
|00003700| 64 20 66 6f 72 20 24 53 | 24 2e 0a 5c 69 74 65 6d |d for $S|$..\item|
|00003710| 20 24 4e 5f 7b 69 7d 24 | 20 69 73 20 70 6f 73 69 | $N_{i}$| is posi|
|00003720| 74 69 76 65 6c 79 20 69 | 6e 76 61 72 69 61 6e 74 |tively i|nvariant|
|00003730| 20 72 65 6c 61 74 69 76 | 65 20 74 6f 20 24 4e 24 | relativ|e to $N$|
|00003740| 20 66 6f 72 20 24 69 3d | 30 2c 31 24 2c 0a 69 2e | for $i=|0,1$,.i.|
|00003750| 65 2e 2c 20 67 69 76 65 | 6e 0a 24 78 20 5c 69 6e |e., give|n.$x \in|
|00003760| 20 4e 5f 7b 69 7d 24 20 | 61 6e 64 20 24 78 20 5c | N_{i}$ |and $x \|
|00003770| 63 64 6f 74 20 5b 30 2c | 74 5d 20 5c 73 75 62 73 |cdot [0,|t] \subs|
|00003780| 65 74 20 4e 24 2c 20 74 | 68 65 6e 20 24 78 20 5c |et N$, t|hen $x \|
|00003790| 63 64 6f 74 20 5b 30 2c | 74 5d 20 5c 73 75 62 73 |cdot [0,|t] \subs|
|000037a0| 65 74 0a 4e 5f 7b 69 7d | 24 2e 0a 5c 69 74 65 6d |et.N_{i}|$..\item|
|000037b0| 20 24 4e 5f 7b 30 7d 24 | 20 69 73 20 61 6e 20 65 | $N_{0}$| is an e|
|000037c0| 78 69 74 20 73 65 74 20 | 66 6f 72 20 24 4e 5f 7b |xit set |for $N_{|
|000037d0| 31 7d 24 2c 20 69 2e 65 | 2e 20 69 66 20 24 78 20 |1}$, i.e|. if $x |
|000037e0| 5c 69 6e 20 4e 5f 7b 31 | 7d 24 2c 0a 24 78 20 5c |\in N_{1|}$,.$x \|
|000037f0| 63 64 6f 74 20 5b 30 2c | 20 5c 69 6e 66 74 79 20 |cdot [0,| \infty |
|00003800| 29 20 5c 6e 6f 74 5c 73 | 75 62 73 65 74 20 4e 5f |) \not\s|ubset N_|
|00003810| 7b 31 7d 24 2c 20 74 68 | 65 6e 20 74 68 65 72 65 |{1}$, th|en there|
|00003820| 20 69 73 20 61 20 24 54 | 20 5c 67 65 71 20 30 24 | is a $T| \geq 0$|
|00003830| 20 73 75 63 68 0a 74 68 | 61 74 20 24 78 20 5c 63 | such.th|at $x \c|
|00003840| 64 6f 74 20 5b 30 2c 54 | 5d 20 5c 73 75 62 73 65 |dot [0,T|] \subse|
|00003850| 74 20 4e 5f 7b 31 7d 24 | 20 61 6e 64 20 24 78 20 |t N_{1}$| and $x |
|00003860| 5c 63 64 6f 74 20 54 20 | 5c 69 6e 20 4e 5f 7b 30 |\cdot T |\in N_{0|
|00003870| 7d 24 2e 0a 5c 65 6e 64 | 7b 72 6f 6d 61 6e 6e 75 |}$..\end|{romannu|
|00003880| 6d 7d 7d 0a 5c 65 6e 64 | 7b 64 65 66 69 6e 69 74 |m}}.\end|{definit|
|00003890| 69 6f 6e 7d 0a 0a 5c 73 | 75 62 73 65 63 74 69 6f |ion}..\s|ubsectio|
|000038a0| 6e 7b 4e 75 6d 65 72 69 | 63 61 6c 20 65 78 70 65 |n{Numeri|cal expe|
|000038b0| 72 69 6d 65 6e 74 73 7d | 20 57 65 20 63 6f 6e 64 |riments}| We cond|
|000038c0| 75 63 74 65 64 20 6e 75 | 6d 65 72 69 63 61 6c 20 |ucted nu|merical |
|000038d0| 65 78 70 65 72 69 6d 65 | 6e 74 73 20 0a 69 6e 20 |experime|nts .in |
|000038e0| 63 6f 6d 70 75 74 69 6e | 67 20 69 6e 65 78 61 63 |computin|g inexac|
|000038f0| 74 20 4e 65 77 74 6f 6e | 20 73 74 65 70 73 20 66 |t Newton| steps f|
|00003900| 6f 72 20 64 69 73 63 72 | 65 74 69 7a 61 74 69 6f |or discr|etizatio|
|00003910| 6e 73 20 6f 66 20 61 20 | 20 0a 7b 5c 65 6d 20 6d |ns of a | .{\em m|
|00003920| 6f 64 69 66 69 65 64 20 | 42 72 61 74 75 20 70 72 |odified |Bratu pr|
|00003930| 6f 62 6c 65 6d 7d 2c 20 | 67 69 76 65 6e 20 62 79 |oblem}, |given by|
|00003940| 20 20 0a 5c 62 65 67 69 | 6e 7b 65 71 6e 61 72 72 | .\begi|n{eqnarr|
|00003950| 61 79 7d 20 0a 7b 5c 64 | 73 20 5c 44 65 6c 74 61 |ay} .{\d|s \Delta|
|00003960| 20 77 20 2b 20 63 20 65 | 5e 77 20 2b 20 64 7b 20 | w + c e|^w + d{ |
|00003970| 7b 5c 70 61 72 74 69 61 | 6c 20 77 7d 5c 6f 76 65 |{\partia|l w}\ove|
|00003980| 72 7b 5c 70 61 72 74 69 | 61 6c 20 78 7d 20 7d 20 |r{\parti|al x} } |
|00003990| 7d 20 0a 26 3d 26 7b 5c | 64 73 20 66 20 5c 71 75 |} .&=&{\|ds f \qu|
|000039a0| 61 64 20 7b 5c 72 6d 20 | 69 6e 7d 5c 20 44 2c 20 |ad {\rm |in}\ D, |
|000039b0| 7d 5c 6e 6f 6e 75 6d 62 | 65 72 5c 5c 5b 2d 31 2e |}\nonumb|er\\[-1.|
|000039c0| 35 65 78 5d 0a 5c 6c 61 | 62 65 6c 7b 62 72 61 74 |5ex].\la|bel{brat|
|000039d0| 75 7d 20 5c 5c 5b 2d 31 | 2e 35 65 78 5d 0a 7b 5c |u} \\[-1|.5ex].{\|
|000039e0| 64 73 20 77 20 7d 26 3d | 26 7b 5c 64 73 20 30 20 |ds w }&=|&{\ds 0 |
|000039f0| 5c 71 75 61 64 20 7b 5c | 72 6d 20 6f 6e 7d 5c 20 |\quad {\|rm on}\ |
|00003a00| 5c 70 61 72 74 69 61 6c | 20 44 20 2c 20 7d 20 5c |\partial| D , } \|
|00003a10| 6e 6f 6e 75 6d 62 65 72 | 0a 5c 65 6e 64 7b 65 71 |nonumber|.\end{eq|
|00003a20| 6e 61 72 72 61 79 7d 20 | 0a 77 68 65 72 65 20 24 |narray} |.where $|
|00003a30| 63 24 20 61 6e 64 20 24 | 64 24 20 61 72 65 20 63 |c$ and $|d$ are c|
|00003a40| 6f 6e 73 74 61 6e 74 73 | 2e 20 54 68 65 20 61 63 |onstants|. The ac|
|00003a50| 74 75 61 6c 20 42 72 61 | 74 75 20 70 72 6f 62 6c |tual Bra|tu probl|
|00003a60| 65 6d 20 68 61 73 20 24 | 64 3d 30 24 20 61 6e 64 |em has $|d=0$ and|
|00003a70| 20 20 0a 24 66 20 5c 65 | 71 75 69 76 30 24 2e 20 | .$f \e|quiv0$. |
|00003a80| 49 74 20 70 72 6f 76 69 | 64 65 73 20 61 20 73 69 |It provi|des a si|
|00003a90| 6d 70 6c 69 66 69 65 64 | 20 6d 6f 64 65 6c 20 6f |mplified| model o|
|00003aa0| 66 20 6e 6f 6e 6c 69 6e | 65 61 72 20 64 69 66 66 |f nonlin|ear diff|
|00003ab0| 75 73 69 6f 6e 20 20 0a | 70 68 65 6e 6f 6d 65 6e |usion .|phenomen|
|00003ac0| 61 2c 20 65 2e 67 2e 2c | 20 69 6e 20 63 6f 6d 62 |a, e.g.,| in comb|
|00003ad0| 75 73 74 69 6f 6e 20 61 | 6e 64 20 73 65 6d 69 63 |ustion a|nd semic|
|00003ae0| 6f 6e 64 75 63 74 6f 72 | 73 2c 20 61 6e 64 20 68 |onductor|s, and h|
|00003af0| 61 73 20 62 65 65 6e 20 | 0a 63 6f 6e 73 69 64 65 |as been |.conside|
|00003b00| 72 65 64 20 62 79 20 47 | 6c 6f 77 69 6e 73 6b 69 |red by G|lowinski|
|00003b10| 2c 20 4b 65 6c 6c 65 72 | 2c 20 61 6e 64 20 52 68 |, Keller|, and Rh|
|00003b20| 65 69 6e 68 61 72 64 74 | 20 5c 63 69 74 65 7b 47 |einhardt| \cite{G|
|00003b30| 6c 6f 4b 52 38 35 7d 2c | 20 0a 61 73 20 77 65 6c |loKR85},| .as wel|
|00003b40| 6c 20 61 73 20 62 79 20 | 61 20 6e 75 6d 62 65 72 |l as by |a number|
|00003b50| 20 6f 66 20 6f 74 68 65 | 72 20 69 6e 76 65 73 74 | of othe|r invest|
|00003b60| 69 67 61 74 6f 72 73 3b | 20 73 65 65 20 5c 63 69 |igators;| see \ci|
|00003b70| 74 65 7b 47 6c 6f 4b 52 | 38 35 7d 20 0a 61 6e 64 |te{GloKR|85} .and|
|00003b80| 20 74 68 65 20 72 65 66 | 65 72 65 6e 63 65 73 20 | the ref|erences |
|00003b90| 74 68 65 72 65 69 6e 2e | 20 53 65 65 20 61 6c 73 |therein.| See als|
|00003ba0| 6f 20 70 72 6f 62 6c 65 | 6d 20 33 20 62 79 20 47 |o proble|m 3 by G|
|00003bb0| 6c 6f 77 69 6e 73 6b 69 | 20 61 6e 64 20 20 4b 65 |lowinski| and Ke|
|00003bc0| 6c 6c 65 72 20 20 0a 61 | 6e 64 20 70 72 6f 62 6c |ller .a|nd probl|
|00003bd0| 65 6d 20 37 20 62 79 20 | 4d 69 74 74 65 6c 6d 61 |em 7 by |Mittelma|
|00003be0| 6e 6e 20 69 6e 20 74 68 | 65 20 63 6f 6c 6c 65 63 |nn in th|e collec|
|00003bf0| 74 69 6f 6e 20 6f 66 20 | 6e 6f 6e 6c 69 6e 65 61 |tion of |nonlinea|
|00003c00| 72 20 6d 6f 64 65 6c 20 | 0a 70 72 6f 62 6c 65 6d |r model |.problem|
|00003c10| 73 20 61 73 73 65 6d 62 | 6c 65 64 20 62 79 20 4d |s assemb|led by M|
|00003c20| 6f 72 5c 27 65 20 5c 63 | 69 74 65 7b 4d 6f 72 65 |or\'e \c|ite{More|
|00003c30| 7d 2e 20 54 68 65 20 6d | 6f 64 69 66 69 65 64 20 |}. The m|odified |
|00003c40| 70 72 6f 62 6c 65 6d 20 | 20 0a 28 5c 72 65 66 7b |problem | .(\ref{|
|00003c50| 62 72 61 74 75 7d 29 20 | 68 61 73 20 62 65 65 6e |bratu}) |has been|
|00003c60| 20 75 73 65 64 20 61 73 | 20 61 20 74 65 73 74 20 | used as| a test |
|00003c70| 70 72 6f 62 6c 65 6d 20 | 66 6f 72 20 69 6e 65 78 |problem |for inex|
|00003c80| 61 63 74 20 4e 65 77 74 | 6f 6e 20 0a 6d 65 74 68 |act Newt|on .meth|
|00003c90| 6f 64 73 20 62 79 20 42 | 72 6f 77 6e 20 61 6e 64 |ods by B|rown and|
|00003ca0| 20 53 61 61 64 20 5c 63 | 69 74 65 7b 42 72 6f 77 | Saad \c|ite{Brow|
|00003cb0| 6e 2d 53 61 61 64 31 7d | 2e 20 20 0a 20 0a 49 6e |n-Saad1}|. . .In|
|00003cc0| 20 6f 75 72 20 65 78 70 | 65 72 69 6d 65 6e 74 73 | our exp|eriments|
|00003cd0| 2c 20 77 65 20 74 6f 6f | 6b 20 24 44 20 3d 20 5b |, we too|k $D = [|
|00003ce0| 30 2c 31 5d 5c 74 69 6d | 65 73 5b 30 2c 31 5d 24 |0,1]\tim|es[0,1]$|
|00003cf0| 2c 20 24 66 20 5c 65 71 | 75 69 76 30 24 2c 20 0a |, $f \eq|uiv0$, .|
|00003d00| 24 63 3d 64 3d 31 30 24 | 2c 20 61 6e 64 20 64 69 |$c=d=10$|, and di|
|00003d10| 73 63 72 65 74 69 7a 65 | 64 20 28 5c 72 65 66 7b |scretize|d (\ref{|
|00003d20| 62 72 61 74 75 7d 29 20 | 75 73 69 6e 67 20 74 68 |bratu}) |using th|
|00003d30| 65 20 75 73 75 61 6c 20 | 73 65 63 6f 6e 64 2d 6f |e usual |second-o|
|00003d40| 72 64 65 72 20 0a 63 65 | 6e 74 65 72 65 64 20 64 |rder .ce|ntered d|
|00003d50| 69 66 66 65 72 65 6e 63 | 65 73 20 6f 76 65 72 20 |ifferenc|es over |
|00003d60| 61 20 24 31 30 30 5c 74 | 69 6d 65 73 31 30 30 24 |a $100\t|imes100$|
|00003d70| 20 6d 65 73 68 20 6f 66 | 20 65 71 75 61 6c 6c 79 | mesh of| equally|
|00003d80| 20 0a 73 70 61 63 65 64 | 20 70 6f 69 6e 74 73 20 | .spaced| points |
|00003d90| 69 6e 20 24 44 24 2e 20 | 49 6e 20 5c 67 6d 72 65 |in $D$. |In \gmre|
|00003da0| 73 28 24 6d 24 29 2c 20 | 77 65 20 74 6f 6f 6b 20 |s($m$), |we took |
|00003db0| 24 6d 3d 31 30 24 20 61 | 6e 64 20 75 73 65 64 20 |$m=10$ a|nd used |
|00003dc0| 66 61 73 74 20 20 0a 50 | 6f 69 73 73 6f 6e 20 72 |fast .P|oisson r|
|00003dd0| 69 67 68 74 20 70 72 65 | 63 6f 6e 64 69 74 69 6f |ight pre|conditio|
|00003de0| 6e 69 6e 67 20 61 73 20 | 69 6e 20 74 68 65 20 65 |ning as |in the e|
|00003df0| 78 70 65 72 69 6d 65 6e | 74 73 20 69 6e 20 5c 53 |xperimen|ts in \S|
|00003e00| 32 2e 20 54 68 65 20 63 | 6f 6d 70 75 74 69 6e 67 |2. The c|omputing|
|00003e10| 20 20 0a 65 6e 76 69 72 | 6f 6e 6d 65 6e 74 20 77 | .envir|onment w|
|00003e20| 61 73 20 61 73 20 64 65 | 73 63 72 69 62 65 64 20 |as as de|scribed |
|00003e30| 69 6e 20 5c 53 32 2e 20 | 41 6c 6c 20 63 6f 6d 70 |in \S2. |All comp|
|00003e40| 75 74 69 6e 67 20 77 61 | 73 20 64 6f 6e 65 20 20 |uting wa|s done |
|00003e50| 0a 69 6e 20 64 6f 75 62 | 6c 65 20 70 72 65 63 69 |.in doub|le preci|
|00003e60| 73 69 6f 6e 2e 20 20 0a | 20 0a 20 0a 0a 5c 62 65 |sion. .| . ..\be|
|00003e70| 67 69 6e 7b 66 69 67 75 | 72 65 7d 5b 68 74 5d 20 |gin{figu|re}[ht] |
|00003e80| 0a 5c 76 73 70 61 63 65 | 7b 32 2e 35 69 6e 7d 0a |.\vspace|{2.5in}.|
|00003e90| 5c 63 61 70 74 69 6f 6e | 7b 7b 5c 72 6d 20 4c 6f |\caption|{{\rm Lo|
|00003ea0| 67 7d 24 5f 7b 31 30 7d | 24 20 6f 66 20 74 68 65 |g}$_{10}|$ of the|
|00003eb0| 20 72 65 73 69 64 75 61 | 6c 20 6e 6f 72 6d 20 76 | residua|l norm v|
|00003ec0| 65 72 73 75 73 20 74 68 | 65 20 6e 75 6d 62 65 72 |ersus th|e number|
|00003ed0| 20 6f 66 0a 7b 5c 72 6d | 20 47 4d 52 45 53 24 28 | of.{\rm| GMRES$(|
|00003ee0| 6d 29 24 7d 20 69 74 65 | 72 61 74 69 6f 6e 73 20 |m)$} ite|rations |
|00003ef0| 66 6f 72 20 74 68 65 20 | 66 69 6e 69 74 65 20 64 |for the |finite d|
|00003f00| 69 66 66 65 72 65 6e 63 | 65 20 6d 65 74 68 6f 64 |ifferenc|e method|
|00003f10| 73 2e 7d 20 0a 5c 6c 61 | 62 65 6c 7b 64 69 66 66 |s.} .\la|bel{diff|
|00003f20| 7d 20 0a 5c 65 6e 64 7b | 66 69 67 75 72 65 7d 20 |} .\end{|figure} |
|00003f30| 20 0a 20 0a 0a 20 0a 49 | 6e 20 74 68 65 20 66 69 | . .. .I|n the fi|
|00003f40| 72 73 74 20 73 65 74 20 | 6f 66 20 65 78 70 65 72 |rst set |of exper|
|00003f50| 69 6d 65 6e 74 73 2c 20 | 77 65 20 61 6c 6c 6f 77 |iments, |we allow|
|00003f60| 65 64 20 65 61 63 68 20 | 6d 65 74 68 6f 64 20 74 |ed each |method t|
|00003f70| 6f 20 20 0a 72 75 6e 20 | 66 6f 72 20 24 34 30 24 |o .run |for $40$|
|00003f80| 20 7b 5c 67 6d 72 65 73 | 6d 7d 20 69 74 65 72 61 | {\gmres|m} itera|
|00003f90| 74 69 6f 6e 73 2c 20 73 | 74 61 72 74 69 6e 67 20 |tions, s|tarting |
|00003fa0| 77 69 74 68 20 7a 65 72 | 6f 20 61 73 20 74 68 65 |with zer|o as the|
|00003fb0| 20 69 6e 69 74 69 61 6c | 20 20 0a 61 70 70 72 6f | initial| .appro|
|00003fc0| 78 69 6d 61 74 65 20 73 | 6f 6c 75 74 69 6f 6e 2c |ximate s|olution,|
|00003fd0| 20 61 66 74 65 72 20 77 | 68 69 63 68 20 74 68 65 | after w|hich the|
|00003fe0| 20 6c 69 6d 69 74 20 6f | 66 20 72 65 73 69 64 75 | limit o|f residu|
|00003ff0| 61 6c 20 6e 6f 72 6d 20 | 20 0a 72 65 64 75 63 74 |al norm | .reduct|
|00004000| 69 6f 6e 20 68 61 64 20 | 62 65 65 6e 20 72 65 61 |ion had |been rea|
|00004010| 63 68 65 64 2e 20 54 68 | 65 20 72 65 73 75 6c 74 |ched. Th|e result|
|00004020| 73 20 61 72 65 20 73 68 | 6f 77 6e 20 69 6e 20 46 |s are sh|own in F|
|00004030| 69 67 2e 7e 5c 72 65 66 | 7b 64 69 66 66 7d 2e 20 |ig.~\ref|{diff}. |
|00004040| 20 0a 49 6e 20 46 69 67 | 2e 7e 5c 72 65 66 7b 64 | .In Fig|.~\ref{d|
|00004050| 69 66 66 7d 2c 20 74 68 | 65 20 74 6f 70 20 63 75 |iff}, th|e top cu|
|00004060| 72 76 65 20 77 61 73 20 | 70 72 6f 64 75 63 65 64 |rve was |produced|
|00004070| 20 62 79 20 6d 65 74 68 | 6f 64 20 46 44 31 2e 20 | by meth|od FD1. |
|00004080| 0a 54 68 65 20 73 65 63 | 6f 6e 64 20 63 75 72 76 |.The sec|ond curv|
|00004090| 65 20 66 72 6f 6d 20 74 | 68 65 20 74 6f 70 20 69 |e from t|he top i|
|000040a0| 73 20 61 63 74 75 61 6c | 6c 79 20 61 20 73 75 70 |s actual|ly a sup|
|000040b0| 65 72 70 6f 73 69 74 69 | 6f 6e 20 6f 66 20 20 0a |erpositi|on of .|
|000040c0| 74 68 65 20 63 75 72 76 | 65 73 20 70 72 6f 64 75 |the curv|es produ|
|000040d0| 63 65 64 20 62 79 20 6d | 65 74 68 6f 64 73 20 45 |ced by m|ethods E|
|000040e0| 48 41 32 20 61 6e 64 20 | 46 44 32 3b 20 74 68 65 |HA2 and |FD2; the|
|000040f0| 20 74 77 6f 20 63 75 72 | 76 65 73 20 61 72 65 20 | two cur|ves are |
|00004100| 0a 76 69 73 75 61 6c 6c | 79 20 69 6e 64 69 73 74 |.visuall|y indist|
|00004110| 69 6e 67 75 69 73 68 61 | 62 6c 65 2e 20 53 69 6d |inguisha|ble. Sim|
|00004120| 69 6c 61 72 6c 79 2c 20 | 74 68 65 20 74 68 69 72 |ilarly, |the thir|
|00004130| 64 20 63 75 72 76 65 20 | 66 72 6f 6d 20 20 0a 74 |d curve |from .t|
|00004140| 68 65 20 74 6f 70 20 69 | 73 20 61 20 73 75 70 65 |he top i|s a supe|
|00004150| 72 70 6f 73 69 74 69 6f | 6e 20 6f 66 20 74 68 65 |rpositio|n of the|
|00004160| 20 63 75 72 76 65 73 20 | 70 72 6f 64 75 63 65 64 | curves |produced|
|00004170| 20 62 79 20 6d 65 74 68 | 6f 64 73 20 45 48 41 34 | by meth|ods EHA4|
|00004180| 20 0a 61 6e 64 20 46 44 | 34 2c 20 61 6e 64 20 74 | .and FD|4, and t|
|00004190| 68 65 20 66 6f 75 72 74 | 68 20 63 75 72 76 65 20 |he fourt|h curve |
|000041a0| 66 72 6f 6d 20 74 68 65 | 20 74 6f 70 2c 20 77 68 |from the| top, wh|
|000041b0| 69 63 68 20 6c 69 65 73 | 20 62 61 72 65 6c 79 20 |ich lies| barely |
|000041c0| 61 62 6f 76 65 20 20 0a | 74 68 65 20 62 6f 74 74 |above .|the bott|
|000041d0| 6f 6d 20 63 75 72 76 65 | 2c 20 69 73 20 61 20 73 |om curve|, is a s|
|000041e0| 75 70 65 72 70 6f 73 69 | 74 69 6f 6e 20 6f 66 20 |uperposi|tion of |
|000041f0| 74 68 65 20 63 75 72 76 | 65 73 20 70 72 6f 64 75 |the curv|es produ|
|00004200| 63 65 64 20 62 79 20 20 | 0a 6d 65 74 68 6f 64 73 |ced by |.methods|
|00004210| 20 45 48 41 36 20 61 6e | 64 20 46 44 36 2e 20 54 | EHA6 an|d FD6. T|
|00004220| 68 65 20 62 6f 74 74 6f | 6d 20 63 75 72 76 65 20 |he botto|m curve |
|00004230| 77 61 73 20 70 72 6f 64 | 75 63 65 64 20 62 79 20 |was prod|uced by |
|00004240| 6d 65 74 68 6f 64 20 41 | 2e 20 0a 20 0a 49 6e 20 |method A|. . .In |
|00004250| 74 68 65 20 73 65 63 6f | 6e 64 20 73 65 74 20 6f |the seco|nd set o|
|00004260| 66 20 65 78 70 65 72 69 | 6d 65 6e 74 73 2c 20 6f |f experi|ments, o|
|00004270| 75 72 20 70 75 72 70 6f | 73 65 20 77 61 73 20 74 |ur purpo|se was t|
|00004280| 6f 20 61 73 73 65 73 73 | 20 74 68 65 20 0a 72 65 |o assess| the .re|
|00004290| 6c 61 74 69 76 65 20 61 | 6d 6f 75 6e 74 20 6f 66 |lative a|mount of|
|000042a0| 20 63 6f 6d 70 75 74 61 | 74 69 6f 6e 61 6c 20 77 | computa|tional w|
|000042b0| 6f 72 6b 20 72 65 71 75 | 69 72 65 64 20 62 79 20 |ork requ|ired by |
|000042c0| 74 68 65 20 6d 65 74 68 | 6f 64 73 20 20 0a 77 68 |the meth|ods .wh|
|000042d0| 69 63 68 20 75 73 65 20 | 68 69 67 68 65 72 2d 6f |ich use |higher-o|
|000042e0| 72 64 65 72 20 64 69 66 | 66 65 72 65 6e 63 69 6e |rder dif|ferencin|
|000042f0| 67 20 74 6f 20 72 65 61 | 63 68 20 63 6f 6d 70 61 |g to rea|ch compa|
|00004300| 72 61 62 6c 65 20 6c 65 | 76 65 6c 73 20 20 0a 6f |rable le|vels .o|
|00004310| 66 20 72 65 73 69 64 75 | 61 6c 20 6e 6f 72 6d 20 |f residu|al norm |
|00004320| 72 65 64 75 63 74 69 6f | 6e 2e 20 57 65 20 63 6f |reductio|n. We co|
|00004330| 6d 70 61 72 65 64 20 70 | 61 69 72 73 20 6f 66 20 |mpared p|airs of |
|00004340| 6d 65 74 68 6f 64 73 20 | 45 48 41 32 20 20 0a 61 |methods |EHA2 .a|
|00004350| 6e 64 20 46 44 32 2c 20 | 45 48 41 34 20 61 6e 64 |nd FD2, |EHA4 and|
|00004360| 20 46 44 34 2c 20 61 6e | 64 20 45 48 41 36 20 61 | FD4, an|d EHA6 a|
|00004370| 6e 64 20 46 44 36 20 62 | 79 20 6f 62 73 65 72 76 |nd FD6 b|y observ|
|00004380| 69 6e 67 20 69 6e 20 65 | 61 63 68 20 6f 66 20 0a |ing in e|ach of .|
|00004390| 32 30 20 74 72 69 61 6c | 73 20 74 68 65 20 6e 75 |20 trial|s the nu|
|000043a0| 6d 62 65 72 20 6f 66 20 | 7b 5c 67 6d 72 65 73 6d |mber of |{\gmresm|
|000043b0| 7d 20 69 74 65 72 61 74 | 69 6f 6e 73 2c 20 6e 75 |} iterat|ions, nu|
|000043c0| 6d 62 65 72 20 6f 66 20 | 24 46 24 2d 65 76 61 6c |mber of |$F$-eval|
|000043d0| 75 61 74 69 6f 6e 73 2c | 20 20 0a 61 6e 64 20 72 |uations,| .and r|
|000043e0| 75 6e 20 74 69 6d 65 20 | 72 65 71 75 69 72 65 64 |un time |required|
|000043f0| 20 62 79 20 65 61 63 68 | 20 6d 65 74 68 6f 64 20 | by each| method |
|00004400| 74 6f 20 72 65 64 75 63 | 65 20 74 68 65 20 72 65 |to reduc|e the re|
|00004410| 73 69 64 75 61 6c 20 6e | 6f 72 6d 20 0a 62 79 20 |sidual n|orm .by |
|00004420| 61 20 66 61 63 74 6f 72 | 20 6f 66 20 24 5c 65 24 |a factor| of $\e$|
|00004430| 2c 20 77 68 65 72 65 20 | 66 6f 72 20 65 61 63 68 |, where |for each|
|00004440| 20 70 61 69 72 20 6f 66 | 20 6d 65 74 68 6f 64 73 | pair of| methods|
|00004450| 20 24 5c 65 24 20 77 61 | 73 20 63 68 6f 73 65 6e | $\e$ wa|s chosen|
|00004460| 20 20 0a 74 6f 20 62 65 | 20 73 6f 6d 65 77 68 61 | .to be| somewha|
|00004470| 74 20 67 72 65 61 74 65 | 72 20 74 68 61 6e 20 74 |t greate|r than t|
|00004480| 68 65 20 6c 69 6d 69 74 | 69 6e 67 20 72 61 74 69 |he limit|ing rati|
|00004490| 6f 20 6f 66 20 66 69 6e | 61 6c 20 74 6f 20 20 0a |o of fin|al to .|
|000044a0| 69 6e 69 74 69 61 6c 20 | 72 65 73 69 64 75 61 6c |initial |residual|
|000044b0| 20 6e 6f 72 6d 73 20 6f | 62 74 61 69 6e 61 62 6c | norms o|btainabl|
|000044c0| 65 20 62 79 20 74 68 65 | 20 6d 65 74 68 6f 64 73 |e by the| methods|
|000044d0| 2e 20 49 6e 20 74 68 65 | 73 65 20 74 72 69 61 6c |. In the|se trial|
|000044e0| 73 2c 20 20 0a 74 68 65 | 20 69 6e 69 74 69 61 6c |s, .the| initial|
|000044f0| 20 61 70 70 72 6f 78 69 | 6d 61 74 65 20 73 6f 6c | approxi|mate sol|
|00004500| 75 74 69 6f 6e 73 20 77 | 65 72 65 20 6f 62 74 61 |utions w|ere obta|
|00004510| 69 6e 65 64 20 62 79 20 | 67 65 6e 65 72 61 74 69 |ined by |generati|
|00004520| 6e 67 20 72 61 6e 64 6f | 6d 20 20 0a 63 6f 6d 70 |ng rando|m .comp|
|00004530| 6f 6e 65 6e 74 73 20 61 | 73 20 69 6e 20 74 68 65 |onents a|s in the|
|00004540| 20 73 69 6d 69 6c 61 72 | 20 65 78 70 65 72 69 6d | similar| experim|
|00004550| 65 6e 74 73 20 69 6e 20 | 5c 53 32 2e 20 57 65 20 |ents in |\S2. We |
|00004560| 6e 6f 74 65 20 74 68 61 | 74 20 66 6f 72 20 65 76 |note tha|t for ev|
|00004570| 65 72 79 20 20 0a 6d 65 | 74 68 6f 64 2c 20 74 68 |ery .me|thod, th|
|00004580| 65 20 6e 75 6d 62 65 72 | 73 20 6f 66 20 7b 5c 67 |e number|s of {\g|
|00004590| 6d 72 65 73 6d 7d 20 69 | 74 65 72 61 74 69 6f 6e |mresm} i|teration|
|000045a0| 73 20 61 6e 64 20 24 46 | 24 2d 65 76 61 6c 75 61 |s and $F|$-evalua|
|000045b0| 74 69 6f 6e 73 20 72 65 | 71 75 69 72 65 64 20 20 |tions re|quired |
|000045c0| 0a 62 65 66 6f 72 65 20 | 74 65 72 6d 69 6e 61 74 |.before |terminat|
|000045d0| 69 6f 6e 20 64 69 64 20 | 6e 6f 74 20 76 61 72 79 |ion did |not vary|
|000045e0| 20 61 74 20 61 6c 6c 20 | 6f 76 65 72 20 74 68 65 | at all |over the|
|000045f0| 20 32 30 20 74 72 69 61 | 6c 73 2e 20 54 68 65 20 | 20 tria|ls. The |
|00004600| 7b 5c 67 6d 72 65 73 6d | 7d 0a 69 74 65 72 61 74 |{\gmresm|}.iterat|
|00004610| 69 6f 6e 20 63 6f 75 6e | 74 73 2c 20 6e 75 6d 62 |ion coun|ts, numb|
|00004620| 65 72 73 20 6f 66 20 24 | 46 24 2d 65 76 61 6c 75 |ers of $|F$-evalu|
|00004630| 61 74 69 6f 6e 73 2c 20 | 61 6e 64 20 6d 65 61 6e |ations, |and mean|
|00004640| 73 20 61 6e 64 20 73 74 | 61 6e 64 61 72 64 20 20 |s and st|andard |
|00004650| 0a 64 65 76 69 61 74 69 | 6f 6e 73 20 6f 66 20 74 |.deviati|ons of t|
|00004660| 68 65 20 72 75 6e 20 74 | 69 6d 65 73 20 61 72 65 |he run t|imes are|
|00004670| 20 67 69 76 65 6e 20 69 | 6e 20 54 61 62 6c 65 20 | given i|n Table |
|00004680| 5c 72 65 66 7b 64 69 66 | 66 73 74 61 74 73 7d 2e |\ref{dif|fstats}.|
|00004690| 20 0a 20 0a 5c 62 65 67 | 69 6e 7b 74 61 62 6c 65 | . .\beg|in{table|
|000046a0| 7d 0a 5c 63 61 70 74 69 | 6f 6e 7b 53 74 61 74 69 |}.\capti|on{Stati|
|000046b0| 73 74 69 63 73 20 6f 76 | 65 72 20 24 32 30 24 20 |stics ov|er $20$ |
|000046c0| 74 72 69 61 6c 73 20 6f | 66 20 7b 5c 72 6d 20 47 |trials o|f {\rm G|
|000046d0| 4d 52 45 53 24 28 6d 29 | 24 7d 20 69 74 65 72 61 |MRES$(m)|$} itera|
|000046e0| 74 69 6f 6e 20 6e 75 6d | 62 65 72 73 2c 20 20 0a |tion num|bers, .|
|000046f0| 24 46 24 2d 65 76 61 6c | 75 61 74 69 6f 6e 73 2c |$F$-eval|uations,|
|00004700| 20 61 6e 64 20 72 75 6e | 20 74 69 6d 65 73 20 72 | and run| times r|
|00004710| 65 71 75 69 72 65 64 20 | 74 6f 20 72 65 64 75 63 |equired |to reduc|
|00004720| 65 20 74 68 65 20 72 65 | 73 69 64 75 61 6c 20 6e |e the re|sidual n|
|00004730| 6f 72 6d 20 62 79 20 20 | 0a 61 20 66 61 63 74 6f |orm by |.a facto|
|00004740| 72 20 6f 66 20 24 5c 65 | 24 2e 20 46 6f 72 20 65 |r of $\e|$. For e|
|00004750| 61 63 68 20 6d 65 74 68 | 6f 64 2c 20 74 68 65 20 |ach meth|od, the |
|00004760| 6e 75 6d 62 65 72 20 6f | 66 20 7b 5c 72 6d 20 47 |number o|f {\rm G|
|00004770| 4d 52 45 53 24 28 6d 29 | 24 7d 20 69 74 65 72 61 |MRES$(m)|$} itera|
|00004780| 74 69 6f 6e 73 20 20 0a | 61 6e 64 20 24 46 24 2d |tions .|and $F$-|
|00004790| 65 76 61 6c 75 61 74 69 | 6f 6e 73 20 77 61 73 20 |evaluati|ons was |
|000047a0| 74 68 65 20 73 61 6d 65 | 20 69 6e 20 65 76 65 72 |the same| in ever|
|000047b0| 79 20 74 72 69 61 6c 2e | 7d 20 0a 0a 5c 62 65 67 |y trial.|} ..\beg|
|000047c0| 69 6e 7b 63 65 6e 74 65 | 72 7d 20 7b 5c 66 6f 6f |in{cente|r} {\foo|
|000047d0| 74 6e 6f 74 65 73 69 7a | 65 0a 5c 62 65 67 69 6e |tnotesiz|e.\begin|
|000047e0| 7b 74 61 62 75 6c 61 72 | 7d 7b 7c 63 7c 63 7c 63 |{tabular|}{|c|c|c|
|000047f0| 7c 63 7c 63 7c 63 7c 7d | 20 5c 68 6c 69 6e 65 20 ||c|c|c|}| \hline |
|00004800| 20 0a 26 26 20 4e 75 6d | 62 65 72 20 6f 66 20 26 | .&& Num|ber of &|
|00004810| 20 4e 75 6d 62 65 72 20 | 6f 66 20 26 20 4d 65 61 | Number |of & Mea|
|00004820| 6e 20 52 75 6e 20 54 69 | 6d 65 20 26 20 53 74 61 |n Run Ti|me & Sta|
|00004830| 6e 64 61 72 64 20 5c 5c | 20 0a 4d 65 74 68 6f 64 |ndard \\| .Method|
|00004840| 20 26 20 24 5c 65 24 20 | 26 20 49 74 65 72 61 74 | & $\e$ |& Iterat|
|00004850| 69 6f 6e 73 20 26 20 24 | 46 24 2d 45 76 61 6c 75 |ions & $|F$-Evalu|
|00004860| 61 74 69 6f 6e 73 26 20 | 28 53 65 63 6f 6e 64 73 |ations& |(Seconds|
|00004870| 29 20 26 20 44 65 76 69 | 61 74 69 6f 6e 20 5c 5c |) & Devi|ation \\|
|00004880| 20 5c 68 6c 69 6e 65 20 | 0a 5c 6c 6f 77 65 72 2e | \hline |.\lower.|
|00004890| 33 65 78 5c 68 62 6f 78 | 7b 45 48 41 32 7d 20 26 |3ex\hbox|{EHA2} &|
|000048a0| 20 5c 6c 6f 77 65 72 2e | 33 65 78 5c 68 62 6f 78 | \lower.|3ex\hbox|
|000048b0| 7b 24 31 30 5e 7b 2d 31 | 30 7d 24 7d 20 26 20 5c |{$10^{-1|0}$} & \|
|000048c0| 6c 6f 77 65 72 2e 33 65 | 78 5c 68 62 6f 78 7b 32 |lower.3e|x\hbox{2|
|000048d0| 36 7d 20 26 20 20 0a 5c | 6c 6f 77 65 72 2e 33 65 |6} & .\|lower.3e|
|000048e0| 78 5c 68 62 6f 78 7b 33 | 32 7d 20 26 20 5c 6c 6f |x\hbox{3|2} & \lo|
|000048f0| 77 65 72 2e 33 65 78 5c | 68 62 6f 78 7b 34 37 2e |wer.3ex\|hbox{47.|
|00004900| 31 32 7d 20 26 20 5c 6c | 6f 77 65 72 2e 33 65 78 |12} & \l|ower.3ex|
|00004910| 5c 68 62 6f 78 7b 2e 31 | 30 34 38 7d 20 5c 5c 20 |\hbox{.1|048} \\ |
|00004920| 0a 46 44 32 20 26 20 24 | 31 30 5e 7b 2d 31 30 7d |.FD2 & $|10^{-10}|
|00004930| 24 20 26 20 32 36 20 26 | 20 35 38 20 26 20 35 33 |$ & 26 &| 58 & 53|
|00004940| 2e 37 39 20 26 20 2e 31 | 38 32 39 20 5c 5c 20 5c |.79 & .1|829 \\ \|
|00004950| 68 6c 69 6e 65 20 20 0a | 5c 6c 6f 77 65 72 2e 33 |hline .|\lower.3|
|00004960| 65 78 5c 68 62 6f 78 7b | 45 48 41 34 7d 20 26 20 |ex\hbox{|EHA4} & |
|00004970| 5c 6c 6f 77 65 72 2e 33 | 65 78 5c 68 62 6f 78 7b |\lower.3|ex\hbox{|
|00004980| 24 31 30 5e 7b 2d 31 32 | 7d 24 7d 20 26 20 5c 6c |$10^{-12|}$} & \l|
|00004990| 6f 77 65 72 2e 33 65 78 | 5c 68 62 6f 78 7b 33 30 |ower.3ex|\hbox{30|
|000049a0| 7d 20 26 20 20 0a 5c 6c | 6f 77 65 72 2e 33 65 78 |} & .\l|ower.3ex|
|000049b0| 5c 68 62 6f 78 7b 34 32 | 7d 20 26 20 5c 6c 6f 77 |\hbox{42|} & \low|
|000049c0| 65 72 2e 33 65 78 5c 68 | 62 6f 78 7b 35 36 2e 37 |er.3ex\h|box{56.7|
|000049d0| 36 7d 20 26 20 5c 6c 6f | 77 65 72 2e 33 65 78 5c |6} & \lo|wer.3ex\|
|000049e0| 68 62 6f 78 7b 2e 31 38 | 35 35 7d 20 5c 5c 20 20 |hbox{.18|55} \\ |
|000049f0| 0a 46 44 34 20 26 20 24 | 31 30 5e 7b 2d 31 32 7d |.FD4 & $|10^{-12}|
|00004a00| 24 20 26 20 33 30 20 26 | 20 31 33 32 20 26 20 38 |$ & 30 &| 132 & 8|
|00004a10| 31 2e 33 35 20 26 20 2e | 33 37 33 30 20 5c 5c 20 |1.35 & .|3730 \\ |
|00004a20| 5c 68 6c 69 6e 65 20 20 | 0a 5c 6c 6f 77 65 72 2e |\hline |.\lower.|
|00004a30| 33 65 78 5c 68 62 6f 78 | 7b 45 48 41 36 7d 20 26 |3ex\hbox|{EHA6} &|
|00004a40| 20 5c 6c 6f 77 65 72 2e | 33 65 78 5c 68 62 6f 78 | \lower.|3ex\hbox|
|00004a50| 7b 24 31 30 5e 7b 2d 31 | 32 7d 24 7d 20 26 20 5c |{$10^{-1|2}$} & \|
|00004a60| 6c 6f 77 65 72 2e 33 65 | 78 5c 68 62 6f 78 7b 33 |lower.3e|x\hbox{3|
|00004a70| 30 7d 20 26 20 20 0a 5c | 6c 6f 77 65 72 2e 33 65 |0} & .\|lower.3e|
|00004a80| 78 5c 68 62 6f 78 7b 34 | 38 7d 20 26 20 5c 6c 6f |x\hbox{4|8} & \lo|
|00004a90| 77 65 72 2e 33 65 78 5c | 68 62 6f 78 7b 35 38 2e |wer.3ex\|hbox{58.|
|00004aa0| 35 36 7d 20 26 20 5c 6c | 6f 77 65 72 2e 33 65 78 |56} & \l|ower.3ex|
|00004ab0| 5c 68 62 6f 78 7b 2e 31 | 39 35 32 7d 20 5c 5c 20 |\hbox{.1|952} \\ |
|00004ac0| 0a 46 44 36 20 26 20 24 | 31 30 5e 7b 2d 31 32 7d |.FD6 & $|10^{-12}|
|00004ad0| 24 20 26 20 33 30 20 26 | 20 31 39 38 20 26 20 31 |$ & 30 &| 198 & 1|
|00004ae0| 30 30 2e 36 20 26 20 2e | 33 32 37 38 20 5c 5c 20 |00.6 & .|3278 \\ |
|00004af0| 5c 68 6c 69 6e 65 20 20 | 0a 5c 65 6e 64 7b 74 61 |\hline |.\end{ta|
|00004b00| 62 75 6c 61 72 7d 7d 0a | 5c 65 6e 64 7b 63 65 6e |bular}}.|\end{cen|
|00004b10| 74 65 72 7d 20 0a 5c 6c | 61 62 65 6c 7b 64 69 66 |ter} .\l|abel{dif|
|00004b20| 66 73 74 61 74 73 7d 20 | 0a 5c 65 6e 64 7b 74 61 |fstats} |.\end{ta|
|00004b30| 62 6c 65 7d 20 20 0a 0a | 49 6e 20 6f 75 72 20 66 |ble} ..|In our f|
|00004b40| 69 72 73 74 20 73 65 74 | 20 6f 66 20 65 78 70 65 |irst set| of expe|
|00004b50| 72 69 6d 65 6e 74 73 2c | 20 77 65 20 74 6f 6f 6b |riments,| we took|
|00004b60| 20 24 63 3d 64 3d 31 30 | 24 20 61 6e 64 20 75 73 | $c=d=10|$ and us|
|00004b70| 65 64 20 72 69 67 68 74 | 20 0a 70 72 65 63 6f 6e |ed right| .precon|
|00004b80| 64 69 74 69 6f 6e 69 6e | 67 20 77 69 74 68 20 61 |ditionin|g with a|
|00004b90| 20 66 61 73 74 20 50 6f | 69 73 73 6f 6e 20 73 6f | fast Po|isson so|
|00004ba0| 6c 76 65 72 20 66 72 6f | 6d 20 7b 5c 66 69 73 68 |lver fro|m {\fish|
|00004bb0| 70 61 63 6b 7d 20 20 0a | 5c 63 69 74 65 7b 53 77 |pack} .|\cite{Sw|
|00004bc0| 61 72 7a 74 72 61 75 62 | 65 72 2d 53 77 65 65 74 |arztraub|er-Sweet|
|00004bd0| 7d 2c 20 77 68 69 63 68 | 20 69 73 20 76 65 72 79 |}, which| is very|
|00004be0| 20 65 66 66 65 63 74 69 | 76 65 20 66 6f 72 20 74 | effecti|ve for t|
|00004bf0| 68 65 73 65 20 20 0a 66 | 61 69 72 6c 79 20 73 6d |hese .f|airly sm|
|00004c00| 61 6c 6c 20 76 61 6c 75 | 65 73 20 6f 66 20 24 63 |all valu|es of $c|
|00004c10| 24 20 61 6e 64 20 24 64 | 24 2e 20 57 65 20 66 69 |$ and $d|$. We fi|
|00004c20| 72 73 74 20 73 74 61 72 | 74 65 64 20 65 61 63 68 |rst star|ted each|
|00004c30| 20 6d 65 74 68 6f 64 20 | 0a 77 69 74 68 20 7a 65 | method |.with ze|
|00004c40| 72 6f 20 61 73 20 74 68 | 65 20 69 6e 69 74 69 61 |ro as th|e initia|
|00004c50| 6c 20 61 70 70 72 6f 78 | 69 6d 61 74 65 20 73 6f |l approx|imate so|
|00004c60| 6c 75 74 69 6f 6e 20 61 | 6e 64 20 61 6c 6c 6f 77 |lution a|nd allow|
|00004c70| 65 64 20 69 74 20 20 0a | 74 6f 20 72 75 6e 20 66 |ed it .|to run f|
|00004c80| 6f 72 20 34 30 20 7b 5c | 67 6d 72 65 73 6d 7d 20 |or 40 {\|gmresm} |
|00004c90| 69 74 65 72 61 74 69 6f | 6e 73 2c 20 61 66 74 65 |iteratio|ns, afte|
|00004ca0| 72 20 77 68 69 63 68 20 | 74 68 65 20 6c 69 6d 69 |r which |the limi|
|00004cb0| 74 20 6f 66 20 72 65 73 | 69 64 75 61 6c 20 20 0a |t of res|idual .|
|00004cc0| 6e 6f 72 6d 20 72 65 64 | 75 63 74 69 6f 6e 20 68 |norm red|uction h|
|00004cd0| 61 64 20 62 65 65 6e 20 | 72 65 61 63 68 65 64 2e |ad been |reached.|
|00004ce0| 20 46 69 67 75 72 65 20 | 5c 72 65 66 7b 70 64 65 | Figure |\ref{pde|
|00004cf0| 70 7d 20 73 68 6f 77 73 | 20 70 6c 6f 74 73 20 20 |p} shows| plots |
|00004d00| 0a 6f 66 20 74 68 65 20 | 6c 6f 67 61 72 69 74 68 |.of the |logarith|
|00004d10| 6d 20 6f 66 20 74 68 65 | 20 45 75 63 6c 69 64 65 |m of the| Euclide|
|00004d20| 61 6e 20 6e 6f 72 6d 20 | 6f 66 20 74 68 65 20 72 |an norm |of the r|
|00004d30| 65 73 69 64 75 61 6c 20 | 76 65 72 73 75 73 20 20 |esidual |versus |
|00004d40| 0a 74 68 65 20 6e 75 6d | 62 65 72 20 6f 66 20 7b |.the num|ber of {|
|00004d50| 5c 67 6d 72 65 73 6d 7d | 20 69 74 65 72 61 74 69 |\gmresm}| iterati|
|00004d60| 6f 6e 73 20 66 6f 72 20 | 74 68 65 20 74 68 72 65 |ons for |the thre|
|00004d70| 65 20 6d 65 74 68 6f 64 | 73 2e 20 57 65 20 6e 6f |e method|s. We no|
|00004d80| 74 65 0a 74 68 61 74 20 | 69 6e 20 20 46 69 67 2e |te.that |in Fig.|
|00004d90| 7e 5c 72 65 66 7b 70 64 | 65 70 7d 20 61 6e 64 20 |~\ref{pd|ep} and |
|00004da0| 69 6e 20 61 6c 6c 20 6f | 74 68 65 72 20 66 69 67 |in all o|ther fig|
|00004db0| 75 72 65 73 20 62 65 6c | 6f 77 2c 20 74 68 65 20 |ures bel|ow, the |
|00004dc0| 70 6c 6f 74 74 65 64 0a | 72 65 73 69 64 75 61 6c |plotted.|residual|
|00004dd0| 20 6e 6f 72 6d 73 20 77 | 65 72 65 20 6e 6f 74 20 | norms w|ere not |
|00004de0| 74 68 65 20 76 61 6c 75 | 65 73 20 6d 61 69 6e 74 |the valu|es maint|
|00004df0| 61 69 6e 65 64 20 62 79 | 20 7b 5c 67 6d 72 65 73 |ained by| {\gmres|
|00004e00| 6d 7d 2c 20 62 75 74 20 | 72 61 74 68 65 72 0a 77 |m}, but |rather.w|
|00004e10| 65 72 65 20 63 6f 6d 70 | 75 74 65 64 20 61 73 20 |ere comp|uted as |
|00004e20| 61 63 63 75 72 61 74 65 | 6c 79 20 61 73 20 70 6f |accurate|ly as po|
|00004e30| 73 73 69 62 6c 65 20 60 | 60 66 72 6f 6d 20 73 63 |ssible `|`from sc|
|00004e40| 72 61 74 63 68 2e 27 27 | 20 20 54 68 61 74 20 69 |ratch.''| That i|
|00004e50| 73 2c 20 0a 61 74 20 65 | 61 63 68 20 7b 5c 67 6d |s, .at e|ach {\gm|
|00004e60| 72 65 73 6d 7d 20 69 74 | 65 72 61 74 69 6f 6e 2c |resm} it|eration,|
|00004e70| 20 74 68 65 20 63 75 72 | 72 65 6e 74 20 61 70 70 | the cur|rent app|
|00004e80| 72 6f 78 69 6d 61 74 65 | 20 73 6f 6c 75 74 69 6f |roximate| solutio|
|00004e90| 6e 20 77 61 73 20 0a 66 | 6f 72 6d 65 64 20 61 6e |n was .f|ormed an|
|00004ea0| 64 20 69 74 73 20 70 72 | 6f 64 75 63 74 20 77 69 |d its pr|oduct wi|
|00004eb0| 74 68 20 74 68 65 20 63 | 6f 65 66 66 69 63 69 65 |th the c|oefficie|
|00004ec0| 6e 74 20 6d 61 74 72 69 | 78 20 77 61 73 20 73 75 |nt matri|x was su|
|00004ed0| 62 74 72 61 63 74 65 64 | 20 0a 66 72 6f 6d 20 74 |btracted| .from t|
|00004ee0| 68 65 20 72 69 67 68 74 | 2d 68 61 6e 64 20 73 69 |he right|-hand si|
|00004ef0| 64 65 2c 20 61 6c 6c 20 | 69 6e 20 64 6f 75 62 6c |de, all |in doubl|
|00004f00| 65 20 70 72 65 63 69 73 | 69 6f 6e 2e 20 20 0a 49 |e precis|ion. .I|
|00004f10| 74 20 77 61 73 20 69 6d | 70 6f 72 74 61 6e 74 20 |t was im|portant |
|00004f20| 74 6f 20 63 6f 6d 70 75 | 74 65 20 74 68 65 20 72 |to compu|te the r|
|00004f30| 65 73 69 64 75 61 6c 20 | 6e 6f 72 6d 73 20 69 6e |esidual |norms in|
|00004f40| 20 74 68 69 73 20 77 61 | 79 20 62 65 63 61 75 73 | this wa|y becaus|
|00004f50| 65 20 20 0a 74 68 65 20 | 76 61 6c 75 65 73 20 6d |e .the |values m|
|00004f60| 61 69 6e 74 61 69 6e 65 | 64 20 62 79 20 7b 5c 67 |aintaine|d by {\g|
|00004f70| 6d 72 65 73 6d 7d 20 62 | 65 63 6f 6d 65 20 69 6e |mresm} b|ecome in|
|00004f80| 63 72 65 61 73 69 6e 67 | 6c 79 20 75 6e 74 72 75 |creasing|ly untru|
|00004f90| 73 74 77 6f 72 74 68 79 | 20 20 0a 61 73 20 74 68 |stworthy| .as th|
|00004fa0| 65 20 6c 69 6d 69 74 73 | 20 6f 66 20 72 65 73 69 |e limits| of resi|
|00004fb0| 64 75 61 6c 20 6e 6f 72 | 6d 20 72 65 64 75 63 74 |dual nor|m reduct|
|00004fc0| 69 6f 6e 20 61 72 65 20 | 6e 65 61 72 65 64 3b 20 |ion are |neared; |
|00004fd0| 73 65 65 20 5c 63 69 74 | 65 7b 57 61 6c 6b 65 72 |see \cit|e{Walker|
|00004fe0| 38 38 7d 2e 20 20 0a 49 | 74 20 69 73 20 73 65 65 |88}. .I|t is see|
|00004ff0| 6e 20 69 6e 20 46 69 67 | 2e 7e 5c 72 65 66 7b 70 |n in Fig|.~\ref{p|
|00005000| 64 65 70 7d 20 74 68 61 | 74 20 41 6c 67 6f 72 69 |dep} tha|t Algori|
|00005010| 74 68 6d 20 45 48 41 20 | 61 63 68 69 65 76 65 64 |thm EHA |achieved|
|00005020| 20 20 0a 74 68 65 20 73 | 61 6d 65 20 75 6c 74 69 | .the s|ame ulti|
|00005030| 6d 61 74 65 20 6c 65 76 | 65 6c 20 6f 66 20 72 65 |mate lev|el of re|
|00005040| 73 69 64 75 61 6c 20 6e | 6f 72 6d 20 72 65 64 75 |sidual n|orm redu|
|00005050| 63 74 69 6f 6e 20 61 73 | 20 74 68 65 20 46 44 50 |ction as| the FDP|
|00005060| 20 20 0a 6d 65 74 68 6f | 64 20 61 6e 64 20 72 65 | .metho|d and re|
|00005070| 71 75 69 72 65 64 20 6f | 6e 6c 79 20 61 20 66 65 |quired o|nly a fe|
|00005080| 77 20 6d 6f 72 65 20 7b | 5c 67 6d 72 65 73 6d 7d |w more {|\gmresm}|
|00005090| 20 69 74 65 72 61 74 69 | 6f 6e 73 20 74 6f 20 64 | iterati|ons to d|
|000050a0| 6f 0a 73 6f 2e 20 20 0a | 20 0a 5c 62 65 67 69 6e |o.so. .| .\begin|
|000050b0| 7b 66 69 67 75 72 65 7d | 5b 74 5d 20 0a 5c 76 73 |{figure}|[t] .\vs|
|000050c0| 70 61 63 65 7b 33 69 6e | 7d 0a 5c 63 61 70 74 69 |pace{3in|}.\capti|
|000050d0| 6f 6e 7b 7b 5c 72 6d 20 | 4c 6f 67 7d 24 5f 7b 31 |on{{\rm |Log}$_{1|
|000050e0| 30 7d 24 20 6f 66 20 74 | 68 65 20 72 65 73 69 64 |0}$ of t|he resid|
|000050f0| 75 61 6c 20 6e 6f 72 6d | 20 76 65 72 73 75 73 20 |ual norm| versus |
|00005100| 74 68 65 20 6e 75 6d 62 | 65 72 20 6f 66 20 0a 7b |the numb|er of .{|
|00005110| 5c 72 6d 20 47 4d 52 45 | 53 7d 24 28 6d 29 24 20 |\rm GMRE|S}$(m)$ |
|00005120| 69 74 65 72 61 74 69 6f | 6e 73 20 66 6f 72 20 24 |iteratio|ns for $|
|00005130| 63 3d 64 3d 31 30 24 20 | 77 69 74 68 20 66 61 73 |c=d=10$ |with fas|
|00005140| 74 20 50 6f 69 73 73 6f | 6e 0a 70 72 65 63 6f 6e |t Poisso|n.precon|
|00005150| 64 69 74 69 6f 6e 69 6e | 67 2e 20 53 6f 6c 69 64 |ditionin|g. Solid|
|00005160| 20 63 75 72 76 65 3a 20 | 41 6c 67 6f 72 69 74 68 | curve: |Algorith|
|00005170| 6d 20 7b 5c 72 6d 20 45 | 48 41 7d 3b 20 64 6f 74 |m {\rm E|HA}; dot|
|00005180| 74 65 64 20 0a 63 75 72 | 76 65 3a 20 7b 5c 72 6d |ted .cur|ve: {\rm|
|00005190| 20 46 44 50 7d 20 6d 65 | 74 68 6f 64 3b 20 64 61 | FDP} me|thod; da|
|000051a0| 73 68 65 64 20 63 75 72 | 76 65 3a 20 7b 5c 72 6d |shed cur|ve: {\rm|
|000051b0| 20 46 53 50 7d 20 6d 65 | 74 68 6f 64 2e 7d 20 0a | FSP} me|thod.} .|
|000051c0| 5c 6c 61 62 65 6c 7b 70 | 64 65 70 7d 20 0a 5c 65 |\label{p|dep} .\e|
|000051d0| 6e 64 7b 66 69 67 75 72 | 65 7d 20 0a 20 0a 20 0a |nd{figur|e} . . .|
|000051e0| 49 6e 20 6f 75 72 20 73 | 65 63 6f 6e 64 20 73 65 |In our s|econd se|
|000051f0| 74 20 6f 66 20 65 78 70 | 65 72 69 6d 65 6e 74 73 |t of exp|eriments|
|00005200| 2c 20 77 65 20 74 6f 6f | 6b 20 24 63 3d 64 3d 31 |, we too|k $c=d=1|
|00005210| 30 30 24 20 61 6e 64 20 | 63 61 72 72 69 65 64 20 |00$ and |carried |
|00005220| 6f 75 74 20 20 0a 74 72 | 69 61 6c 73 20 61 6e 61 |out .tr|ials ana|
|00005230| 6c 6f 67 6f 75 73 20 74 | 6f 20 74 68 6f 73 65 20 |logous t|o those |
|00005240| 69 6e 20 74 68 65 20 66 | 69 72 73 74 20 73 65 74 |in the f|irst set|
|00005250| 20 61 62 6f 76 65 2e 20 | 4e 6f 20 70 72 65 63 6f | above. |No preco|
|00005260| 6e 64 69 74 69 6f 6e 69 | 6e 67 20 20 0a 77 61 73 |nditioni|ng .was|
|00005270| 20 75 73 65 64 20 69 6e | 20 74 68 65 73 65 20 65 | used in| these e|
|00005280| 78 70 65 72 69 6d 65 6e | 74 73 2c 20 62 6f 74 68 |xperimen|ts, both|
|00005290| 20 62 65 63 61 75 73 65 | 20 77 65 20 77 61 6e 74 | because| we want|
|000052a0| 65 64 20 74 6f 20 63 6f | 6d 70 61 72 65 20 0a 74 |ed to co|mpare .t|
|000052b0| 68 65 20 6d 65 74 68 6f | 64 73 20 77 69 74 68 6f |he metho|ds witho|
|000052c0| 75 74 20 70 72 65 63 6f | 6e 64 69 74 69 6f 6e 69 |ut preco|nditioni|
|000052d0| 6e 67 20 61 6e 64 20 62 | 65 63 61 75 73 65 20 74 |ng and b|ecause t|
|000052e0| 68 65 20 66 61 73 74 20 | 20 0a 50 6f 69 73 73 6f |he fast | .Poisso|
|000052f0| 6e 20 70 72 65 63 6f 6e | 64 69 74 69 6f 6e 69 6e |n precon|ditionin|
|00005300| 67 20 75 73 65 64 20 69 | 6e 20 74 68 65 20 66 69 |g used i|n the fi|
|00005310| 72 73 74 20 73 65 74 20 | 6f 66 20 65 78 70 65 72 |rst set |of exper|
|00005320| 69 6d 65 6e 74 73 20 69 | 73 20 0a 6e 6f 74 20 63 |iments i|s .not c|
|00005330| 6f 73 74 20 65 66 66 65 | 63 74 69 76 65 20 66 6f |ost effe|ctive fo|
|00005340| 72 20 74 68 65 73 65 20 | 6c 61 72 67 65 20 76 61 |r these |large va|
|00005350| 6c 75 65 73 20 6f 66 20 | 24 63 24 20 61 6e 64 20 |lues of |$c$ and |
|00005360| 24 64 24 2e 20 57 65 20 | 66 69 72 73 74 20 20 0a |$d$. We |first .|
|00005370| 61 6c 6c 6f 77 65 64 20 | 65 61 63 68 20 6d 65 74 |allowed |each met|
|00005380| 68 6f 64 20 74 6f 20 72 | 75 6e 20 66 6f 72 20 36 |hod to r|un for 6|
|00005390| 30 30 20 7b 5c 67 6d 72 | 65 73 6d 7d 20 69 74 65 |00 {\gmr|esm} ite|
|000053a0| 72 61 74 69 6f 6e 73 2c | 20 0a 73 74 61 72 74 69 |rations,| .starti|
|000053b0| 6e 67 20 77 69 74 68 20 | 7a 65 72 6f 20 61 73 20 |ng with |zero as |
|000053c0| 74 68 65 20 69 6e 69 74 | 69 61 6c 20 61 70 70 72 |the init|ial appr|
|000053d0| 6f 78 69 6d 61 74 65 20 | 73 6f 6c 75 74 69 6f 6e |oximate |solution|
|000053e0| 2c 20 61 66 74 65 72 20 | 77 68 69 63 68 20 20 0a |, after |which .|
|000053f0| 74 68 65 20 6c 69 6d 69 | 74 20 6f 66 20 72 65 73 |the limi|t of res|
|00005400| 69 64 75 61 6c 20 6e 6f | 72 6d 20 72 65 64 75 63 |idual no|rm reduc|
|00005410| 74 69 6f 6e 20 68 61 64 | 20 62 65 65 6e 20 72 65 |tion had| been re|
|00005420| 61 63 68 65 64 2e 20 0a | 20 0a 0a 5c 73 65 63 74 |ached. .| ..\sect|
|00005430| 69 6f 6e 2a 7b 41 63 6b | 6e 6f 77 6c 65 64 67 6d |ion*{Ack|nowledgm|
|00005440| 65 6e 74 73 7d 0a 54 68 | 65 20 61 75 74 68 6f 72 |ents}.Th|e author|
|00005450| 20 74 68 61 6e 6b 73 20 | 74 68 65 20 61 6e 6f 6e | thanks |the anon|
|00005460| 79 6d 6f 75 73 20 61 75 | 74 68 6f 72 73 20 77 68 |ymous au|thors wh|
|00005470| 6f 73 65 20 77 6f 72 6b | 20 6c 61 72 67 65 6c 79 |ose work| largely|
|00005480| 0a 63 6f 6e 73 74 69 74 | 75 74 65 73 20 74 68 69 |.constit|utes thi|
|00005490| 73 20 73 61 6d 70 6c 65 | 20 66 69 6c 65 2e 20 48 |s sample| file. H|
|000054a0| 65 20 61 6c 73 6f 20 74 | 68 61 6e 6b 73 20 74 68 |e also t|hanks th|
|000054b0| 65 20 49 4e 46 4f 2d 54 | 65 58 20 6d 61 69 6c 69 |e INFO-T|eX maili|
|000054c0| 6e 67 0a 6c 69 73 74 20 | 66 6f 72 20 74 68 65 20 |ng.list |for the |
|000054d0| 76 61 6c 75 61 62 6c 65 | 20 69 6e 64 69 72 65 63 |valuable| indirec|
|000054e0| 74 20 61 73 73 69 73 74 | 61 6e 63 65 20 68 65 20 |t assist|ance he |
|000054f0| 72 65 63 65 69 76 65 64 | 2e 0a 20 0a 20 0a 5c 62 |received|.. . .\b|
|00005500| 65 67 69 6e 7b 74 68 65 | 62 69 62 6c 69 6f 67 72 |egin{the|bibliogr|
|00005510| 61 70 68 79 7d 7b 31 30 | 7d 20 0a 5c 62 69 62 69 |aphy}{10|} .\bibi|
|00005520| 74 65 6d 7b 62 73 7d 20 | 7b 5c 73 63 20 52 2e 7e |tem{bs} |{\sc R.~|
|00005530| 41 2e 20 42 72 75 61 6c | 64 69 20 61 6e 64 20 42 |A. Brual|di and B|
|00005540| 2e 7e 4c 2e 20 53 68 61 | 64 65 72 7d 2c 20 0a 7b |.~L. Sha|der}, .{|
|00005550| 5c 65 6d 20 4f 6e 20 73 | 69 67 6e 2d 6e 6f 6e 73 |\em On s|ign-nons|
|00005560| 69 6e 67 75 6c 61 72 20 | 6d 61 74 72 69 63 65 73 |ingular |matrices|
|00005570| 20 61 6e 64 20 74 68 65 | 20 63 6f 6e 76 65 72 73 | and the| convers|
|00005580| 69 6f 6e 20 6f 66 20 74 | 68 65 0a 70 65 72 6d 61 |ion of t|he.perma|
|00005590| 6e 65 6e 74 20 69 6e 74 | 6f 20 74 68 65 20 64 65 |nent int|o the de|
|000055a0| 74 65 72 6d 69 6e 61 6e | 74 7d 2c 20 69 6e 20 41 |terminan|t}, in A|
|000055b0| 70 70 6c 69 65 64 20 47 | 65 6f 6d 65 74 72 79 20 |pplied G|eometry |
|000055c0| 61 6e 64 0a 44 69 73 63 | 72 65 74 65 20 4d 61 74 |and.Disc|rete Mat|
|000055d0| 68 65 6d 61 74 69 63 73 | 2c 20 54 68 65 20 56 69 |hematics|, The Vi|
|000055e0| 63 74 6f 72 20 4b 6c 65 | 65 20 46 65 73 74 73 63 |ctor Kle|e Festsc|
|000055f0| 68 72 69 66 74 2c 20 50 | 2e 0a 47 72 69 74 7a 6d |hrift, P|..Gritzm|
|00005600| 61 6e 6e 20 61 6e 64 20 | 42 2e 20 53 74 75 72 6d |ann and |B. Sturm|
|00005610| 66 65 6c 73 2c 20 65 64 | 73 2e 2c 20 41 6d 65 72 |fels, ed|s., Amer|
|00005620| 69 63 61 6e 20 4d 61 74 | 68 65 6d 61 74 69 63 61 |ican Mat|hematica|
|00005630| 6c 0a 53 6f 63 69 65 74 | 79 2c 20 50 72 6f 76 69 |l.Societ|y, Provi|
|00005640| 64 65 6e 63 65 2c 20 52 | 49 2c 20 31 39 39 31 2c |dence, R|I, 1991,|
|00005650| 20 70 70 2e 20 31 31 37 | 2d 2d 31 33 34 2e 0a 20 | pp. 117|--134.. |
|00005660| 0a 5c 62 69 62 69 74 65 | 6d 7b 64 6a 64 7d 20 7b |.\bibite|m{djd} {|
|00005670| 5c 73 63 20 4a 2e 20 44 | 72 65 77 2c 20 43 2e 7e |\sc J. D|rew, C.~|
|00005680| 52 2e 20 4a 6f 68 6e 73 | 6f 6e 2c 20 61 6e 64 20 |R. Johns|on, and |
|00005690| 50 2e 20 76 61 6e 20 64 | 65 6e 20 44 72 69 65 73 |P. van d|en Dries|
|000056a0| 73 63 68 65 7d 2c 20 0a | 7b 5c 65 6d 20 53 74 72 |sche}, .|{\em Str|
|000056b0| 6f 6e 67 20 66 6f 72 6d | 73 20 6f 66 20 6e 6f 6e |ong form|s of non|
|000056c0| 73 69 6e 67 75 6c 61 72 | 69 74 79 7d 2c 20 4c 69 |singular|ity}, Li|
|000056d0| 6e 65 61 72 20 41 6c 67 | 65 62 72 61 20 41 70 70 |near Alg|ebra App|
|000056e0| 6c 2e 2c 0a 31 36 32 20 | 28 31 39 39 32 29 2c 20 |l.,.162 |(1992), |
|000056f0| 74 6f 20 61 70 70 65 61 | 72 2e 20 0a 20 0a 5c 62 |to appea|r. . .\b|
|00005700| 69 62 69 74 65 6d 7b 67 | 7d 20 7b 5c 73 63 20 50 |ibitem{g|} {\sc P|
|00005710| 2e 7e 4d 2e 20 47 69 62 | 73 6f 6e 7d 2c 20 7b 5c |.~M. Gib|son}, {\|
|00005720| 65 6d 20 43 6f 6e 76 65 | 72 73 69 6f 6e 20 6f 66 |em Conve|rsion of|
|00005730| 20 74 68 65 20 70 65 72 | 6d 61 6e 65 6e 74 20 69 | the per|manent i|
|00005740| 6e 74 6f 20 74 68 65 20 | 0a 64 65 74 65 72 6d 69 |nto the |.determi|
|00005750| 6e 61 6e 74 7d 2c 20 50 | 72 6f 63 2e 20 41 6d 65 |nant}, P|roc. Ame|
|00005760| 72 2e 20 4d 61 74 68 2e | 20 53 6f 63 2e 2c 20 32 |r. Math.| Soc., 2|
|00005770| 37 20 28 31 39 37 31 29 | 2c 0a 70 70 2e 7e 34 37 |7 (1971)|,.pp.~47|
|00005780| 31 2d 2d 34 37 36 2e 0a | 20 0a 5c 62 69 62 69 74 |1--476..| .\bibit|
|00005790| 65 6d 7b 6b 6c 6d 7d 20 | 0a 7b 5c 73 63 20 56 2e |em{klm} |.{\sc V.|
|000057a0| 7e 4b 6c 65 65 2c 20 52 | 2e 7e 4c 61 64 6e 65 72 |~Klee, R|.~Ladner|
|000057b0| 2c 20 61 6e 64 20 52 2e | 7e 4d 61 6e 62 65 72 7d |, and R.|~Manber}|
|000057c0| 2c 20 7b 5c 69 74 0a 53 | 69 67 6e 73 6f 6c 76 61 |, {\it.S|ignsolva|
|000057d0| 62 69 6c 69 74 79 20 72 | 65 76 69 73 69 74 65 64 |bility r|evisited|
|000057e0| 7d 2c 20 4c 69 6e 65 61 | 72 20 41 6c 67 65 62 72 |}, Linea|r Algebr|
|000057f0| 61 20 41 70 70 6c 2e 2c | 20 35 39 0a 28 31 39 38 |a Appl.,| 59.(198|
|00005800| 34 29 2c 20 70 70 2e 7e | 31 33 31 2d 2d 31 35 37 |4), pp.~|131--157|
|00005810| 2e 0a 20 0a 5c 62 69 62 | 69 74 65 6d 7b 6d 7d 20 |.. .\bib|item{m} |
|00005820| 0a 7b 5c 73 63 20 4b 2e | 20 4d 75 72 6f 74 61 7d |.{\sc K.| Murota}|
|00005830| 2c 20 4c 55 2d 7b 5c 65 | 6d 20 64 65 63 6f 6d 70 |, LU-{\e|m decomp|
|00005840| 6f 73 69 74 69 6f 6e 20 | 6f 66 20 61 20 6d 61 74 |osition |of a mat|
|00005850| 72 69 78 20 77 69 74 68 | 0a 65 6e 74 72 69 65 73 |rix with|.entries|
|00005860| 20 6f 66 20 64 69 66 66 | 65 72 65 6e 74 20 6b 69 | of diff|erent ki|
|00005870| 6e 64 73 7d 2c 20 4c 69 | 6e 65 61 72 20 41 6c 67 |nds}, Li|near Alg|
|00005880| 65 62 72 61 20 41 70 70 | 6c 2e 2c 20 34 39 0a 28 |ebra App|l., 49.(|
|00005890| 31 39 38 33 29 2c 20 70 | 70 2e 7e 32 37 35 2d 2d |1983), p|p.~275--|
|000058a0| 32 38 33 2e 0a 0a 5c 62 | 69 62 69 74 65 6d 7b 41 |283...\b|ibitem{A|
|000058b0| 78 65 6c 73 73 6f 6e 7d | 0a 7b 5c 73 63 20 4f 2e |xelsson}|.{\sc O.|
|000058c0| 7e 41 78 65 6c 73 73 6f | 6e 7d 2c 20 7b 5c 65 6d |~Axelsso|n}, {\em|
|000058d0| 20 43 6f 6e 6a 75 67 61 | 74 65 20 67 72 61 64 69 | Conjuga|te gradi|
|000058e0| 65 6e 74 20 74 79 70 65 | 20 6d 65 74 68 6f 64 73 |ent type| methods|
|000058f0| 20 66 6f 72 20 75 6e 73 | 79 6d 6d 65 74 72 69 63 | for uns|ymmetric|
|00005900| 20 61 6e 64 0a 20 20 69 | 6e 63 6f 6e 73 69 73 74 | and. i|nconsist|
|00005910| 65 6e 74 20 73 79 73 74 | 65 6d 73 20 6f 66 20 6c |ent syst|ems of l|
|00005920| 69 6e 65 61 72 20 65 71 | 75 61 74 69 6f 6e 73 7d |inear eq|uations}|
|00005930| 2c 20 4c 69 6e 65 61 72 | 20 41 6c 67 65 62 72 61 |, Linear| Algebra|
|00005940| 20 41 70 70 6c 2e 2c 20 | 32 39 20 28 31 39 38 30 | Appl., |29 (1980|
|00005950| 29 2c 0a 20 20 70 70 2e | 7e 31 2d 2d 31 36 2e 0a |),. pp.|~1--16..|
|00005960| 0a 5c 62 69 62 69 74 65 | 6d 7b 42 72 6f 77 6e 2d |.\bibite|m{Brown-|
|00005970| 53 61 61 64 31 7d 0a 7b | 5c 73 63 20 50 2e 7e 4e |Saad1}.{|\sc P.~N|
|00005980| 2e 20 42 72 6f 77 6e 20 | 61 6e 64 20 59 2e 7e 53 |. Brown |and Y.~S|
|00005990| 61 61 64 7d 2c 20 7b 5c | 65 6d 20 48 79 62 72 69 |aad}, {\|em Hybri|
|000059a0| 64 20 7b 4b 7d 72 79 6c | 6f 76 20 6d 65 74 68 6f |d {K}ryl|ov metho|
|000059b0| 64 73 20 66 6f 72 20 6e | 6f 6e 6c 69 6e 65 61 72 |ds for n|onlinear|
|000059c0| 0a 20 20 73 79 73 74 65 | 6d 73 20 6f 66 20 65 71 |. syste|ms of eq|
|000059d0| 75 61 74 69 6f 6e 73 7d | 2c 20 53 49 41 4d 20 4a |uations}|, SIAM J|
|000059e0| 2e 20 53 63 69 2e 20 53 | 74 61 74 69 73 74 2e 20 |. Sci. S|tatist. |
|000059f0| 43 6f 6d 70 75 74 2e 2c | 20 31 31 20 28 31 39 39 |Comput.,| 11 (199|
|00005a00| 30 29 2c 20 0a 20 20 70 | 70 2e 7e 34 35 30 2d 2d |0), . p|p.~450--|
|00005a10| 34 38 31 2e 0a 0a 5c 62 | 69 62 69 74 65 6d 7b 44 |481...\b|ibitem{D|
|00005a20| 45 53 7d 0a 7b 5c 73 63 | 20 52 2e 7e 53 2e 20 44 |ES}.{\sc| R.~S. D|
|00005a30| 65 6d 62 6f 2c 20 53 2e | 7e 43 2e 20 45 69 73 65 |embo, S.|~C. Eise|
|00005a40| 6e 73 74 61 74 2c 20 61 | 6e 64 20 54 2e 7e 53 74 |nstat, a|nd T.~St|
|00005a50| 65 69 68 61 75 67 7d 2c | 20 7b 5c 65 6d 20 49 6e |eihaug},| {\em In|
|00005a60| 65 78 61 63 74 20 7b 4e | 7d 65 77 74 6f 6e 0a 20 |exact {N|}ewton. |
|00005a70| 20 6d 65 74 68 6f 64 73 | 7d 2c 20 53 49 41 4d 20 | methods|}, SIAM |
|00005a80| 4a 2e 20 4e 75 6d 65 72 | 2e 20 41 6e 61 6c 2e 2c |J. Numer|. Anal.,|
|00005a90| 20 31 39 20 28 31 39 38 | 32 29 2c 20 70 70 2e 7e | 19 (198|2), pp.~|
|00005aa0| 34 30 30 2d 2d 34 30 38 | 2e 0a 0a 5c 62 69 62 69 |400--408|...\bibi|
|00005ab0| 74 65 6d 7b 45 45 53 7d | 0a 7b 5c 73 63 20 53 2e |tem{EES}|.{\sc S.|
|00005ac0| 7e 43 2e 20 45 69 73 65 | 6e 73 74 61 74 2c 20 48 |~C. Eise|nstat, H|
|00005ad0| 2e 7e 43 2e 20 45 6c 6d | 61 6e 2c 20 61 6e 64 20 |.~C. Elm|an, and |
|00005ae0| 4d 2e 7e 48 2e 20 53 63 | 68 75 6c 74 7a 7d 2c 20 |M.~H. Sc|hultz}, |
|00005af0| 7b 5c 65 6d 20 56 61 72 | 69 61 74 69 6f 6e 61 6c |{\em Var|iational|
|00005b00| 0a 20 20 69 74 65 72 61 | 74 69 76 65 20 6d 65 74 |. itera|tive met|
|00005b10| 68 6f 64 73 20 66 6f 72 | 20 6e 6f 6e 73 79 6d 6d |hods for| nonsymm|
|00005b20| 65 74 72 69 63 20 73 79 | 73 74 65 6d 73 20 6f 66 |etric sy|stems of|
|00005b30| 20 6c 69 6e 65 61 72 20 | 65 71 75 61 74 69 6f 6e | linear |equation|
|00005b40| 73 7d 2c 20 53 49 41 4d | 20 4a 2e 0a 20 20 4e 75 |s}, SIAM| J.. Nu|
|00005b50| 6d 65 72 2e 20 41 6e 61 | 6c 2e 2c 20 32 30 20 28 |mer. Ana|l., 20 (|
|00005b60| 31 39 38 33 29 2c 20 70 | 70 2e 7e 33 34 35 2d 2d |1983), p|p.~345--|
|00005b70| 33 35 37 2e 0a 0a 5c 62 | 69 62 69 74 65 6d 7b 45 |357...\b|ibitem{E|
|00005b80| 6c 6d 61 6e 7d 0a 7b 5c | 73 63 20 48 2e 7e 43 2e |lman}.{\|sc H.~C.|
|00005b90| 20 45 6c 6d 61 6e 7d 2c | 20 7b 5c 65 6d 20 49 74 | Elman},| {\em It|
|00005ba0| 65 72 61 74 69 76 65 20 | 6d 65 74 68 6f 64 73 20 |erative |methods |
|00005bb0| 66 6f 72 20 6c 61 72 67 | 65 2c 20 73 70 61 72 73 |for larg|e, spars|
|00005bc0| 65 2c 20 6e 6f 6e 73 79 | 6d 6d 65 74 72 69 63 0a |e, nonsy|mmetric.|
|00005bd0| 20 20 73 79 73 74 65 6d | 73 20 6f 66 20 6c 69 6e | system|s of lin|
|00005be0| 65 61 72 20 65 71 75 61 | 74 69 6f 6e 73 7d 2c 20 |ear equa|tions}, |
|00005bf0| 50 68 2e 44 2e 20 74 68 | 65 73 69 73 2c 20 44 65 |Ph.D. th|esis, De|
|00005c00| 70 61 72 74 6d 65 6e 74 | 20 6f 66 20 43 6f 6d 70 |partment| of Comp|
|00005c10| 75 74 65 72 0a 20 20 53 | 63 69 65 6e 63 65 2c 20 |uter. S|cience, |
|00005c20| 59 61 6c 65 20 55 6e 69 | 76 65 72 73 69 74 79 2c |Yale Uni|versity,|
|00005c30| 20 4e 65 77 20 48 61 76 | 65 6e 2c 20 43 54 2c 20 | New Hav|en, CT, |
|00005c40| 31 39 38 32 2e 0a 0a 5c | 62 69 62 69 74 65 6d 7b |1982...\|bibitem{|
|00005c50| 47 6c 6f 4b 52 38 35 7d | 0a 7b 5c 73 63 20 52 2e |GloKR85}|.{\sc R.|
|00005c60| 7e 47 6c 6f 77 69 6e 73 | 6b 69 2c 20 48 2e 7e 42 |~Glowins|ki, H.~B|
|00005c70| 2e 20 4b 65 6c 6c 65 72 | 2c 20 61 6e 64 20 4c 2e |. Keller|, and L.|
|00005c80| 7e 52 68 65 69 6e 68 61 | 72 74 7d 2c 20 7b 5c 65 |~Rheinha|rt}, {\e|
|00005c90| 6d 20 43 6f 6e 74 69 6e | 75 61 74 69 6f 6e 2d 63 |m Contin|uation-c|
|00005ca0| 6f 6e 6a 75 67 61 74 65 | 0a 20 20 67 72 61 64 69 |onjugate|. gradi|
|00005cb0| 65 6e 74 20 6d 65 74 68 | 6f 64 73 20 66 6f 72 20 |ent meth|ods for |
|00005cc0| 74 68 65 20 6c 65 61 73 | 74 2d 73 71 75 61 72 65 |the leas|t-square|
|00005cd0| 73 20 73 6f 6c 75 74 69 | 6f 6e 20 6f 66 20 6e 6f |s soluti|on of no|
|00005ce0| 6e 6c 69 6e 65 61 72 20 | 62 6f 75 6e 64 61 72 79 |nlinear |boundary|
|00005cf0| 20 76 61 6c 75 65 0a 20 | 20 70 72 6f 62 6c 65 6d | value. | problem|
|00005d00| 73 7d 2c 20 53 49 41 4d | 20 4a 2e 20 53 63 69 2e |s}, SIAM| J. Sci.|
|00005d10| 20 53 74 61 74 69 73 74 | 2e 20 43 6f 6d 70 75 74 | Statist|. Comput|
|00005d20| 2e 2c 20 36 20 28 31 39 | 38 35 29 2c 20 70 70 2e |., 6 (19|85), pp.|
|00005d30| 7e 37 39 33 2d 2d 38 33 | 32 2e 0a 0a 5c 62 69 62 |~793--83|2...\bib|
|00005d40| 69 74 65 6d 7b 47 6f 6c | 75 62 2d 56 61 6e 4c 6f |item{Gol|ub-VanLo|
|00005d50| 61 6e 7d 0a 7b 5c 73 63 | 20 47 2e 7e 48 2e 20 47 |an}.{\sc| G.~H. G|
|00005d60| 6f 6c 75 62 20 61 6e 64 | 20 43 2e 7e 46 2e 20 56 |olub and| C.~F. V|
|00005d70| 61 6e 7e 4c 6f 61 6e 7d | 2c 20 7b 5c 65 6d 20 4d |an~Loan}|, {\em M|
|00005d80| 61 74 72 69 78 20 43 6f | 6d 70 75 74 61 74 69 6f |atrix Co|mputatio|
|00005d90| 6e 73 7d 2c 20 0a 20 20 | 53 65 63 6f 6e 64 20 65 |ns}, . |Second e|
|00005da0| 64 2e 2c 20 54 68 65 20 | 4a 6f 68 6e 73 20 20 48 |d., The |Johns H|
|00005db0| 6f 70 6b 69 6e 73 20 55 | 6e 69 76 65 72 73 69 74 |opkins U|niversit|
|00005dc0| 79 20 50 72 65 73 73 2c | 20 42 61 6c 74 69 6d 6f |y Press,| Baltimo|
|00005dd0| 72 65 2c 20 4d 44 2c 20 | 20 31 39 38 39 2e 0a 0a |re, MD, | 1989...|
|00005de0| 5c 62 69 62 69 74 65 6d | 7b 4d 6f 72 65 7d 0a 7b |\bibitem|{More}.{|
|00005df0| 5c 73 63 20 4a 2e 7e 4a | 2e 20 4d 6f 72 5c 27 65 |\sc J.~J|. Mor\'e|
|00005e00| 7d 2c 20 7b 5c 65 6d 20 | 41 20 63 6f 6c 6c 65 63 |}, {\em |A collec|
|00005e10| 74 69 6f 6e 20 6f 66 20 | 6e 6f 6e 6c 69 6e 65 61 |tion of |nonlinea|
|00005e20| 72 20 6d 6f 64 65 6c 20 | 70 72 6f 62 6c 65 6d 73 |r model |problems|
|00005e30| 7d 2c 20 69 6e 0a 20 20 | 43 6f 6d 70 75 74 61 74 |}, in. |Computat|
|00005e40| 69 6f 6e 61 6c 20 53 6f | 6c 75 74 69 6f 6e 73 20 |ional So|lutions |
|00005e50| 6f 66 20 4e 6f 6e 6c 69 | 6e 65 61 72 20 53 79 73 |of Nonli|near Sys|
|00005e60| 74 65 6d 73 20 6f 66 20 | 45 71 75 61 74 69 6f 6e |tems of |Equation|
|00005e70| 73 2c 20 45 2e 7e 4c 2e | 20 41 6c 6c 67 6f 77 65 |s, E.~L.| Allgowe|
|00005e80| 72 20 61 6e 64 0a 20 20 | 4b 2e 7e 47 65 6f 72 67 |r and. |K.~Georg|
|00005e90| 2c 20 65 64 73 2e 2c 20 | 4c 65 63 74 75 72 65 73 |, eds., |Lectures|
|00005ea0| 20 69 6e 20 41 70 70 6c | 69 65 64 20 4d 61 74 68 | in Appl|ied Math|
|00005eb0| 65 6d 61 74 69 63 73 2c | 20 56 6f 6c 2e 20 32 36 |ematics,| Vol. 26|
|00005ec0| 2c 20 41 6d 65 72 69 63 | 61 6e 0a 20 20 4d 61 74 |, Americ|an. Mat|
|00005ed0| 68 65 6d 61 74 69 63 61 | 6c 20 53 6f 63 69 65 74 |hematica|l Societ|
|00005ee0| 79 2c 20 50 72 6f 76 69 | 64 65 6e 63 65 2c 20 52 |y, Provi|dence, R|
|00005ef0| 49 2c 20 31 39 39 30 2c | 20 70 70 2e 7e 37 32 33 |I, 1990,| pp.~723|
|00005f00| 2d 2d 37 36 32 2e 0a 0a | 5c 62 69 62 69 74 65 6d |--762...|\bibitem|
|00005f10| 7b 53 61 61 64 7d 0a 7b | 5c 73 63 20 59 2e 7e 53 |{Saad}.{|\sc Y.~S|
|00005f20| 61 61 64 7d 2c 20 7b 5c | 65 6d 20 4b 72 79 6c 6f |aad}, {\|em Krylo|
|00005f30| 76 20 73 75 62 73 70 61 | 63 65 20 6d 65 74 68 6f |v subspa|ce metho|
|00005f40| 64 73 20 66 6f 72 20 73 | 6f 6c 76 69 6e 67 20 6c |ds for s|olving l|
|00005f50| 61 72 67 65 20 75 6e 73 | 79 6d 6d 65 74 72 69 63 |arge uns|ymmetric|
|00005f60| 0a 20 20 6c 69 6e 65 61 | 72 20 73 79 73 74 65 6d |. linea|r system|
|00005f70| 73 7d 2c 20 4d 61 74 68 | 2e 20 43 6f 6d 70 2e 2c |s}, Math|. Comp.,|
|00005f80| 20 33 37 20 28 31 39 38 | 31 29 2c 20 70 70 2e 7e | 37 (198|1), pp.~|
|00005f90| 31 30 35 2d 2d 31 32 36 | 2e 0a 0a 5c 62 69 62 69 |105--126|...\bibi|
|00005fa0| 74 65 6d 7b 53 61 61 64 | 2d 53 63 68 75 6c 74 7a |tem{Saad|-Schultz|
|00005fb0| 7d 0a 7b 5c 73 63 20 59 | 2e 7e 53 61 61 64 20 61 |}.{\sc Y|.~Saad a|
|00005fc0| 6e 64 20 4d 2e 7e 48 2e | 20 53 63 68 75 6c 74 7a |nd M.~H.| Schultz|
|00005fd0| 7d 2c 20 7b 5c 65 6d 20 | 7b 5c 72 6d 20 47 4d 52 |}, {\em |{\rm GMR|
|00005fe0| 45 53 7d 3a 20 41 20 67 | 65 6e 65 72 61 6c 69 7a |ES}: A g|eneraliz|
|00005ff0| 65 64 20 6d 69 6e 69 6d | 61 6c 0a 20 20 72 65 73 |ed minim|al. res|
|00006000| 69 64 75 61 6c 20 6d 65 | 74 68 6f 64 20 66 6f 72 |idual me|thod for|
|00006010| 20 73 6f 6c 76 69 6e 67 | 20 6e 6f 6e 73 79 6d 6d | solving| nonsymm|
|00006020| 65 74 72 69 63 20 6c 69 | 6e 65 61 72 20 73 79 73 |etric li|near sys|
|00006030| 74 65 6d 73 7d 2c 20 53 | 49 41 4d 20 4a 2e 20 53 |tems}, S|IAM J. S|
|00006040| 63 69 2e 20 53 74 61 74 | 69 73 74 2e 0a 20 20 43 |ci. Stat|ist.. C|
|00006050| 6f 6d 70 75 74 2e 2c 20 | 37 20 28 31 39 38 36 29 |omput., |7 (1986)|
|00006060| 2c 20 70 70 2e 7e 38 35 | 36 2d 2d 38 36 39 2e 0a |, pp.~85|6--869..|
|00006070| 0a 5c 62 69 62 69 74 65 | 6d 7b 53 77 61 72 7a 74 |.\bibite|m{Swarzt|
|00006080| 72 61 75 62 65 72 2d 53 | 77 65 65 74 7d 0a 7b 5c |rauber-S|weet}.{\|
|00006090| 73 63 20 50 2e 7e 4e 2e | 20 53 77 61 72 7a 74 72 |sc P.~N.| Swarztr|
|000060a0| 61 75 62 65 72 20 61 6e | 64 20 52 2e 7e 41 2e 20 |auber an|d R.~A. |
|000060b0| 53 77 65 65 74 7d 2c 20 | 7b 5c 65 6d 20 45 66 66 |Sweet}, |{\em Eff|
|000060c0| 69 63 69 65 6e 74 20 7b | 5c 72 6d 20 46 4f 52 54 |icient {|\rm FORT|
|000060d0| 52 41 4e 7d 0a 20 20 73 | 75 62 70 72 6f 67 72 61 |RAN}. s|ubprogra|
|000060e0| 6d 73 20 66 6f 72 20 74 | 68 65 20 73 6f 6c 75 74 |ms for t|he solut|
|000060f0| 69 6f 6e 20 6f 66 20 65 | 6c 6c 69 70 74 69 63 20 |ion of e|lliptic |
|00006100| 70 61 72 74 69 61 6c 20 | 64 69 66 66 65 72 65 6e |partial |differen|
|00006110| 74 69 61 6c 20 65 71 75 | 61 74 69 6f 6e 73 7d 2c |tial equ|ations},|
|00006120| 20 41 43 4d 0a 20 20 54 | 72 61 6e 73 2e 20 4d 61 | ACM. T|rans. Ma|
|00006130| 74 68 2e 20 53 6f 66 74 | 77 61 72 65 2c 20 35 20 |th. Soft|ware, 5 |
|00006140| 28 31 39 37 39 29 2c 20 | 70 70 2e 7e 33 35 32 2d |(1979), |pp.~352-|
|00006150| 2d 33 36 34 2e 0a 0a 5c | 62 69 62 69 74 65 6d 7b |-364...\|bibitem{|
|00006160| 57 61 6c 6b 65 72 38 38 | 7d 0a 7b 5c 73 63 20 48 |Walker88|}.{\sc H|
|00006170| 2e 7e 46 2e 20 57 61 6c | 6b 65 72 7d 2c 20 7b 5c |.~F. Wal|ker}, {\|
|00006180| 65 6d 20 49 6d 70 6c 65 | 6d 65 6e 74 61 74 69 6f |em Imple|mentatio|
|00006190| 6e 20 6f 66 20 74 68 65 | 20 7b 5c 72 6d 20 47 4d |n of the| {\rm GM|
|000061a0| 52 45 53 7d 20 6d 65 74 | 68 6f 64 20 75 73 69 6e |RES} met|hod usin|
|000061b0| 67 0a 20 20 7b 48 7d 6f | 75 73 65 68 6f 6c 64 65 |g. {H}o|useholde|
|000061c0| 72 20 74 72 61 6e 73 66 | 6f 72 6d 61 74 69 6f 6e |r transf|ormation|
|000061d0| 73 7d 2c 20 53 49 41 4d | 20 4a 2e 20 53 63 69 2e |s}, SIAM| J. Sci.|
|000061e0| 20 53 74 61 74 69 73 74 | 2e 20 43 6f 6d 70 75 74 | Statist|. Comput|
|000061f0| 2e 2c 20 39 20 28 31 39 | 38 38 29 2c 0a 20 20 70 |., 9 (19|88),. p|
|00006200| 70 2e 7e 31 35 32 2d 2d | 31 36 33 2e 0a 0a 5c 62 |p.~152--|163...\b|
|00006210| 69 62 69 74 65 6d 7b 57 | 61 6c 6b 65 72 38 39 7d |ibitem{W|alker89}|
|00006220| 0a 5c 73 61 6d 65 61 75 | 74 68 6f 72 2c 20 7b 5c |.\sameau|thor, {\|
|00006230| 65 6d 20 49 6d 70 6c 65 | 6d 65 6e 74 61 74 69 6f |em Imple|mentatio|
|00006240| 6e 73 20 6f 66 0a 20 20 | 74 68 65 20 7b 5c 72 6d |ns of. |the {\rm|
|00006250| 20 47 4d 52 45 53 7d 20 | 6d 65 74 68 6f 64 7d 2c | GMRES} |method},|
|00006260| 20 43 6f 6d 70 75 74 65 | 72 20 50 68 79 73 2e 20 | Compute|r Phys. |
|00006270| 43 6f 6d 6d 2e 2c 20 35 | 33 20 28 31 39 38 39 29 |Comm., 5|3 (1989)|
|00006280| 2c 20 20 70 70 2e 7e 33 | 31 31 2d 2d 33 32 30 2e |, pp.~3|11--320.|
|00006290| 0a 0a 5c 65 6e 64 7b 74 | 68 65 62 69 62 6c 69 6f |..\end{t|hebiblio|
|000062a0| 67 72 61 70 68 79 7d 20 | 0a 0a 5c 65 6e 64 7b 64 |graphy} |..\end{d|
|000062b0| 6f 63 75 6d 65 6e 74 7d | 20 0a 0a |ocument}| .. |
+--------+-------------------------+-------------------------+--------+--------+